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Equation (3) on page 858 should read as follows: 


x= &(p)ixy + We P)Xiy tooet WrikPXik tees l. 


On the same page, all the x's in Fig. 1 should be lower case (xo, Xqs Xoo» and x) rather than capitals, 

















AN ESTIMATE OF THE MAXIMUM DEVIATION 


FROM A GIVEN TRAJECTORY 
E. A. Barbashin 


Sverdlovsk 


Translated from Avtormatika i Telemekhanika, Vol. 21, No. 10, 


pp. 1341-1351, October, 1960 
Original art icte submitted April 15, 1960 


The method proposed in [1] is used to estimate the maximum deviation of motion from a given 


trajcctory. 


A method of choosing control functions and vectors is introduced for decreasing this devia- 
SY tion. The case of an infinite time interval is considered. The relation between the problem 
Pa investigated and that of the accumulation of disturbances is studied. 


We consider the system of differential equations 


dx 2 -- 
ae = Ste) eat+ SY burr (d) 
k=1 k=1 
(1) 


(i=1,2,... 2, m<a), 


where aj, (t) are continuous functions on the interval 
to St St, + T, and the bj, are constants such that the 
“control” vectors by (bik, bzk,..-» bak) (kK = 1, 2, ..., m) 
are linearly independent. 

As in [1], we pose the problem of finding control 
functions u,(t), such that for te St =t, + T the relation 


xj (t) = fi (t) (i = 1, 2,..., 2) (2) 
is satisfied, where xj;(t) is the solution of the system (1) 
and f ;(t) is a given system of functions with xj(t9) = 
=f; (te). When m < n, this problem is, as a rule, insol- 
uble, and so we try to find functions u;,(t) [2] such that 
the expression 


||z|| = max max 
i 


X(t)— fill 3 
max IQ —-hOL © 
will be a minimum, 

Following the ideas of [1], we will not try to find 
control functions guaranteeing that the value of Ii zll is 


exactly a minimum, but will obtain an estimate of the 
form 


izk SAllyil, 


Where A is a constant and Ilyll is the deviation from 
zero in some sense of another quantity y. We replace 


the difficult problem of finding minll z ll by the simpler 
problem of finding min iyi, The loss in accuracy of 
the approximation due to this procedure is not great for 
small values of liyll, while the decrease in difficulty in 
the numerical calculations is considerable. 


In Section 1 we show that the method of choosing 
the control functions and vectors proposed earlier en- 
sures a decrease, not only in the mean square value of 
the deviation, but also a decrease in the maximum de- 
viation, It must be taken into consideration that the 
value of a control function at any given instant is chosen 
independently of the length of the interval over which 
the approximation is being carried out, and also indepen- 
dently of the value of this function at any other instant 
of time. 


In Section 2, the same method of choosing the con- 
trol functions [see the system (10)) is used to obtain 
other methods of estimating the maximum deviation. 
Although these other methods have no special advan- 
tages over the estimate of Theorem 1, they are of in- 
terest in several cases. In particular, Theorem 3 (and 
also Theorems 1 and 2) solves the problem of the accum- 
ulation of disturbances in 2 dynamic system. From this 
point of view, the results obtained are completely anal- 
agous to those obtained by B. V. Bulgakov in (3, 4]. We 
also note that the character of any of the three estimates 
of lizii depends on the function space in which we con- 
sider y. An explanation of this phenomenon is given in 
the Appendix. 


We consider the case T = @ in Section 3, It is 
natural to consider this case in connection with questions 
in the theory of stability. It is unfortunately impossible 
in the present paper to investigate this problem in the 
light of the ideas in the detailed investigation by Mas- 
sera and Sch®8ffer (5). 
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1. The Basic Estimate of the Maximum 
of the Deviation and the Choice of Con- 








trol Functions 
In the system (1) we make the change of variables 





a= %— fil(t) (i = 1, 2,..-, 7), 


and obtain the new system 


n 


“a4= > Aix (t)Z_ + ui (t) (i =1,2,..., 2), (4) 
k=1 


where 


& n. . 
yi(t) = >) Outs (t) + >) aie OP ) — A) 
k=1 k=1 
(i = 1, 2,..., ”). 
(5) 
The solution of the system (4) that satisfies the con- 
ditions 2j(t,) = 0 is 


ton 
z(t) = \> wi (t,t) yx (t) dt (i= 1,2,...,m), (6) 


t. k=1 


where Wj,(t, 7) are the elements of the fundamental 
matrix W(t,T) defined by the equation dW/ dt = A(t)W 
and the condition W(t,T ) = I [I is the unit matrix, and 
A(t) is the matrix with its elements, coefficients ajk 
(t)]}. 

In what follows we will use the notation 


||z|| = max max {z(t)|, 
i tet<t,47 


y*(t)= > y? (0), 


i=1 


w;? (t,t) = > wi; (t, T) (i = 1, 2,.00 2). 
k=1 


Let yj(t) be functions which are for the moment 


te+T 
arbitrary, but such that the integral 
lo 
The following theorem opens up the possibility of 
estimating lizll in terms of the mean square value of 
the quantity y(t). 
Theorem 1. There exists a minimum number M, 
for which 


t+T 


zli<@( \ var) 


'/, 


(7) 


y? (t) dt-exists. 


The number M can be obtained from the formula 


t t/, 
M=m 2(t, t)dt). 
techie T (\ (¢, *) ) (8) 


If the functions a;},(t) are constants, then 


t,+T 
M = max ( wi (to + T, x) dx ) 


te 


"I, 


(9) 


Theorem 1 evidently also solves the problem of the 
accumulation of disturbances when the mean square of 
the external forces is bounded. In fact, if 


te+T 
\ v@)ar=Z, 
to 
then the formula (7) gives the maximum deviation, 
equal to ML. From the proof of Theorem 1 given in the 
Appendix. it is obvious what form the functions yj (t) 


must take to give the maximum deviation. 
It follows from Theorem 1 that ll zll can be de- 


creased by decreasing the quantity ||y|| = 


- #3 y? (x) dt ' 
& 


The method of choosing the control functions uj (t) 
to ensure that lizli is a minimum was given in Section 


2 of [1]. In this paper it was shown that the control func- 


tions are determined by a system of equations of the 
form 


™m 

S} (Bx, bi) wi (2) = (7), by) (= 1, 2,-0-5 my (10) 

i=1 
where(6j,, bj) is the scalar product of the vectors by; and 
bj, andr (t)= 1, (t), r2(t),..., In (t) are vector functions 
satisfying the equation 


n 


nif@=A— by Aix (t) fx (2) 


k=1 


(i = 1, 2,..., m). 
(11) 


If uj(t) is chosen according to (10), then we will 
have 


. I (by, be,---, Dy» ¥ (2) 
y* (t) orig b ye (2) — Ths, Bac-++s De) » (12) 


i=1 





where the numerator and denominator in the fraction on 
the right-hand side are the Gramme determinants corres- 
ponding to the systems of vectors. 








a 


tan 


Se ah & 


1 
i 





Thus, we have 


H? = |ly|? = 


‘ tet? 
=T (bb; b,) ) I’ (b;, bg,..., Dm, 7 (t)) at. 





(13) 
It should be noted that if the functions u,(t) are 
chosen according to (10), then 


>) eee (t) = 
k=1 


(b,, b,)... (b,, Dm), b, 
= Fi bab) | (Pm: bi)... (ms Bm), Bm 
(r, b,)...(f, Dm), 0 





In Section 3 of [1], a recommendation is given for 
the choice of the control vectors by, bz, ..., bm, to give 
a minimum value for il yll if the functions uz, (t) are 
chosen according to (10). The choice of the control 
functions and the control vectors given in [1] thus con- 
tributes,not only to decreasing the mean square devia- 
tion, but also to decreasing the maximum deviation. 


2. Other Estimates of iizil 
In this section we give other estimates of li zil which 
are independent of the estimate given inSection 1. Let 
yi(t) be an arbitrary function such that the integral 
tT 
\ y (t) dt exists. 
to 
theorem 2. The inequality 





te+T 
lz\|<N \ y(eyae 


t. 


(14) 


holds, where N = max max max wj(t,T). 

i tpStStg+Ty to STSt 

If all the eigenvalues Aj (t) of the matrix A(t) + 
+ A%t) (A(t) is the matrix having the coefficients a;}(t) 
as elements, A *(t) is its adjoint) are nonpositive for any 
value of t in the interval ty) St St + T, then N = 1, 
and 
tefT 
Wii \ y(eae. 


Formula (14) gives a new estimate of iizil. In 
connection with the choice of control functions, it fol- 
lows from Section 2 of [1] that Theorem 2 gives no new 
information in choosing these functions, except that 
which follows from the system (10) investigated in the 
present paper. In fact, when th. control functions are 
chosen according to (10), the minimum value of not 
t+T 

\ y*(t)dt, but also that of the 
t, 


only the .ntegral 


tT 
integral \ y(t)dt is ensured. 
t, 


Theorem 3, There exists a least number C, such 
that 


lz} SC max y(?). 


(15) 
teSt<ttT 


This number is given by the formula 


t 


C=max max { wy (t, t)drt. 


(16) 
t tyct<t,4+T Fy 


If the functions aj, (t) are constant, then 


to+T 


C = max wi(to-+- 7, t)dt. (17) 
i 


t 
In addition to the basic problem, Theorem 3 also 
solves the problem concerning the accumulation of dis- 


turbances in the system for bounded external vector for- 
ces. If the functions yj (t) satisfy the inequality 1/* (t) = 


n 
_ > yi (th QR, 
i=1 
ula (15) that li ziiSRC, while,from the proof of Theorem 


3 given in the Appendix, it follows that the upper bound 
of RC is attained for 


then it follows from Form- 


where t and i are chosen so that the maximum (16) is 
attained. If the functions aj},(t) are constant, then we 
have 
tT 
| z || <— R max \ wi (to + 7, 1) dt. 
i , 
te 
The result of the solution in the present paper and 
in [1] differ from that in [3] in that here the sum of the 
squares of the external influences yj (t) is bounded, and 
not each separate influence, 
We now return to the main problem. Since y;(t) = 
= >) dinux (t) — ri (O), where the functions rj (t) are 
k=1 
given by Formula (11), then according to [1) y(t) = 


n 
- (> ye (t))'* gives the distance from a point A lying 
i=1 


on the curve zj =r; (t) to a point B in the plane Q gener- 
ated by the vectors b,, by, . . . . By. At any given in- 
stant t, the distance AB will be a minimum if the func- 
tions u; (t) are chosen according to (10), In this case 

p (t)= y*(t)= AB® are calculated from Formula (12),while, 
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a 


for the indicated choice of the control functions the 
quantity p(t) will be a minimum. In order to use Theo- 
rem 3, we must calculate the greatest value of the 
function p (t) for ta<t<tg+T. This can obviously be 
done by elementary methods. In fact, since p(t)=-[r(t), 
y(t)], where g(t) is the vector with projections r;(t) and 
y(t) is the vector with projections y,(t), we have p(t} 
—(r,y)-(y. 2), 

But, according to [6] (p. 205), 


(b,, b,), it il (b,, Dm), b, | 


0-4 “: a, “ih sos be 6 oad: 


1 
¥()= =F 1h, Bm), --. 5 (Dy Bin), Dm 





where [ =I'(by,b, . . . , bry) is Gramme determinant 
for the vectors by,b;,... . Dm. 
From (18) we see that 


(y, r) site 


(b,, by), ..., (by, by), (by, vr) 
pap snoody a Get ade De es 
=T | (b,, bn), --., (Dm, Bm), (bm, vr) |= @ ¥)- 

(b,, r), ..-, (bm, r), (Fr, ¥) 


A necessary condition for a maximum is therefore 


(r, y) = 0; (19) 


and this condition is satisfied if <0, or if # is orthogo- 
nal toy. We will not consider tie case y=0, since in 


n 
this case y? = > yi? = 0 , and there is consequently 
i=) 
no maximum. 
To obtain a necessary condition, we introduce the 
parameter s, denoting the arc length of the curve z; = 
=r; (t). We have 


5 =—(@, y)—(, ¥) = 


=— (S455, y)—P(b, by... bm, ¥). 


The equality(, y) = TQ, b., sees bm: t) is ob- 
tained directly from (18). 


When we take into consideration that 4 "y} = 0, 


we can write the sufficient condition for a maximum in 
the form 
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’ bm, r) > 0. (20) 


(Sx. y)s*+T(b, by, ... 


In order to find the greatest value of p (t),we must 
calculate p(t) at the points t, and t+ T,*and also at the 
maximum points, and choose the greatest value of p (t) 
from these values, 


3. The Approximate Trajectory in the 





Case T= 

We now consider the possibility of approximating 
to the given trajectory for an infinite time interval. We 
first of all set an upper limit to the absolute value of 
the elements of the fundamental matrix, and assume 
((7), p. 366) that 





| wi, (t, t)| << Be#t—*) 
(21) 
» BR, kami, 2, ..., MR). 


(i= 4, 2,... 


This condition is equivalent to assuming that any 
solution of the system 


n 
x; = > Aix, (t) Ly ure, go asc > (22) 
k=1 
satisfies the inequality 


| 21 (t)|< B( 2s ay? (t0))* ema, 


k=1 


In fact, an arbitrary kth column wy, (t, T), Wek(t.T), 


. ++» Wok(t,7) of the matrix W(t,T) satisfies the system 
(22) if wy, (7, T)=0 for i#k and wi (7, 7)=1. 

The condition (21) is always satisfied if the matrix 
A(t)=||a;,, [| is constant and its eigenvalues have nega- 
tive real parts. 

It follows from Theorems 1, 2, and 3, that when 
the conditions (21) are satisfied, 


co 
n 


Iz <B (= 


/, 
se \ v*(n)dr)", (23) 


* 
o 


2 || <BV n \ y(x)dr, (24) 


og 


Wzi<B* sup yi). (25) 
t «<t<oo 
It is obvious that in order that the solution of the 
system (4) be bounded, it is sufficient that one of the 
conditions co 
\ y? (t) dt < co, (26) 








-_> -—+ | -—« 


oD ~w 


Law | 








00 5. The proposed method is most effective for small 





( y(t) dt < ©, (27) min||y||, and when the conditions (21) are satisfied. 
" APPENDIX 
Proof of Theorem 1. From the relation (6) we have 
sup y(t) < co (28) 
ty<t<co t 
be satisfied. ja |< \ w, (t, T) y (t) dt, (29) 
I. G. Malkin showed in [8] that it follows from the i, 


boundedness of all solutions of the system (2); for arbi- . ‘ z 
trary yj(t) satisfying condition (28), that condition (21) is and w on we use the Bunyakovskii- Schwartz inequality, 
satisfied, we obtain 


We now consider what conditions must be imposed : P 
on the functions fo) in order that the inequalities (26), : : ‘at. s 
(27), and (28) be satisfied when the control functions are I4OIS (\w (¢, 7) ar) ({ # (t) et ) Ss 
chosen according to (10). : > 

Since from (12) we have 





t te+T 
“ls ‘Ys 
S 2(t, t) d 2 (t)dt), 
een re <Q erener)"(] rena) 
y* (i) = ik... ay (30) 
This inequality becomes an equality when the func- 


these conditions must obviously have the form tion y, (rT) has the form 


t pn es! f O<t<t 
| T(by, by, ..., Bs F(Q))AE< 00, (26% nO= ea fo , 


: y,(t)=0 for tI t<te+T. 


oo 





\ VI (b,, iio.zatsI, Din, r (t)) dt <0, (27') From (30) we obtain 
te t ‘ t+T ; 
sup I'(b,, ba, ..., Bm, r(é))<90. (28%) max | z;(¢)|< max (| w,* (t, T) at) ‘( j y* (tT) dt) ; 
t<t<oo te<t<te+T te<t<i,+T i. i, 
Here the vector function f(t) is given by formula from which the formulas (7) and (8) follow, 
(11). It is also easy to show that the choice of the sys- Now let aj, (t) be constant; we will show that 
tem of control vectors according to [1] will be possible : 
ro ) F (t, toe) = Jue, t)dt is a monotonically increasing 
if the condition ( r? (t) dt < ©° is satisfied. te 
if function of t. Actually, in this case ([9], p. 173) we have 


Wik =(t, T)=Wi,(t=r). If we let tg<t,<t,, we can write 
Conclusions 


1. The choice of control functions according to (10) ty 
makes it possible to decrease not only the mean square F (ts, t) = \w? (tp —t) dt= 
deviation, but also the maximum deviation. a 

2. The values at any given instant of time of con- 





trol functions chosen according to (10) do not depend on tette—ts te 
the length of the time interval over which the approxi- = w,* (tg— T) dt + \ w,* (ta—-T) dt. 
mation is being made, or on the values of the control te tole 
functions at any other instant of time. 
3. The choice of control functions according to (10) After the change of variable n =T~t,+t,, in the 
yields the possibility of selecting an optimal system of second integral we have 
control vectors by using the method described in Section tettsa—ty ty 
3 of [1}. F (te, te) = wi (— tae + | w2(s—n) dn = 
4. The proposed method makes it possible to avoid i, i, 
difficulties in numerical calculation, and permits asim- totte—ts 
ple geometric interpretation and the application of mod- = w,? (te— T) at + F (ts, to). 
elling for the solution of the problem. te 
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ey + eee 


— 





Since the first integral is positive, we obtain the de- 
sired result. 
Formula (9) follows from the fact that F(t, t),as a 
monotonically increasing function, attains its maximum 
for t=to+ T. 
We note that in the case of variable aj;,(t), the 
conclusion that F(t, ts) is increasing is no longer true. 
Similarly, consider the equation z=-tantz. The gen- 
eral solution of this equation is z = c cost, w(t,T ) = 
t 

cos t/cos T, F(t, 0) = [wee rier = 4 sin 2t. The func- 
0 

tion F(t, 0) is evidently not monotonic in the interval 


o ot es t< e+e if € is sufficiently small. 


Proof of Theorem 2. We conclude from the inequal- 
ity (29) that 





t 
}2;(t)| < max w,(t, T) \yimer; 
octet i, 


when we increase the upper limit of integration to t9+T, 
we obtain 


te+T 
, d , 
Ini s mee 4S T) i y(t) dt 


and (14) follows. 
We now consider the case when the eigenvalues of 
the matrix A(t) + A®%({t) are nonpositive for allt=t, In 


this case, the quadratic form >) a,, (t) 2,7, will also be 
‘k= 

nonpositive, It is known ([9], p. 172) that the matrix 

W(t, T) is equal to F(t) Fl (r), where F(t) is the funda- 

mental matrix of the system x= A(t)x. It is thus easy to 

obtain the relation 


dw 


from which it follows that the rows(wWj;, Wie... - 
the matrix W(t,T) satisfy the system 


» Win) of 


n 


dw 
z= — 2 tm (T) wi (=A, Qu... 0). 


n 
Recalling that w,*(t, t) = >} w;,? (t, T), we have 
k=1 


n n 
5 s(t, 1) =—2)) Wi, b Onyx (T) Yim = 


k=1 m=1 


=—2 > Ong (T) %,%im = 9, 


k,m=1 
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according to the conditions of the theorem. Thus, 


max w,*(t, tT) = w,2(t, t) = 1. It should be noted that 
<< 


the condition that the eigenvalues of the matrix 
A(t) + A® (t) be nonpositive can be replaced by the 
Sylvester criterion, or by a modification of this criterion 
([6], p. 250). 

Proof of Theorem 3. It follows from (29) that 





ln(1< ( w, (t, T)dt max y(t), 
i, ti<t<t 


and so 
1 


max |2,(t)|<— max Jc, ae max y(t), 


te<t<t,4+T to<t<te+T ‘ eSt<te+-T 
(31) 
which yields 
|zi|< Cmax y(t), 
te<t<i,+T 
where 


t 


|» (t, t) dt. 
i, 


C=muax max 
i t<t<t.4T 


The fact that C is the least of the numbers satisfy- 
ing the inequality follows, because the right-hand side of 
the inequality (31) is attained for 


w,; (t, T) 
y; (t) =D) * 


Reasoning similar to that used in the proof of 
Theorem 1 can also be used to prove the validity of 
Formula (17). 

Note, We will consider the relation (6) as the appli- 
cation of an operator W, which transforms the vector 
function y(t) into the vector function z(t), In the pres- 
ent paper z(t) is an element of the space C, (t», tg+ T), 
in which the norm is defined by the formula |; z || = 
=P sig Ol - 

For each of the three theorems, the function y(t) 
lies in a different space. 

In the first theorem, y(t) lies in the space 


te+T 
' 
Len (to, to+ T) with the norm ll y Il =( t y (0) at ) A 


te 
In the second theorem, y(t) lies in the space Lyy(to, tg+T) 


te T 
with norm || y|| = t y (t)dt , and in the third theorem 
i, 
y(t) lies in the space Cp," (to, tg+ T) with norm 


yll= max y(t). 
te<t<te+T 











In the first and third theorems we calculate the 
norm of the operator W, i.e., the least value of the num- 
ber A for which the inequality ||z||<Al]y|| is valid. 

We note that the problem concerning the accumu- 
lation of disturbances in an arbitrary trajectory is essen- 
tially the problem of determining the norm of the opera- 
tor W. 


Example, We consider the system of differential 
equations 


4y= — 21+ du (t) 


(by? + b,* = 1). (32 
%2.=> — 322 — bau (t) 
The fundamental matrix W(t, T) of the system x,= 


= —X;, X= —3x, has the form 


e~'—*), 9 
W (t, Tt) = 








0 : e 3tt—*) 


Let the curve x,=e~™*, x,=e"* be given, By 


choosing u(t), we will try to find an approximation to 
this curve in the interval 0<t=T. It follows from (11) 
that r,(t)= —e"™*, r,(t)= e"**. The system (10) yields 
u(t)=e"™" (b,— b,). By using the formulas in (12) we 
obtain 


y? (t) = e—* [2 — (bg — by)?), 


and from (13), 


2 — (b2— bi)? 
4 


T 
H? = \ y*(t) dt = (i —e*?). 
0 


In order to use Theorem 1, we must calculate the 
number M from (9). We have 


7 
1 
M}*= ( e~AT—*) dy = t= e*T), 
0 


T 
1 
M;? = \ eT) ay = = (1 —e*7), 
0 


and se 


M = max vz (1 —e~*T), Vi (i — e~*T) ). 


We thus have || z|| =< MH, and for T= © we have 


yey 2 ba — bi) 
H 2V2 . 


If we now carry out the analysis according to Sec- 
tion 3 of [1], or if we investigate directly the result we 
have obtained, we can verify that for b,= 2/2, and by= 
= —/5/2, we can obtain an accurate representation of 
the given trajectory. 

We now consider the method of estimation of Sec- 
tion 2 of the present paper, According to Theorem 2 we 
have N=1, 

On the other hand, 

T 


H,= ( y (t)dt= 


0 





oo" 


siinicngilliaial 2— (,— bi? 
= \ eo *“Y2—(@.—hypdi= V Ms a) (4 —e27), 
From Theorem 2 we obtain 
2—(.— bi)? 
2 <n = VOOM _ematy 


and from T= 


V2 — (b:— bi)* 
ll - \\ S TF ae 


In order to apply the result of the third theorem, we 
calculate 


owe 


T 
C= max eT) 4, \ e~T—*) ax} = 
0 


--3T 


1 
=max a—e", zi —e )}- 


Moreover, 
max y(t) = V2 — (b,—6))?. 
0<t<T 
Thus, according to (15), we obtain 


| 21S CV 2—(b:— by)’ 


and for T= 
Wei < V2—@.— bi). 


As we can see, for T= the best estimate of ||z|| is 
obtained from Theorem 1, and the worst from Theorem 3. 
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The stable modes of an airplane wing are found with an analog computer for the case where 
there are no nonlinear elements in the aileron control system. It is shown that the presence 

in the control circuit of a gap in certain parameters leads to the excitation Of hard oscillations 
in the system, The probability of exciting growing wing oscillations over a given period of 
time in which the air flow is turbulent is calculated, 


An airplane wing in flight is a complex nonuniform 
oscillatory system with distributed parameters, An ac- 
curate solution of the wing stability problem is practical- 
ly impossible, even with the help of modern computers, 
At present, only approximate methods are known for 
solving this problem to one or another degree of accura- 
cy. 

Replacement of the nonuniform distributed oscilla- 
tory system with a simplified lumped- parameter system, 
described by a finite number of ordinary differential 
equations, is the basis of all the approximate methods, 

One of the most widely used replacement methods 
is the method of [1], which is derived from the work of 
Raleigh and Pitts, and has been developed further by 
Bubnov and Galerkin. 

In the application of the Bubnov-Galerkin method 
to the simplest, but orie of the most practically impor- 
tant cases, the airplane wing, in a steady flow of air, 
can be considered as an oscillatory system with two de- 
grees of freedom, in which certa:n parameters depend 
on the flow velocity. This paper is devoted to a study 
of the oscillations and the stability of such a system 
when nonlinear elements and statistical influences 
caused by the turbulent state of the atmosphere have a 
substantial effect. 

In the majority of works in which the flutter of a 
wing with an aileron is examined (see, for example, 
[2-4]), the system is assumed to be linear, and the ai- 
leron control stick to be either fixed or free. In actual 
systems, there are always nonlinear elements as a re- 
sult of gaps, dry friction, etc. Therefore, it is quite in- 
teresting to analyze those cases in which the presence 
of these nonlinear elements cannot be overlooked be- 
cause they can have a substantial effect,not only on the 
character of the movement, but also on the stability of 
the different states of movement. An example of such 
an effect was computed in [5]. 


Nonlinear flutter problems, like all problems of 
this type, can be solved analytically in an approximate 
manner by means of the harmonic linearization method 
(Krylov-Bogolyubov method), if it is known that nearly 
harmonic modes of oscillation predominate in the sys- 
tem. The Krylov-Bogolyubov method, as it applies to 
nonlinear flutter problems, is examined, for example, 
in [7]. The use of a mathematical analog for studying 
airplane wing flutter is presented in [4, 5, 6], and oscil- 
latory systems with nonlinear hardness are examined in 
the latter two, 

The influence of turbulence on an airplane was also 
examined in a number of papers (for example, in [8)). 
In these papers, the airplane was represented as a linear 
oscillatory system, and the problem was solved by the 
ordinary theory of stationary random processes in linear 
systems, 

The presence of nonlinear elements in a system can 
lead to the appearance of periodic modes (both stable 
and unstable), with which linear theory cannot deal, 
The existence of unstable periodic modes corresponding 
to unstable limited cycles means that the excitation of 
hard oscillations in the system is possible, The latter 
case is especially interesting, both theoretically and 
practically, 

In this case, the airplane can indeed fly safely only 
in quiet ir, since dangerous oscillations of the wings 
cannot arise in the absence of disturbance (the system 
is stable in the presence of small oscillations). If ran- 
dom jolts, caused by atmospheric turbulence, act on the 
wing, dangerous growing oscillations can appear after 
some time. Therefore, the statistical evaluation of the 
probability of such oscillations arising over a given 
period of time is an important problem, whose solution 
will be presented later, 

An approximate analytic solution of this problem 
for an airplane wing with an aileron is possible under 
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the condition that the oscillations in the system are 
nearly harmonic, that is, the complex system being 
examined can be replaced with some equivalent oscil- 
latory system with one degree of freedom, whose spec- 
tral characteristic approximates that of the whole sys- 
tem in the region of the most important frequencies, 
Then the Fokker- Planck equation can be used to solve 
the problem in a way analogous to that done in [9, 10]. 
However, there is quite an amount of computation even 
with the simplifications mentioned, and the study of the 
system is exceedingly complicated over the whole range 
of the parameters v and n, where v is the flight velocity, 
and n is the relative frequency of the rotary aileron os- 
cillations. Therefore, the limits of the modes interest- 
ing to us were first determined by an analog method, 
and then the amplitudes of the limited cycles and the 
excitation probabilities were calculated for the values 
of the parameters obtained. 


1. Basic Equations of the Bending Oscil- 





lations of a Wing with an Aileron 





The Bubnov-Galerkin method gives the following 
equations for bending-aileron oscillations of a canti- 
lever wing (Fig. 1): 


Mh + 853 + Mojh = Fi (h, h, 8, 8, v) + Fe(v, t), 


Soh + 133 +- 1,38 = WM, (h, h, 3, 3, v) + M, (x, t), 
(1) 


where M, Ig, Sg, and wg are the average parameters 
of the wing and aileron, respectively, characterizing the 
wing mass, moment of inertia and static moment of the 
aileron, the natural frequency of wing bending oscilla- 
tions, and the natural frequency of aileron rotary oscil- 
lations, determined by the elasticity of the control levers 
and cables; F, and M, are the average aerodynamic force 
and moment acting on the wing and aileron which are 
moving in air with a velocity v; F, and Mp are the ran- 
dom aerodynamic force and moment caused by the tur- 
bulent disturbances of the vertical component of flight 
velocity, 





Fig. 1. Equivalent diagram of the 
bending oscillations of a wing with 
an aileron. 
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According to the theory based on the hypothesis of 
stationariness (see [11]), these aerodynamic forces and 
moments are expressed in the following manner; 


OITA ea, AR 
Fy = ap L{—F 8a bh BB, 


F, = 2npLbvu, 
“2 a 2 
M,= np L,{— 857 —t 5 b—e pRB} 


M, = npb* L, dvu, 


where pp is the air density at the given altitude; L, the 
wing span; Lg, the aileron span; b, the semichord of the 
wing; v, the horizontal flight velocity; u, the random 
vertical velocity appearing as a result of the turbulent 
state of the atmosphere; a, y, 5,6, €—coefficients 
which depend on the wing geometry. 

The numerical values of the parameters basically 
correspond to the example given in [3]: 


M = 38.3 kg; Ig = 0.0529 kgm’; Sg = 0,175 kgm ; 


@, = 30 radhee; §=96L = 4,24m; La=1,9m; 
b=0.762m; 2p =3.105 kgfm® « = 0.3458; 
v= 1.088; 8 = 0,007; £ = 0.0117; e = 0.0104. 


For convenience in analysis we will introduce the 
dimensionless coordinates, velocity and time: 


Then Eq. (1) can be written in the form 


d?x ~dzx — —dy 
ae + 2+ ay vic TF Gy, 0°Y + A307) + 


tT 


d?y d*y x eae 
+ uF Hh), qa ty + an vp? y -} 


—dy —dzx ar (2) 
+ 4y20v > eg A230 T + uaa /e(t), 
where 
Wp 1 4 
Bn il pr ly 
Op, Ai ( ) @2bM 2 @), 


/2(t) = 7) M; (=)- 


The random functions f,(r) and f,(r) can be con- 
sidered stationary and normal random processes which 
can be characterized sufficiently by giving their spec- 
tral densities, which, in their turn, are completely de- 
termined by the spectral density of the vertical compo- 
nent of turbulent air velocity u? 

* We will neglect the effect of fluctuations in the hori- 
zontal velocity, since they are significantly less than 
the flight velocity. 








A graph of the spatial spectrum of the turbulence, 
given in [8], can be used to compute the spectral den- 
sity of the random quantity uy, As emphasized in this 
work, the values of spectral density shown on the graph 
are very approximate, since they depend strongly on the 
meteorological conditions, the local terrain, and the 
flight altitude. For example, the intensity of the turbu- 
lence increases sharply during a.storm, and may differ 
from the value given on the graph by orders of magni- 
tude. 


2. Basic Equations of the Bending Oscil- 
lations of a Wing with the Aileron Con- 











trolled by a Hydraulic Amplifier 

The presence of special aileron control mechanisms 
has a substantial effect on the oscillations of the entire 
system, Let us examine the schematic diagram of a 
hydraulic amplifier aileron control as shown in Fig. 2. 
The dotted lines designate a parasitic feedback path be- 
tween the movement of the wing and system of aileron 
control levers and cables, and the hydraulic amplifier. 
Moreover, we will assume that a gap 2 can exist inthe 
piston valve control system of the hydraulic amplifier 1 
The deflection of the piston valve is composed of the 
displacement given by the pilot [the “controlling input" 
q(t)], and some additional displacement xh, due to the 
parasitic feedback loop and caused by the effect of wing 
oscillations on the aileron control drive. For simplicity, 
we will consider xh to be proportional to the vertical 
displacement of the wing h. The quantity x will be 
called the feedback coefficient. 





In the analysis, we will assume that the controlling 
input to the piston valve is constant, or that it changes 
very slowly in comparison to the oscillation period of the 
wing. Moreover, we will assume that the feedback co- 
efficient k can have only two values; +1 and ~1, 

If we limit ourselves to the linear portion of the 
servomotor characteristic, the oscillations in the system 
containing the hydraulic amplifier, for the case when 
the pilot exerts no controlling action, are described by 
the equations 


d? og - 
Sa tot anne + aye y + 


—d d* 
+ ayyv Fs + Ose Ts = f(t), 


(3) 
i2 - —d 
po + n?(y—a4z) + Og, 0 y + Ogg 5 + 


+ 0330 5 4} Ou Fa = f,(t), 3 = a, (9 (xx) — 2}. 
Here, z=r/b is the dimensionless coordinate of the 
servomotor shaft displacement 3, and the function ¢(K x) 
determines the dependence between the displacements 
of the elements on the other side of the gap. Let us 
suppose that the gap size a can fluctuate within the limits 
a=0—0.01. The coefficients a, and a, are, respectively, 
equal to 2 and 3, The remaining coefficients are deter- 
mined by the parameters of the wing and aileron given 


above, 
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A block diagram of the system described by Eqs.(3) 
is shown in Fig.3. The block with the transfer function 
Y(p) includes the wing and the aileron, that is, it rep- 
resents an oscillatory system with two degrees of free- 
dom; € (t) designates the turbulent input to the system. 


3. Determining the Stable Regions of the 





System in the Absence of Random Disturb- 





ances 

As has already been shown, at first,an electronic 
analog of the system is studied on a type MN-7 analog 
computer, 

In order to check the correctness of the analog, and 
to clarify the effect of the servomotor on system stabil- 
ity, a curve was taken of the dependence of the critical 
speed of the plane, at which flutter arises, on the rela- 
tive aileron frequency n with the aileron control stick 
fixed. In this case, the system is described by Eqs. (2), 
without the terms on the right. 

The stable region is located outside the curve shown 
in Fig. 4 by the dashed line, The results agree with the 
data in [2-4]. 

If the control stick is not fixed, the unstable region 
changes substantially, due to the presence of parasitic 
feedback in the system. We note that the concept of an 
unstable region independent of the oscillation ampli- 
tude has meaning only for a linear system. Therefore, 
let us first examine how the unstable region changes in 
the case when there is no gap in the aileron control sys- 
tem. It appears that the unstable region expands (Fig.4, 
curve 1) if the feedback is positive, which corresponds to 
k = +1, but if the feedback is negative (x = -1), the sys- 
tem becomes more stable, and the unstable region dimin- 
ishes (Fig. 4, curve 2). Insofar as the system is linear in 
the absence of a gap, only soft oscillations can be ex- 
cited in it, 

The stable regions are found, in this case, on an 
analog Constructed according to the block diagram 
shown in Fig. 5, in correspondence with Eqs. (3), but the 
gap is set equal to zero. 

The first input to the summing amplifier 2 is con- 
nected either to the integrator output 15 (x = +1), or 
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Fig. 3. Block diagram of the wing-aileron- 
hydraulic-amplifier system. 











to the inverter output 11 (« = -1), depending on the 
sign of the feedback. 

The presence of a gap makes it possible for hard 
oscillations to be excited in the system being analyzed, 
and for stable limited cycles to appear. 

The performance of the system was studied as a 
function of the aileron relative frequency n for a fixed 
flight velocity, V=0.8. The following phenomena were 
observed: Soft oscillations were excited in the system 
when ny=n<ng, which corresponds to an unstable region 
if the aileron control stick is fixed. If positive feedback 
existed, the oscillations grew toward infinity, but a 
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Fig. 5. Block diagram of an analog of the 
wing-~aileron~hydraulic-amplifier system. 
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stable limited cycle, whose amplitude A,., varied as 
a function of n, appeared if the feedback was negative; 
this is shown in Fig. 6 by the solid line, The system is 
stable for all other values of n, if the feedback is nega- 
tive. If the feedback is positive, then for nj <n <n, and 
Ny =n =n}, where nj and nj correspond to the limits of 
stability in the absence of a gap, hard oscillations are 
excited in the system. The dashed lines in Fig. 6 show 
the amplitude of the minimum initial disturbance, 
which must be imparted to the wing so that it can ac- 
complish unattenuated oscillations, as a function of n. 

It is evident that this quantity characterizes the 
amplitude of the unstable limited cycle. 


4. Measuring the Probability of System 





Excitation in the Presence of Turbulence 








The same method as in [12] was used to measure 
the probability of system excitation as a function of 
time (for given parameters and amplitude of the fluctu- 
ating disturbances), For the same set of parameters, up 
to fifty time intervals were measured during which the 
system remained unexcited, The average excitation 
time tjy and probability p(r) that the system willbe 
excited at time T‘=>T were computed from these meas- 
urements, 

Strictly speaking, the noise of the special generator 
fed to the analog must have a spectrum similar to that 
of the turbulence. However, the spectral components of 
the noise in a narrow band near the free oscillation fre- 
quency of the system play the main role in system ex- 
citation, If the free oscillation frequency did not de- 
pend on the amplitude, then it would be possible for a 
fixed flight velocity ¥V and aileron frequency n to put 
a random input €(r), with practically any spectrum, into 
the system being studied, providing that spectrum had 
the same density as the turbulence at the free oscilla- 
tion frequency of the system, w4, 
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Fig. 6. Limited cycle amplitude as a 


function of the aileron relative fre- 
quency (experimental curves), 
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In reality, a weak dependence of the free oscilla- 
tion frequency on the amplitude is observed in the non- 
linear system being analyzed, This can be seen, for 
example, from Table 1, in which the oscillation fre- 
quencies measured at different wing displacements are 
shown for a number of values of n. 

Here the following designations are introduced; w@ 
is the system oscillation frequency corresponding to very 
small amplitudes; w (A,) is the oscillation frequency 
near an unstable limited cycle having an amplitude Aj; 
w(eo) is the frequency corresponding to very large ampli- 
tudes, 

In view of the weak dependence of the oscillation 
frequency on the amplitude, and, consequently, the 
comparatively small change in the spectral density of 
the noise as the amplitude varies, noise having a uni- 
form spectrum in the region of the system operating fre- 
quency, and whose spectral density agrees with that of 
the actual disturbance at some average frequency Sj, ay, 
can be put into the system. The amplitude of this aver- 
age spectral density as a function of n is also shown in 
Table 1. It is evident that the spectral density of the 
noise fed into the system must be varied in correspond- 
ence with the data of Table 1, as the relative aileron 
frequency n changes, In practice, this spectral density 
change is accomplished in the following manner, The 
signal from an ultralow-frequency noise generator (see 
[12}) is fed into the analog through one of the amplifiers 
whose amplification factor is varied as n changes in cor- 
respondence with the change in spectral density. An 
adequately accurate relation between the amplitudes of 
the noise spectral density is ensured as n changes, in 
this way. The absolute value of the spectral density is 
measured with an accuracy on the order of 25% by the 
method described in [12], The error in this measure- 
ment has no effect, since the absolute value of the tur- 
bulence, as has already been shown, is given very ap- 
proximately, and may change substantially. depending 
on the flight conditions, 

The results of measurements of the average system 
excitation time are presented in Table 1 for a number 
of values of n; the feedback in this case is positive, It 
was necessary to increase the amplitude of the turbu- 
lence spectral density by a factor of four in comparison 
with the data of [8] so that the measurement time would 
not be too long, It has been shown that as n increases 
within the limits nz=n=<n}, the average excitation time 
grows sharply,both as a result of the increase in the am- 
plitude A, of the unstable limited cycle, and because of 
the decrease in the equivalent intensity of the turbu- 
lence, Therefore, the average excitation time in this 
region could be measured at only one value, n= 1,1, 

It is easy to compute the probability of system ex- 
citation as a function of time from the data obtained in 
measuring the average time. The results of the calcula- 
tions for the case when n= 0,55 are shown as a graph in 
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TABLE 1 








n 10° A, @ (0) w (A;) @ (20) 10°S'y, av Tay 
0.5 2.10 | 0.888 0.871 0.873 9.72 180 
0.55 | 1.44] 0.892 0.882 0.894 5.56 {32 
1.4 4.74} 1.125 1.100 1.105 3.82 33 
1.2 2.20 1.4 1.123 1.155 3.53 —- 
1.5 3.60 | 1.295 1.162 1.221 3.06 — 





























Fig. 7. The exponential function exp(-t/Tgy) is shown 
by dashed lines on the same figure. It can be seen from 
the graph that the distribution observed does not have 
the exponential form predicted by theory for the excita- 
tion probability p(r) in the case of sufficiently low noise. 
This difference can be explained by the fact that the 
noise acting on the system is high, so that the time that 
the mapping point spends in the phase space near the 
boundaries of excitation is comparable with its average 
dwell time in the region of the unexcited states, 

One can conclude from this that the methods for 
theoretically computing the excitation probability, 
described in [9, 10] for the case of low noise, do not 
give accurate quantitative results, However, with the 
aid of these methods, one can obtain a theoretical esti- 
mate of the order of the average system excitation time, 
and this has practical value. 


5. Theoretical Calculation of Airplane 





Wing Stability in the Presence of Turbu- 





lence 


It has already been indicated that an analytic solu- 
tion of the stability problem of a complex system 
affected by fluctuations is possible if certain restrictions 
are imposed on the system parameters, namely, oscil- 
lation in the system must be nearly harmonic. The re- 
sults of studies on analogs have shown that this require- 
ment is met sufficiently well over a wide range of 
change in the parameters. 
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Fig. 7. Probability of system excitation for 
n=0,55. 
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The idea of the method presented in this paper is as 
follows, The restriction mentioned above means that 
the solution of Eqs, (3) must be put in the form of almost 
harmonic oscillations with slowly changing amplitudes; 


x= Acosat, y = Bcos(at + q,), 


z = C cos (at + @z). 


Then Eqs. (3) can be linearized in the ordinary man- 
ner by using the harmonic linearization method (the 
Krylov-Bogolyubov method). We obtain an equivalent 
linear set of equations as a result, and their coefficients 
depend on the amplitude of the wing oscillations A. In 
operational form, this linearized set of equations is 


(p? +- ay vp + 1)a+ 
+ (ay4 p® + ay3 Up + a2") y =f, (7), 
(@oq p® + Gog vp) x + (4) 
+ (p? — ay. vp + ag, v2 + n*) y — ay nz = /2 (1), 
—agxC (u)x + (p+ as)z2=0, 


where C(u) is determined from the function ¢¥ (Kk x) = 
= ¥(K AcoswT) expressed as a complex Fourier series; 


@ (xx) = x (u)r+.... 
u=i1 snl oe , Civ) =C,(u) + jC (uv) signa, 


a pase 


! ; eh 
are sinu -- —u VY 1—u- 
1 bs 


1 ° 
Cy (u) = se 00" 1); 
the term a is a gap size, 

Let us now write the characteristic equation for the 
set obtained, for a given A; it takes the form 


L (p) = p> + by p* + dep’ + bsp? + bp + b,=0. 

(5) 

Equation(5) is an algebraic equation with complex 
coefficients, Since the coefficients of the equation 








— —_— aa -— 








specifically depend on the sign of w, the roots of this 
equation must be complex conjugates of the form p = 
=-6(u) + jw(u). Therefore, one need only find the solu- 
tion of Eq. (5) for positive values of w, for example. 

Let us investigate the behavior of the characteristic 
equation roots as a function of the relative aileron fre- 
quency n at a fixed flight velocity, ¥= 0.8. 

Let us assume that 


p = —4(u) + jo + jm, (u) (6) 


for finding an approximate solution of Eq. (5); we will 
define w» so that when p= jw», the modulus of L(p) will 
be maximum. Since the system has been assumed to be 
narrow-band, the following inequalities must be satis- 
fied; 5(u) << wp,and uw, (u) « up, 

Computation and experimental data show that up 
differs little from unity at all values of relative aileron 
frequency n, and at any amplitudes which are of interest 
to us, Therefore, we can say that w»=1, 

By substituting expression (6) into Eq. (5), and re- 
stricting ourselves to terms of the first order of smallness 
relative to 5 and u4, we get the linear set of equations 
which determine 6(u) and ,(u). 

An analysis of the behavior of the attenuation co- 
efficient 4(u) as a function of n gives results which 
qualitatively agree with experimental data, 

For ny=N<Ng, the system is unstable,no matter how 
small the amplitudes, and the oscillations jump toward 
infinity if there is positive feedback. However, a stable 
limited cycle exists if there is negative feedback. The 
boundaries of the unstable region, n, and ng, are deter- 
mined from the equation 


6(u)|,_._, =9. 
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Fig. 8. Amplitude of the limited cy- 
cles as a function of relative aileron 


frequency (theoretical curves). 








The solution of this equation gives n,= 0,612, ng= 
= 0.930, which differs numerically somewhat from the 
data obtained from the analog (n, = 0.565, ng= 1.06). 
This difference can be explained by the inaccurate cor- 
respondence between the analog system parameters and 
the equation being investigated, and also by the approxi- 
mateness of the theoretical calculations. 

For the case of negative feedback, the amplitude of 
the stable limited cycle can be computed as a function 
of n over the range of values obtained for n, by setting 
6(u) equal to zero for a fixed n. 

The results of such a computation of the stable 
limited cycle amplitude are shown in the graph (Fig. 8, 
solid line). 

The range of n over which the system will have 
hard excitation is determined from the equation 


6(u)|,_, =0. (1 


By solving this equation, one can see that hard ex- 
citation of the system is possible only if the feedback is 
positive, for the values of system parameters which we 
have chosen, It should be mentioned here that the sys- 
tem may have hard oscillations for other parametric 
values even if the feedback is negative, 

The solution of Eq. (7) is nj= 0,514, n= 1,08, for 
k =1, Consequently, the hard oscillation mode occurs 
when n has values which satisfy the inequalities; 
nj=n=<n, and ng=n<nj, 

Setting 5(u) equal to zero for a fixed n, we can find 
the dependence of the unstable limited cycle amplitude 
A; onn, graphically. The results of these calculations 
are shown on the graph (Fig. 8, dashed lines) to a scale 
appropriate for the experiment, and are also given in 
Table 2, A knowledge of the unstable limited cycle 
amplitude is necessary for calculating the average ex- 
citation time of the system. 

The oscillations in a complex system (3) are the 
sum of two “natural” oscillations whose frequency and 
attenuation decrement are a function of amplitude, 
Bending oscillations of the wing predominate in one of 
these oscillations, and rotary oscillations of the aileron 
in the other, It appears that the first “natural” oscilla- 
tion attenuates very slowly, whereas the second attenu- 
ates much more rapidly. Therefore, the approximately 
linearized oscillatory system with two degrees of free- 
dom (4) can be replaced by a certain equivalent oscil- 
latory system with one degree of freedom, whose oscil- 
lations are nearly the same as the bending oscillations 
of the wing, that is, they have the same frequency and 
attenuation decrement, 

We get an approximate equation describing the wing 
oscillations as a result of this replacement: 


(p® +> 26p + @*) xz = w° C(t), (8) 
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TABLE 2 





m |10*A,| (0) | @(A;) |10* 8 (0) 


10°Suay| 10°A-"| Yay | fav 





0.55)3.23/0.95) 0.93 | 1.24 
0.6 |1.3910. 96; 0.94 | 0.29 
0.95)1.10)1.13) 1.12 | 1,22 
1.0 |1.41/1.15} 1.13 | 4.83 
1.05/1.96)1.15) 1.16 | 8.76 
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where €(r) is the equivalent noise, whose intensity de- 
pends on the amplitude of the oscillations A, in the gen- 
eral case, 

Now the Fokker- Planck equation can be used if the 
noise correlation time T,,, is short in comparison to the 
system time constant, which is about 1/5, and the noise 
intensity is small enough so that the following conditions 
are met; 


e. U 9- tT 
(Co) "<A, for 1 > = 
{oS (o) | °*< A, for Toor < * , (9) 


where S (w) is half the spectral density of €(r) at the 
frequency w, 

Setting x= AeWT +9), x= AjwelT + %) and substi- 
tuting them in Eq. (8), we get a shorter equation for Ain 
the usual manner (see [10]): 


Fasiys, . 
A= ee j(wt+¢9) (p* =! w*) x =_ 


= — 4(4) + 5 + 2'(0), (10) 


where ¢°(r) is the equivalent white noise having a corre- 


lation function K(t) = 2A7! 8(t), and A“t= “ S (w). 
Since the oscillation frequency depends on the ampli- 
tude, the intensity of the equivalent noise will also de- 
pend on amplitude. 

The Fokker~ Planck equation takes the following 
form in the given case: 


> 


Ow 0 {| a OF in. 

— te a — Ab (A) + — jw ot aaa (A li). 

at dA} 4A \"S OA (11) 
The stationary solution of this equation correspond- 

ing to null probability flow is expressed in the following 

way, as can be seen by direct substitution: 


A 
w(A) = Anexp {— \ 4.43 (A) d Al; 


the normalizing constant Cg is approximately equal to 
5|4=0 = 5(0) if the following conditions are met: 


1/V 260) <A, 1/VRi<A, (12) 
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47 | 1.05 | 10% 15-104 
06 | 0.97 | 1500 | 470 
76 | 0.63 |6-10° 400 
67 | 0.714 |2-1021 |9-105 
55 | 0.84 |5-104/4-104 
where 


. d 
A, = — gy lase (A)} vars 


It is now easy to obtain an expression for the average 
frequency of system transition through the excitation 
houndary, which is the inverse of the average excitation 
time, by a method analogous to the one given in [10]: 


(19 


= 6 (0) AY a exp {— \n.0(-1) dah. 


Hence, a graph of \ Ad(A) as a function of A for dif- 
ferent values of n must be plotted in order to compute 
the average excitation time of the system. Then the in- 
tegral in the brackets of the exponential in (13) can be 
evaluated from this graph, and the value of K, found. 

Calculations show that the equivalent noise inten- 
sity A”! depends weakly on A, the amplitude of the 
oscillations, Therefore, the average value of \, which 
does not depend on amplitude, but is a function of n, 
must be used when the average excitation time of the 
system is computed from (13), The average noise inten- 
sities for different n are given in Table 2, 

Knowing \, the conditions of applicability of the 
theory [the conditions of (12)] can be checked. It ap- 
pears that these conditions are met very poorly. For 
example, when n= 0.55, we have: 1/.6(0)#1/¥K,~3 -107* 
A,*3-10°", Conditions (12) are met even more poorly 
when the noise intensity is increased fourfold, as was 
done in the experiment, This circumstance has already 
been noted in considering the results of system studies on 
an analog. 

The average system excitation timesT,,y, calculated 
from (13) for various aileron relative frequencies n, are 
given in Table 2, The average excitation times T,,, 
for a fourfold increase in noise intensity are given in the 
same table, The data obtained agree qualitatively with 
the experimental results (see Table 1). 

We should remember that the numerical values of 
the average time obtained from theoretical calculations 





tHere ¢° designates the complex conjugate of ¢, and 
the sign < > designates a statistical average. 
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have meaning only as estimates of the quantities men- 
tioned, 
CONCLUSIONS 

1, The presence of a parasitic positive feedback 
path in the control system causes the unstable region to 
expand, In this case, if there is a gap in the aileroncon- 
trol drive, the system can be excited only under definite 
initial conditions (in the absence of fluctuation disturb- 
ance), for certain parametric values. 

2, The danger of growing oscillations of an airplane 
wing arises as a result of the excitation of hard oscilla- 
tions in flight through turbulent air, 

3. The probability of growing oscillations arising 
in a definite time interval depends on the intensity of 
the fluctuations caused by turbulent disturbance of the 
vertical flight velocity component, and this probability 
decreases sharply with a decrease in this intensity. 

4, The basic role in system excitation is played by 
the spectral components of the turbulence near the natu- 
ral bending oscillation frequency of the wing, wp, in 
other words, the amplitude of the spatial spectral density 
of the turbulence, Sj,(Q), at Q=wp/v. 

5. The theoretical methods presented in this paper 
in combination with a study of the system on an analog 
computer permit one to get statistical estimates for the 
excitation time of dangerous wing oscillations as a func- 
tion of the flight velocity and natural frequency of the 
aileron. 

The authors express gratitude to R. L. Stratonovich 
for much valuable advice and suggestions. 
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A method for studying one class of control systems is examined. The method is based on 
constructing vector diagrams and finding the reflection of the parameter space on the 
amplitude -phase plane. 





| 1, Statement of the Problem. Let the transfer function of an open-loop control system have the 
; form 


M 
NW (p) = — Fess (P). (1) | 


where f(p) is the transfer function of a group of links containing no unknown parameters, 


8 r ( 
v= A=1 v=1 A=1 


: 2 
D(p) = S mDAp). ” 


Here m, and ky are the system parameters to be determined, and the coefficients of the polynomials 
Ay(p), By(p), Cyy (p) and Dy(p) depend on the parameters 7y(u=1, 2,... , t) whose precise value is unknown, 
| Let us call 1 the “unreliable" parameters, 


M (p) = 3} m,A,(p) + >) knBa(p) + SY mera (p), 


The transfer functions of a wide class of control systems can be reduced to the form of (1) and (2). In par- 
ticular, the transfer functions of systems which are controlled by two or more physical quantities have such a 
form—for example, the transfer functions of airplane [1] and ship [2] control systems, In this case, the “reliable” 
parameters are those of the regulator, and the “unreliable” parameters are the aerodynamic and hydrodynamic 
coefficients, 


The problem consists of determining the parameters, m, and ky, which ensure stability and the required 
system quality under the conditions that the “unreliable” parameters 1, may vary between definite limits, 


2. Initial Constructions, Let us first assume that the requirements on system quality are merely that 
the system must have a certain degree of stability [3], and let us consider the function Wé€ + iw), representing 
each of the terms in (2) in vector form. 








Having been given an initial system of parametric values 2,(u=1, 2, ... , t), let us construct the follow- 
ing vectors for each frequency w: 


A, (—e + iw, ly) = Sat, 


v==1 


(3) 


where | Az] =a‘ (i,) R"; arg AS = Y (9 A. a Ne a) are the coefficients of the polynomial A,(p); /?. = 


SS —eeeeeEeE——e————eeee—eeEeeeeeeee 











=@ V 1+ ( - ) 7p = arc lg - in, is the highest degree of the polynomial being examined, 


When the parameters v change, the vector Ay receives the increment 


’ 
ly 


TA, = 38 dat Ay, ” 


where 
ce aX {(ly + Als). (lz + Ah), --- (I, + Al,)] — ev (Inv tas + + + 1) 


ee eee 


v 





aY (Iq. Tay. 1) , 


and A ly is the increment of Ly. 


Considering only the linear members of 6 a? expressed as a series in Al yu We represent Sa), in the form 


t ¥ 
Tv . a \' da, 


The sign "0" indicates that the corresponding quantities are determined from the initial values of the para- 
meters, 


The vectors B(p), Cc rv). Dy (p) can be constructed et and their increments found. . 


Next, having been given the initial parameters Ke , we construct geometric series for each frequency 
w, which define the vectors Df-€+iw) and —M(-€+ oh ctie (Cte in correspondence with formulas (2). The latter 
vector is derived from the vector -M(-€+iw) by extending it (—€ +iw)|times and rotating it by an angle arg f 
(— € + iw), 


We will start these geometric series from point 0, and we will denote the ends by the letters D and M, 
respectively; we will indicate the parameters k, and m,,, on which they depend, at each of the component 
vectors, 

k, — k? 
A change in the parameter k, by an amount 6k) = ss a causes the vector 0M to change by an amount 
a 
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Analogously, a change in m,, causes the vectors 0D and OM to change by the amounts 


Am, D = dm,[m° Dy (p)), 


Ba, M = om, | my A, (p) -| y k5, my Cy (p) Ci (P)|, m 


respectively. 
When my and ky change simultaneously, the increment of 0M is found from the formula 





Ay m,M = Ay, M 4-Bm,M + Bajm, i, (8) 


where 
Ay,m, M = Gh, dm, (km! Cy, (p))- 


The vectors plotted for initial vaiues of the parameters are enclosed in square brackets in expressions (6), 
(7), and (8). 


Hence, a change in a parameter corresponds to a proportional change in the lengths of the vectors represent- 
ing it. If the system stability does not meet requirements on system quality, the construction must be done over 
for ==0, 


It follows from (1) and the constructions described above that the ratio of the segments 0M to 0D equals the 
modulus of the amplitude-phase characteristic of an open-loop system (normal for €=0 or expanded €>0), and 
its argument is the angle x between OD and 0M. The ratio of the lengths of MO and MD and the angle between 
vectors MD andMO0 are equal to the modulus and argument, respectively, of the amplitude-phase characteristic 
of a closed-loop system with respect to the controlling input. The real part of the amplitude-phase characteris- 
tic is equal to the projection of M0 along MD, divided by the length of the segment MD. 

















Fig. 2. 








Example, Let us examine an automatic control system; the structure and transfer functions of the individual 
parts are shown in Fig. 1. The transfer function of this system has the form 


(kip -+ ko) (bsp + bo) + ke (map = Meo) (CaP --.Co) Pt, 


W (p) = — p (dap* + dop + ds) (up + Mo) 





A diagram of the vectors 





By (p) = bp + bo, Bi (p) = pBo (p). Cao (p) = cap + cor Car (p) = PCo0 (P)» 
De (p) = dsp? + dap*+dyp. Dy(p) = pDo(p)’ 





where p=iw, is shown in Fig. 2, The diagram corresponds to a negative value of d, and to positive values of the 
remaining coefficients of the polynomials B(p), C(p), and D{(p). A diagram of vectors Da(p) and D,(p), where 
P=-€ +iw, is shown on the same figure for the sake of illustration. The diagram of vector OD = D{iw) and OM= 
-M(iw)e" “3 is shown in Fig. 3 for the frequencies Wo, 44, W,, and W3, Vectors By. By, Sus and rm are rotated 

n radians in constructing the vectors OM, and vectors Dp and D, are rotated wt, radians in constructing OD. The 


amplitude-phase characteristic corresponding to this system of vectors is shown in the same figure (curve 1). 
3. Analysis and Replotting of the Vector Diagrams with Fixed “Unreliable” 





Parameters. We will assume that the “unreliable” parameters are fixed, Let us examine the basic methods 
for analyzing automatic control systems by means of vector diagrams. 





1, The vector diagrams which have been plotted permit one to determine whether or not the system is 
stable with the given initial parameters, or whether it has a given degree of stability. This may be done on the 
basis of the Nyquist— Mikhailov criterion. The number of roots of the denominator with a positive real part (or 
with a real part greater than~€ ) may be found, for example, on the basis of the Mikhailov criterion from the be- 
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Fig. 3. Vector diagram and amplitude-phase characteristic of the open-loop system. 
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havior of vector 5D. Thus,the system of vectors plotted in Fig. 3 corresponds to an unstable system (the denomi- 
nator of the transfer function has one root with a positive real part, since d;<0), 


2, The vector diagrams make it possible to evaluate the effect of parameters m,, and k, on the shape of 
the amplitude~phase characteristic. Thus, for the system considered in the example, increasing k, at all frequen- 
cies causes angle x to increase, that is, the argument of the amplitude-phase characteristic increases, A de- 
crease in m, leads to the same result. Decreasing ky causes the modulus of the amplitude-phase characteristic to 
decrease, and a certain increase in its argument, etc. 


3. The vector diagrams permit one to see the desirability of changing one or another of the parameters 
being examined as a group, and to determine which parameters should be changed and how, in order to shape the 
amplitude-phase characteristic into the desired form. 


The amount of parametric change is limited, naturally, by the physical properties of each parameter. For 
example, the amplitude—phase characteristic of the system must be shaped so that it passes over the point (1, 
Oi), in order to ensure that it will be stable. The evaluations done inSection 2 permit one to determine what 
can be achieved, in particular, by increasing ky and decreasing Ko. 


4, Having been given 6m,, and 5k, on the basis of the evaluations done in Sections 2 and 3, the vectors 
0p and OM, corresponding to a change in the parameters m,, and ky, can be plotted. To do this, one need only 
find the increments of the respective vectors from formulas (6), (7), and (8), and plot them geometrically from 
points M and D, A number of such constructions must be done to find the optimum combination of parameters. 
Constructions corresponding to 5k g=-0.5 and 5k,=2 are shown in Fig. 3 by dashed lines. Curve 2 in this figure 
shows the amplitude~phase characteristic of the system altered in this manner. 


5. The vector diagrams plotted make it possible to decide whether it is desirable to introduce corrective 
links into the control system (including differentiation and integration). The introduction of such links leads to 
an extension and rotation of the vectors proportional to the parameters which are coefficients for the quantities 
being corrected, Thus, for example, the derivative of Y, could be introduced into Fig. 3 to increase the angle 
X. This would be a vector added to point M, perpendicular to the sum of the vectors proportional to kp». 


6. Stable regions (or regions with a given degree of stability) can be plotted for any two parameters, my 
and ky, with the aid of vector diagrams. The stable region boundary points are found on the basis of a coinci- 
dence of points M and D, 


Figure 4 shows the loci of the stable region boundary points tn the plane of parameters ky and m, for the sys- 
tem corresponding to curve fin Fig. 3. Point N of the stable region boundary is located at the intersection of a 
line parallel to the vector ky and passing through point M*, and a line parallel to the difference m, — m, and 
passing through point D. Here, the upper indices (M or D) indicate that the corresponding vector enters the nu- 
merator or denominator of the transfer function. This construction is based on the fact that if the position of point 
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Fig. 4. Plot of the stable region. 











M remains unchanged, the addition of vector 4-6 to point D leads to the same change in the mutual position of 
points D and M as the addition of & to point D, and 6 to point M, The stable region in Fig. 4 is plotted on the 
coordinates 5k9= Ak/A} and 6m,= 4m,/m{, which are defined in terms of a ratio of the corresponding segments, 


If both variants of the parameter enter expression (2) both linearly and in a product, then determining the 
stable region boundary points consists of solving a certain quadratic equation, whose coefficients are easily deter- 
mined graphically (see Appendix), 


7. By using vector diagrams, it is also possible to analyze control systems qualitatively by plotting the real 
part of the amplitude~phase characteristic of the closed-loop system with to the controlling input [4]. 
The construction consists of finding the ratio of the projection MP of vector on MD to the length of segment 
MD at each frequency w. The effect of the system parameters on the shape of the amplitude-phase character- 
istic of the closed-loop system can be traced directly from the vector diagrams, A plot of the real part of the 
amplitude-phase characteristic of the system considered in the example is shown in Fig. 5 (curve 1). It follows 
from an examination of the vector diagram shown in Fig. 3 that it is sufficient to decrease mg in order to in- 
crease the positive region, for example. Corresponding constructions are shown in Fig. 5; the change in the posi- 
tion of points M, D, and P and the corresponding curve are designated by the index 2. 


4. Construction of the Loci of the Ends of Vectors 0M and Op as the "Unreli- 
able"* Parameters Change within Given Limits. Let the parameter ly assume an increment 
Al, . Then the increment of Ao, in correspondence with (4) and (5), will equal 











Ye +. v 3Y 8 at 1 9a, (9) 
y 0 


vy=1 





The increments 4,,B,, 4 and A, D_,, caused by the increment of 1,,, are found analogously. The 
increment of vectors and OD. F M and a can be found from expressions (2), in which A), By, Cyy, 
and Dy are replaced by their es 





It is evident that a change in J, by an amount nAly corresponds to increments of 0M and 0D equal to 
nd, M and nd, D. When several parameters change simultaneously, the increments of OM and OD are found as 
the geometric sums of the increments caused by the change in each parameter individually. 


We will show how to find the correspondence between the subregion of the “unreliable” parameter space 
characterized by the inequalities 1 - 4h, sl u =< an + Ay (w= 1,2,...,t) and the plane region within which 
the amplitude-phase characteristic can be located. 











Fig. 5. Construction of the real part of the amplitude-phase characteristic of the 
closed-loop system. 


967 





ee 





™ 
% 











Fig. 6. Construction of the region { M} . 
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Fig. 7. Construction of the region within which the amplitude-phase 
characteristic can be located. 


es When 1, changes in the range ty- Aly =l, = ty +Al,, the geometric locus of the ends of the vector 
OM(OD) will be the segment formed by laying out vectors 4),M and -A,, MA, D and -A,,D) from point M(D). 
Here points M* and D* correspond to the value ty =f,,. When the two “unreliable” parameters (ff =s;, 82) change, 
the geometric locus of point M(D) will be a parallelogram with its center at point MD"), and with its sides par- 
allel to A MAD). When several “unreliable” parameters change, the geometric locus of points M and D will 
be certain rectilineal polygons with their centers at M' and D‘ respectively. The rule for constructing these poly- 
gons is as follows. 





We lay out the vectors ApM and ~4, M(H =1, 2,..., t) from some point, Let us call the ends of one of 
the segments obtained a, and by. The ends of the segments encountered in going counterclockwise from a, to by 
will be called ag, as, ... , a, and the ends of segments encountered in going from b, to a; will be called bg, bs, 
+ + « » Dy (see Fig. 6), The region being sought ,{ M} ,will be the interior of a 2t-gon obtained by consecutively 
adding the vectors 4b, + a;b2+. » «+ a,b, + byay+ baag+. . . + Dear. 


The region { D} is constructed analogously. 


Having constructed { M} and { D} at all frequencies, we shift them parallelly, so that their centers coincide 
with points M‘ and D’ on the vector diagrams. We call the angle between 0d and 0m a, and the ratio 0m/0d 
will be called r, The point d belongs to {D}, and point m belongs to {M}. The region within which the am- 











plitude ~ phase characteristic can be located on plane W, where w=wj, does not extend beyond the limits of the 
region characterized by the inequalities (see Fig. 7) 


min r (@;) < | W (i@;) | < max r (@;), 
min @ (w;) < arg W (iw) < max a (@;). 


Similar constructions can be done initially for all the given part of the “unreliable” parameter space. If 
the corresponding amplitude-phase characteristic location region does not provide sufficient bases for judging 
one or more properties of a system, the limits of this region can be made more precise by breaking up the “un- 
reliable” parameter space into subspaces, 


CONCLUSION 
The method presented, based on the construction of vector diagrams, permits one to study a wide class of 
automatic control systems, and to solve stability and control quality problems. 


The method permits one to evaluate the joint effect of both “reliable” and "unreliable" parameters on the 
shape of the amplitude-phase characteristic of an open-loop or closed-loop system, 


The construction of vector diagrams can be used for determining the amplitude-phase characteristic of the 
linear part of nonlinear systems being studied, for example, by the methods given in [5, 6]. 


APPENDIX 


Determining the Boundaries of the Stable Region in the Plane of the Parameters Entering the Numerator and 
Denominator Both Linearly and in the Form of a Product, Let us find a boundary point of the stable region at the 
frequency w in the plane of parameters m, and ky, assuming that A,(iw), By(iw), Dgiw) and Cy iw) are not equal 
to zero [see formula (2)]. 








The increments of 0D and 6M caused by a change in these parameters can be written in the form 


_—— —_— - 


AD =hy, AM = mz + hey + rory, 
according to (7) and (8). ; 


Here x and y are the relative increments of ky and m,, respectively, and nj and hj are vectors corresponding 
to the initial values of these parameters. The problem consists of finding values of x and y which satisfy the 
equality 


OD —OM = nyx + nay + rsry, 
where n, = h, — h,. 


Let us draw a line through point M parallel to i, and lines parallel to iz and ng through point D, The inter- 
section of the latter two lines with a line passing through M will be designated by the letters A and B, respectively. 
We will initially assume that A and B lie on one side of M (see Fig. 8). Moreover, let the directions of vectors 
Dy, Dey and 1 ig be such that it is possible to go from M to D along MEFD for positive values of x and y (ME = Tix, 

EF = TigXy, FD= Tigy). Since triangles ABD and AEF are similar, it follows that 


"/ 
onthe 
~~ 








Fig. 8. Figure used in plotting the stable region. 
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pb=eg ™ @ aio! (10) 


where a, b, c, and d are equal to the lengths of AD, AM, BD, and BM, respectively. 
We get Equation (11) for finding x from proportion (10): 





cng d beng 
a+(ft-Z)e— anyng pone J (11) 


and after solving it, y is determined from the equality 


nx —d (12) 
¥=4"G—ad)ns puss Ts ° 


We will assume that the directions of any two of the vectors Ty, Mz, and Ng are opposite to those of ME, EF, 
and FD, Changing the sign of the relative increments (x and y» if opposite to fy and ME, and iz and FD, or only 
X, if opposite to fi, and ME, and jis and EF, or only y if opposite to fi, and FD, and Tis and EF), we arrive at the 
same expression for (11) and (12). 


If only fis is directed opposite to EF, we get expressions differing from (11) and (12) by the sign of the quan- 
tity c, in a manner analogous to that already examined, from the similarity of triangles ABD and AB'F', By ap- 
propriately choosing the signs of the relative increments, we can arrive at the same situation as when the direc- 
tions of any one or all three of the vectors are opposite to those shown in Fig. 8. 


When A and B lie on different sides of M, x and y are found from proportions differing from (10) only by the 
sign of d, as is clear from Fig. 8, In this case, the relative increments of m, and k, defining the system at the 
boundary of the stable region are found in all cases from relations (11) and (12). The quantity c> 0 if the direc- 
tions of Tig, Tig,and fig agree with the directions of MA, AD, and BD; there are two different directions,and c<0 in 
the remaining cases. The quantity d>0 if points A and B are on one side of point M. The increment of ky is 
positive if the direction of fy agrees with that of MA; the increment of m, is positive if the direction of fi, agrees 
with that of AD. 
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The paper studies certain problems associated with the design of machines that can classify 
objects according to an indentifying feature which is not specified a priori. A number of 
general conditions are indicated for which this problem can be solved by machines with 
"stimulation" (machines which receive certain additional information from the operator). It 
is demonstrated that in a number of cases machines without stimulation can be used for the 
same purpose; a brief analysis is made of the capabilities of one machine of this type which 





has already been constructed (the “Perceptron”). 


1. Introduction 

In designing various types of self- adaptive machines, 
the first problem is to design a machine that classifies 
input situations according to an identifying feature which 
is not known a priori. An example of this kind may be 
found in a high-speed computer which changes or formu- 
lates its program during the solution of the problem. 
After this machine has solved a number of problems ,it 
will have accumulated experience which is manifested 
in the set of programs it has formulated. In solving a 
new problem the machine must first determine whether 
or not this problem is “similar to some other problem 
it has solved previously; if this is so, then it must at- 
tempt to use a program which is already available. 
However, it is not known a priori which problems the 
machine must assume to be “similar.” Thus, the expe- 
rience of the machine must include not only a set of 
programs for solving specific problems, but also a prog- 
ram which is developed during the process of operation 
and makes it possible to separate the entire set of prob- 
lems which are to be solved into classes of similar prob- 
lems. 

Another example is a machine designed for control- 
ling a technological object according to a program which 
is not specified a priori but is developed on the basis of 
the actions of an operator under various conditions, But, 
as a rule, there is no possibility of recording the actions 
of the operator in the memory of the machine under all 
conditions that might conceivably arise—these conditions 
are too numerous. In that case the machine must sepa- 
rate the entire set of different combinations of instrument 
readings into classes corresponding to identical required 
control actions. 





The basic features of the classification problem are 
often illustrated by the means of the following very clear 
example. Assume that photoelements, which act as the 
detectors in the machine, “scan" a certain field. The 
operator (human or machine) “shows” the machine geo- 
metric patterns in arbitrary sequence; these patterns have 
been subdivided by the operator into two classes accord- 
ing to some identifying feature, For example: 1) patterns 
located in the left and right sides of the field; 2) patterns 
located in the upper and lower parts of the field; 3) pat- 
terns with large and small areas; 4) squares and triangies; 
5) patterns that are symmetrical relative to the horizontal 
and vertical axes, etc, Two situations are possible at 
the output of the machine (for example, the appearance 
of "0" or "1"). No separation of the input situations into 
two classes has been fed into the machine a priori. It is 
required that after a certain self- adaptive period the ma- 
chine will be capable of producing a "1" at the output 
each time that the machine is shown a pattern belonging 
to some particular class, and of producing a “0” for the 
appearance of a pattern belonging to a different class, 

Since the machine “does not know” a priori which 
identifying feature has determined the separation of the 
patterns into two classes by the operator, it will err at 
first, no matter what its algorithm for operation may be. 
This means that one and the same situation at the out- 
put of the machine (for example, the appearance of a 
"1") will arise for the appearance of patterns belonging 
to different classes; conversely, different situations at the 
output of the machine may correspond to different pat- 
terns belonging to the same class, The algorithm must 
be such that after a certain number of “viewings” the 
machine begins to operate without error, or at least in a 


971 








— en — 


or er ee 


manner where the probability of error is small. We shall 
assume that when a self-adaptive machine has reached 
this quality of operation, it has developed a “comprehen- 
sion” (in the example given here, the comprehensions 
that are developed are “right” and “left;” “upper” and 
"lower," etc.). 

Showing another pattern to the machine corresponds 
to the appearance of a certain situation at the input of 
the machine (the excitation of certain photoelements 
and the nonexcitation of other photoelements), Each 
specific sitation at the input will sometimes be called 
an'dbject.” 

In the general case, the machine develops a com- 
prehension and then converts the set of input situations into a 
set of output situations. In our example, the set of out- 
put situations consists of two elements (0 and 1). The 
machine must have a series of variable parameters for 
which each value completely determines the transforma- 
tion law carried out by the machine, We shall call each 
set of specific values for these variable parameters a 
“state” of the machine. 

Initially we shall assume that the operator feeds the 
machine a specific signal each time the machine com- 
mits an error (i.e., the operator feeds additional informa- 
tion into the machine). The state of the machine may 
vary according to a definite law (in the general case it 
is not necessarily a deterministic law), depending on the 
input situations and on the additional information re- 
ceived from the operator. Such machines shall be called 
"adaptive machines with stimulation.” 

Below we shall demonstrate that machines with 
stimulation are not the only class of adaptive machine. 
Under certain definite conditions, machines can develop 
comprehensions without using the additional information 
from the operator. 


2. Certain Conditions Which Assure the 





Development of Comprehensions by Ma- 
chines 

In the general case the operator can subdivide the 
set of input situations into arbitrary classes (comprehen- 
sions). In principle, the machine may err until all pos- 
sible objects in the given experiment have appeared, and 
the operator has in some way informed the machine what 
output situation must correspond to each possible object. 
But the method of “total scanning,” which is the only 
one that can solve the problem of developing compre- 
hensions in this general case, is inexpedient if the num- 
ber of objects is sufficiently great,and if it is necessary 
to solve the problem in a limited time with a limited 
memory. 

Evidently, the only method by which total scanning 
can be avoided is to limit a priori the set of allowed 
classes into which the set of objects can be subdivided 
by the operator. 

We shall examine certain examples of the limita- 
tions which can be imposed on the allowed classes of 
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objects and shall examine certain concepts concerning 
the design of machines which develop comprehensions 
for these limitations. For simplicity we shall assume 
that all possible input situations are always subdivided 
into three classes by the operator: a) objects correspond- 
ing to one comprehension; b) objects corresponding to a 
different comprehension; c) objects which do not appear 
in the given experiment. In other words, the machine 
must develop two mutually complementary comprehen- 
sions for the condition that a portion of the objects does 
not appear. 

Example 1. Assume that the input of the machine 
is formed by r detectors whose state is described by the 
vector (Uy, Ug, . . . , Up); assume further that n arbitrary 
functions of the state of the detectors are known: 


% = f, (Uy, Ue,...,Up) (6=1,2,...,n). (1) 


Each object corresponds to a definite point with the 
coordinates X4, Xz, . . . » X, in the space {x}. We shall 
impose the following limitation on the allowed classes 
of objects. In the space {x} it is possible to select two 
regions E and E* such that: a) points corresponding to 
objects associated with one comprehension are confined 
to one,and only one, region, and points corresponding to 
objects associated with different comprehensions are 
confined to different regions; b) there exists a set of 
hyperplanes such that the regions E and E* lie on differ- 
ent sides of any hyperplane in this set. 

We shall turn our attention to the fact that the re- 
gions E and E* are not fixed a priori and may be differ- 
ent for different classes of objects. 

Assume the equation for a hyperplane in {x} is 
written as 


6 = Dy Mati + Ang = 0. (2) 
i=g 


If the hyperplane o=0 passes between the regions E 
and E* it follows that the linear function 


6 = >) Mats + Ants 


t=1 


acquires positive values in one of the regions and nega- 
tive values in the other region. 

We shall study the (n+ 1)-dimensional space { } 
with the coordinates A;,A 9,..., net" We note that 
if a certain point (X14, Xe. oe mm ee 1) in this space 
defines the hyperplane 


n 
$= Dd Ria + Tots =0 
i=} 


in the space { x}, it follows that any point in {\ } 











lying on the straight line connecting the origin with the 
point (Xi, Xe, . « +» Mpg) Will define the same 
hyperplane, Therefore, the set of separating hyper- 
planes in the space { x} will correspond to the set of 
straight lines in the space {A } emanating from the 
origin and filling a certain hypercone. The problem of 
developing comprehensions reduces to transferring the 
current value of the vector A = (A,;,Aq,.. -Am,,) into 
the indicated hypercone. The number of errors during 

a series of trials (a series of “scans’ of different objects) 
can serve as the criterion of the distance between 
and the hypercone. For example a machine classifying 
objects at its input for the condition that the correspond- 
ing regions E and E’ can be separated by a hyperplane 
can operate according to the following algorithm. When 
a new object appears the machine determines the values 
of all Xj from (1) (all x; are known functions); then 
the value of the linear function o is determined for the 
current value (Ay, Aq, . . . An+4), and a signal 0 or 1 
appears at the output of the machine, depending of the 
sign of o. From additional signals received from the 
operator the machine counts up the number of errors for 
the series of trials and alters (for example, randomly or 
in some other fashion) the value of (A, A»... . Anes 
then it again determines the number of errors during the 
series of trials, etc., until the number of errors becomes 
minimal. 

An algorithm of this type is the basis for a machine 
designed by O. Selfridge [1] and called “Pandemonium.” 
The machine consists of several computer (or converter) 
stages (Fig. 1)*. The first stage consists of detectors D 
which perceive the objects. The detector outputs supply 
the computer devices C of the second stage that form 
certain functions x;, X_,... , Xp for the state of the 
detectors which are known a priori. The output from the 
ith computer device is applied to the corresponding am- 
plifier A with a variable gain A;, The outputs from all 
the amplifiers are summed, and the sum 


n 
ee > Mery 
i=1 


is applied to the decision block DB which produces 1 or 
0 depending on the sign of 0. If “Pandemonium” errs, 
then the output applies a signal to the control unit CU. 
The number of errors in the operation of Pandemonium 
over a sufficiently long series of trials F(A;, A, ..., Ap) 
is a function of the coefficients A; The control unit 
varies A in such a way as to minimize the number of 
errors during a series of trials; the variation is performed 
according to the method of steepest descent. Thus, the 
control ‘unit is essentially an external controller with 
"regulating organs" Ay, Az,... , Ap, that minimize the 
function F(Ay, Ag, ... , A,). 

In the example involving a machine which must 
classify geometric patterns projected onto the photo- 
element field xj may describe the state of ith photo- 


element. Assume that when the photoelement is excited 
the corresponding coordinate acquires the value 1, and 
that the value 0 is acquired when it is not excited, 
Then the act of “showing” any pattern to the machine 
corresponds to a point in the space {x} lying on one of 
the vertices of a unit n-dimensional hypercube (n is the 
number of photoelements). It is easy to see that when 
the machine develops such comprehensions as "left" and 
“right,” “upper” and “lower,” “large” and “small,” we 
can find corresponding regions E and E', However, in 
this space it is impossible to develop such comprehen- 
sions as a “square” or a “triangle,” “a pattern symme- 
trical with respect to the vertical axis,” and "a pattern 
symmetrical with respect to the horizontal axis.” 
Example 2. It was assumed above that only one 
space for performing the separation was known to us 
a priori. However, it is also possible to encounter 
the following case. Several spaces are specified 
(several sets of functions for the state of detectors) 
along with the following limitation on the allowed 
classes of objects: two regions E ande’ having the 
properties indicated in the first example are found 
in at least one of these spaces. 

The problem of developing comprehensions reduces 
to a search for the required vector in the spaces which 
are successively formed. The following condition may 
be used to govern the transition to trials in the next 
space: the absence of a variation in a certain current 
value of the vector A in the sampled space such that 
the number of errors during the series of trials will be 
reduced, Since the desired variation in the vector A 
must be determined each time on the basis of a finite 
number of steps, it follows that there may exist a cer- 
tain probability that transition to trials in the next 
space will be performed erroneously. Thus, in the case 
under study the machine may develop a comprehension 
only with a certain probability. 

As we have already noted, in our example the photo- 
element field cannot develop the comprehensions: a 
" square” and a “triangle,” a “pattern symmetrical 
with respect to the horizontal axis, “ and a “pattern 
symmetrical with respect to the vertical axis? if the 
search for the separating hyperplane is performed in a 
space where each coordinate corr-’sponds to the state of 
one photoelement. However, it is possible to organize 
another space {x } in the following fashion. Assume 
that the field of the photoelements forms a matrix and 
that the excitation of each photoelement causes the 
appearance of a pair of numbers (j, k) at the input of 
the machine; here j and k respectively denote the num- 
bers of the row and column which locate the particular 


* The Pandemonium can develop several comprehensions 
simultaneously. In this paper we will describe a simpli- 
fied Pandemonium network designed for developing just 

two mutually complementary comprehensions, 
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photoelement. In each row we take the excited photo- 
element with the smallest column number k(j),,;, and 
the excited photoelement with the largest column num- 
ber k(j)max, and then form the following half-sum for 
each j: 


- 1 
ki) = kW may t EW 


Then we find the maximum half-sum with respect to 
k 
i: 7 k(j), and choose fj = mex as function of 


the input situation, If not a single pair of numbers with 
any j appears at the input of the machine when a cer- 
tain object is “shown,” then we shall assume that the 
value of the corresponding function is ff» = 1. 

If the photoelements form a matrix with n rows 
and m columns,we shall have n functions 


min 


k - 
j 
Analogously, we form m functions 
7 (k) 
Ink = ax] (ew i,2,..., m). (3b) 
k 


We shall study the space {x } with the coordinates 
xj = fj (i=1, 2,..., +m). In this space it is already 
possible to isolate patterns which are symmetrical with 
respect to the horizontal axis, patterns which are symme- 
trical with respect to the vertical axis, squaresf , and 
triangles, 

These examples demonstrate that the manner in 
which a machine with stimulation must operate is deter- 
mined by the limitation imposed on the allowed classes 
of objects,t 


3. Machines without Stimulation 





Recently it has been determined that when certain 
additional limitations are imposed on the allowed 
classes of objects (these limitations are based on the 
statistics which govern the appearance of the object), it 
is possible to design machines which separate objects 
into classes without stimulation. 
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Up to now we have studied the following problem: 
the operator subdivides the entire set of objects into two 
classes; the machine, after a certain adaptive process, 
is required to subdivide these very objects into classes 
which coincide more or less exactly with the classes 
specified by the operator. Here we noted that the 
classes into which the operator subdivides the set of ob- 
jects must satisfy certain limitations, 

Now we shall study a different problem, A set of 
objects appears at the input of the machine; these ob- 
jects are objectively (independently of the desires of the 
operator) subdivided into classes which satisfy the same 
limitations as those given in the preceding problem. 
The machine must subdivide the objects into these 
classes and provide the operator with information on 
objectively existing classes which were not known to the 
operator. 

Both of these problems can be posed with respect 
to a machine without stimulation. And,although in 
principle the second problem is of greater interest, it is 
true that essentially the same machine can be used for 
solving both the first and second problems. 

In order to clarify how either of these problems can 
be solved,we shall study the following simplest example. 
Assume that each object is made to correspond to a cer- 
tain number x; assume also that on the x axis it is pos- 
sible to find two nonintersecting segments [x;, x,] and 
[xq, x] such that the numbers corresponding to objects 
of the same class belong to only one segment, and the 
numbers corresponding to objects of different classes 
belong to different segments. 


t It is assumed that the axes of symmetry for a square 
coincide with one of the rows and one of the columns of 
the photoelement matrix, 


It should be noted that the general problem of develop- 
ing comprehensions is closely associated with the far 
more difficult problem of formulating the space proper, 
as well as with the problem of separating two sets in a 
specified space,which is treated in this paper, This 
problem has evidently not been posed or solved in any 
acceptable form as yet, 
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Assume further that the probability-density distrib- 
ution f(x) for the appearance of objects with the coordi- 
nate x is specified, Then it is possible to establish rela- 
tionships between the numbers xq, Xg, Xy", X_" (Fig. 2) 
and the properties of the probability-density function 
f(x),which are sufficiently general so that when the ma- 
chine operates by comparing the average value X with 
the current value x (Fig, 3),it is capable of developing 
the specified comprehensions without any additional in- 
formation from the operator.** This principle can be 
generalized to apply to multidimensional spaces, 

The first to propose a self-adaptive machine with- 
out stimulation was F, Rosenblatt [2]. His machine 
(which he called the "Perceptron"t Tf ) operates accord- 
ing to a principle that differs from the one cited above. 
The shortcoming of the principle proposed by F. Rosen- 
blatt is that the limitations imposed on the allowed 
classes of objects for which the Perceptron can operate 
are not natural ones and cannot be interpreted in con- 
venient terms, However, since the Perceptron is the 
first constructed machine without stimulation, the sub- 
sequent sections of this paper will analyze the principle 
of operation of the Perceptron and certain problems in 
the theory of this operation. 


4, The Perceptron Networktt 








A Perceptron (Fig. 4) consists of a set of identical 
elements which are conditionally called neurons. A 
certain function x of the state of the detectors is applied 
to the input of each neuron at discrete instants n=0, 1, 
2,..., Corresponding to the instants at which the ob- 
jects appear, Each neuron is an amplifier with a vari- 
able gain A that changes at these same discrete instants. 
The increment of gain during one cycle is determined 
from the formula*** 


AA” = ANH — 4" = (am — ka") p", 


where p is the signal produced at the output of the next 
Perceptron stage, and k is a certain constant. The out- 
put of the neuron at the nth instant is the quantity rx". 

All neurons are subdivided into two groups. Assume 
that there are / neurons in the first group and m neurons 
in the second group. We shall prime those quantities 
which refer to neurons of the second group. The vari- 
able p is the same for all neurons in one group. Thus, 
we obtain! + m relationships: 


fiz} 


ix a 


Fig. 2, 





AM = (af — kA?) p” (= 4, 2, 2.50 


(4a) 
Adj" = (a;"— kjhj") p™ ,m). 
The outputs of all the neurons are summed; the out- 
puts of the neurons in the first group are taken with a 
plus sign, while the neurons in the second group are 
taken with a minus sign, The resulting sum 


l m 
on = DS} Ate? — Dapey 
i= j=t 


(j=1,2,... 


is applied to the input of the decision element DE at 
each instant; the decision element has two outputs p and 
p’ and operates according to the law S 


2 1 o*>0 
hess a" < 0, 

( 4b) 
p"=1—p*. 


Note that if the sequence of values acquired by 
the vector x"=(x,", xf, ..., xfs x;9, oone xP) and the 
initial value of the vector A°=(A4, A}, see» doe Aas 
coos wm) are known, then Eqs. (4a) and (4b) completely 
define the value of the vector \" and the outputs p” and 
p'" of the decision element at all instants n=0, 1, 2,... 
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Fig. 3. 
**For example, a sufficient condition is the validity of 
the inequalities Xe Xs 
a1\ f(e)dz +2, | s(2)d2> 20 
m ss 
i 


x, 
2,\ f (x) dz +- zs f(z) dr <«,. 


ee 


Xe 
Here \ j(z)dz is the probability that an object be- 
*9 
x 
longing to the first class will appear, and \ { (x) dz 


is the probability that an object belonging td the sec- 
ond class will appear, 

tt From the English word " perception” —perception, 
recognition, 


+t The presentation of the Perceptron properties in this 
and subsequent sections differs from the corresponding 
treatment in [2]. 


©**Here and throughout, the superscripts n, n+1, etc., 
indicate a corresponding instant; they are not exponents, 
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We shall study the ! + m)-dimensional space { x } 
with the coordinates x), x9, ...) Xz» Xj» ese» Xn )s 
Each object appearing at the input corresponds to a point 
in {x}, At each instant the equation 


l m 
a= Dy Ain— >) A;'2; = 0 

t=] j=1 
specifies a hyperplane in {x } whose directional cosines 
are determined by the neuron gains (i.e., by the value of 
the vector A"), Since during the process involved in the 
appearance of the object at the Perceptron input A" 
varies in accordance with (4a) and (4b), it follows that 
the hyperplane o=0 is “rotated” in some fashion. 

At each instant n the appearance of p"=1 at the 
Perceptron output corresponds to the appearance of an 
object whose coordinates lie on one side of the hyper- 
plane; the appearance of p’"=1 at the Perceptron output 
corresponds to the appearance of an object whose coor- 
dinates lie on the other side of the hyperplane, If,after 
a sufficiently large number of objects appear,the posi- 
tion of the hyperplane is stabilized (i.e., it merely 
"rocks" near a certain position), then the appearance of 
p=1 or p=0 at the output of the machine will determine 
the subdivision of the output into two classes in a more 
or less stable fashion; only objects whose coordinates are 
located in the “rocking” region will correspond first to 
p=1, then to p=0. Assume (as we have done throughout 
this paper) that the objective difference between the 
classes of objects resides in the fact that each class cor- 
responds to a region in the space {x }; here these re- 
gions are governed by the limitations in Section 2. If 


= 





4 
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we 
» 


Fig. 5, 


976 


the subdivision of objects into two classes by the Percep- 
tron corresponds to the separation of these regions, then 
this will indicate that the machine has developed a com- 
prehension, 

In the analysis below we shall demonstrate that the 
Perceptron is far from being able to develop a compre- 
hension in this sense for all cases, 


5. Certain Problems in the Theory of a 





Perceptron 





We shall derive certain relationships which deter- 
mine the process involved in the time variation of the 
vector A. For simplicity we shall assume that there are 
two neurons with the gains A and A‘ and the inputs x and 
x", However, it is obvious that if we treat x, x’, A, A’ 
as vectors, then the basic results remain valid. 

Assume that the plane {x } contains two regions 
E and E* which can be separated by the straight line 7, 
passing through the origin; then the set of these separat- 
ing straight lines form the angle a, (Fig. 5). As we 
have already noted in Section 2, each separating straight 
line 2, in the plane {) } corresponds to a straight line 
ly passing through the origin; the set of separating lines 
which fills the angle a, on the plane {x } corresponds 
to the set of straight linesin the plane { } that fills the 
two opposing angles &, and a} (a “two-dimensional 
cone"). 

The problem which the machine must solve reduces 
to separating objects whose coordinates lie in the regions 
E and E*, Since all of the straight lines which separate 
these regions lie within the angle a, on the plane {x}, 
it follows that the machine develops a comprehension 
only if the values of the neuron gains A and A‘ lie either 
solely within the angle c, or solely within the angle a\ 
on the plane {A} for most of the time after the adapt- 
ive process has been completed, Under these conditions 
the point (A, A") must not pass from & y, into a, (or 
vice versa after the adaptive process has been com- 
pleted; this holds because if such a transition were to 
occur, the sign of © would be changed, and on the basis 
of (4b) the alternate values p=1 and p=0 would corre- 
spond to the same region E (or E*); this is inadmissible. 

Before we pass on to further derivations we shall 
study the simplest example which will show that in cer- 
tain cases the Perceptron is actually capable of develop- 
ing a comprehension (i.e., it is capable of performing 


























the necessary separation); it is also possible to encounter 
cases when the Perceptron is not capable of developing 
a comprehension. 

Example, We shall study a Perceptron in which 
k=k’'=1, In accordance with (1), A=x and A*=x' at 
each instant; here x is the x coordinate of the nearest 
preceding object for which o =0, and x’ is the X coor- 
dinate of the nearest preceding object for which 0<0, 


We shall superimpose the planes {x } and {A} and 
shall assume for simplicity that the angle &, is symme- 
trical relative to the bisectrix of the first quadrant, 
Then it is not difficult to see from the equation 

o = Ax—A'x'=0 that the anglea y coincides with the 
angle &, when the planes are superimposed (Fig. 6a). 

If the separation has been performed, then the point 

(A, A") must lie only in one of the shaded rectangles, 

It is obvious that if we shift the regions E and E* in such 
a way that the angles &, arid &) are not altered, it is 
possible to obtain shaded rectangles both inside and out- 
side the angle @, (Fig,6,a and b). However, if no rect- 
angle lies inside the angle a), then it is certain that 
separation is impossible. 

We shall now continue our study of the general case. 
We shall assume that the sequence in which objects 
appear is random and that we know only the probability- 
density function ¢(x,x") for the appearance of an object 
with the coordinates (x, x"); the probability that a cer- 
tain object will appear at a certain specified instant 
is not dependent on what has appeared previously. 

Note that in this case the probabilities that the 
point (A", A*") will be shifted to the point (A, \*n+4, 
depend solely on the current value of (A", A") and are 
determined by the function ¢(x, x"); i.e., the process 
consisting of the variation of (A, A*) is a random Mar- 
kov process, 

Assume that we know how to determine f,(A, °) 
(this function is the probability-density function for 
(A, A‘) after n instants}, 

For each appearance of a new object a variation 
can occur either for A only or for A* only, since either 
p=0 or p=1 is mandatory. Therefore, only those points 
(A, A’) for which either A=A or A*=A" can arrive at 
the point (x, *) during one step. Here the point 
(A, *) is transferred to the point (X, X) if an object ap- 
pears with the coordinates (x, x") satisfying the condi- 
tion 0 = Ax—)'x" =0(p=1) and x= Au —(1—k)A; this 
result is valid on the basis of Eq. (4a). Analogously, the 
point (X, A’) is transferred to the point (X, x) if an ob- 
ject appears with the coordinates (x, x") which satisfy 
the condition 6 = Ax—A'x'<0 (p=0) and x"= 
rN —(1— kh"). The probability density for transfer from 
the point (A, *) to the point (X, ®) is therefore equal to 


i (0, %, %) = lo GR—A—WA, zyar, 
KN <(K—(A— kA, 


where the integral is taken over all x’ satisfying the 
condition 


ah <(kK—(A—k) ga. 


~ rhe probability density for transfer from the point 
(A, A’) to the point (A, A’) is equal to 


wa(d’, 2, 0) =\p(e, ¥ —(1 — ke) ’) de, 
th < fi’ —(1— kr‘, 
where the integral is taken over all x satisfying the 


diti *~ oa 
"EK cl a ae, 


The recursion relationship which determines the 
Markov process in the case undet oA is written as 


fas (h, *) = Vin0 X’) pi (a, Ky Xda + 


fn, var’, &, Pyar’, 


where the integrals are taken over all possible values 
of A and A’, and the functions py and #, are determined 
according to the expressions cited above. 

The probability that the point (A, A’) will be located 
within the angle @) after the nth instant is equal to 


P, (Qa, W)€a) = {I In(d, \’)ddad’, (6a) 


@) 
The corresponding probability that the point will 
be located within the angle a’, is equal to 


Pa ((s, ¥)€a3) = (lta, Wy daa’, 

aa 

Here the integrals are taken over the regions be- 
longing to the angles &) and a’) , respectively. 

It is evident that the Perceptron develops the re- 
quired comprehensions if after a sufficiently large num- 
ber of different objects haveappeared, there is established 
a sequence of distributions f,(A, A’) such that either the 
value of Py or the value of Pj, is close to 1 in the pre- 
dominant number of cases. Quantitatively this condi- 
tion can be expressed as follows: for each initial value 
(A°, A*°) we find a value of N such that for all n>N 
either the inequality 

n 


(6b) 


Pn (7a) 
a—Nei >i —* 
or the inequality 
P,, 
(7b) 


s—Net 
will be satisfied. 
It is possible to require that a more rigorous con- 
dition be satisfied: for each (A°, A'°) we must find an 
N such that for all n>N we have either 


“ P,=1 (8a) 
P= 1. (8b) 
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The validity of one of these conditions means that 
in the final analysis a sequence of probability-density 
functions f,( A, A") is established such that the point 
(A, A’) never emerges beyond the limits of either the 
angle &) or the anglea, therefore the probability of 
error in the Perceptron becomes equal to zero. 

It is considerably simpler to predetermine the 
ability of the Perceptron to satisfy one of the conditions 
(8a) or (8b) rather than one of the conditions (7a) or 
(7b). Assume that the adaptive process has been com- 
pleted and that one of the conditions (8a), (8b) is sat- 
isfied, This means that the straight line 0 =0 passes 
rigorously between the regions E and E’, Therefore, in 
view of (4a), the value of \ can change only at those 
instants when (x, x’) € E, and that of A‘, only at 
those instants when (x, x") € E*; or conversely, the 
change is dependent on whether the point (A, A") is 
within the angle «, or within the angle a‘, [i.e., 
whether condition (8a) or (8b) is satisfied], Assume, for 
' example, that condition (8a) is satisfied, and that 
changes in those cases when (x, x") € E (the remaining 
cases are treated analogously), 

We shall now study only those instants s= 0, 1, 2,..+0 
at which an object with the coordinates (x, x") appears, 
instead of studying all instants; therefore, we assume that 
the value of A changes, We shall choose some instant 
after the adaptive process has been completed as the 
initial instant; under these conditions the point (A, A‘) 
is already within the angle a). 


In accordance with the difference equation (4a), 
s—1 


M= >} (4—k)p ai + (4 — ky not tt (9 
j=0 
The variables x) are independent random quantities 
with identical probability-density functions that are 
equal to the conditional probability-density function for 
the condition (x, x’) €E: 
{9 (=, =) ae" 





{5 p(x,2’)drdz' ’ 


since the instants s= 0, 1, 2, . . are by definition (x, x’) € 
€E at these instants. 

Since A is a linear function of independent random 
variables, it follows that 


s—1 
M = >) (4 — ky + (4 — eye, (10) 
j=0 


where the bar above a quantity denotes mathematical 
expectation, The series (10) converges only for the 
condition 0< k < 2; here 

lim K* =% = 42. (11) 


8-cc 


From (9) it is also easy to derive the expression for 
the dispersion of the limiting probability-density func- 
tion; 


lim D (A*) = D(a) = Ee (2). (12) 





Thus, knowing the arrangement of the regions E and 
E’, it is possible to find the angles ay and a> ; knowing 
Y(x, x"), it is possible to determine 4, D(A) and to 
estimate whether \° is located solely within the angle 
a 
r* 











mate whether AS is located solely within the angle @)- 

We should pay attention to the fact that the disper- 
sion D(x) in (12) is computed only for the region E, 
Equation (12) illustrates the part played by the feedback 
loops for the quantities p and p’ (Fig. 4). If these rela- 
tionships do not exist, then the quantity A will change 
for the appearance of objects from both the region E and 
the region E'; the dispersion D(x) will then be computed 
over both regions, This will cause a pronounced in- 
crease in the dispersion of the limiting distribution and 
will mean that the probability of the emergence of 
the point (A, A‘) beyond the limits of the bounded angle 
a, will have increased, 

We shall now pose the following problem: the oper- 
ator has a priori subdivided the objects into two classes 
which satisfy the limitations cited in Section 2, Is it 
possible always to choose the statistic involved in the 
appearance of the objects in such a way that the Percep- 
tron will perform the separation conceived of by the 
operator ? 

An analysis of the Perceptron operation shows that 
for certain additional conditions the Perceptron is capa- 
ble of performing the separation conceived of by the 
operator (i.e., it is able to develop the required compre- 
hensions), However, these conditions are not obvious, 

The example cited above demonstrates the possibil- 
ity that there may be no statistics of object appearance 
for which the Perceptron can perform the required sepa- 
ration, 

In the general case the capacity of the Perceptron 
for developing the required comprehensions is determined 
by its structure (the value of the quantity k and the spe- 
cific form of the function x), by the nature of the ob- 
jects which are "shown" to the machine (the arrange- 
ment of the regions E and E"), and by the statistics gov- 
ering the “showing” process (the form of the probability- 
density function for the appearance of the objects), In 
order to determine this capacity a priori for a specified 
structure and specified "showing" statistics,it is neces- 
sary to perform a computation using formulas (5), (6), 
(7); this process in itself requires the use of mathemati- 


cal machines, LTERATURE CITED 
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The paper studies a network for summing the fine and coarse control signals; this network 
can be used to control a servodrive while simultaneously monitoring the accuracy of auto- 


matic tracking with only one null indicator, 


In existing servodrives semiautomatic control and 
monitoring of the automatic tracking accuracy are 
achieved with two null indicators that are connected in 
parallel with the selsyn receivers (or other detecting ele- 
ments), 

The null indicators have different sensitivities 
with respect to the mismatch angle, and their pointers 
deviate in direct proportion to the envelope correspond- 
ing to the voltages from the fine and coarse detector 
elements: 


U,,, sin ® U ,, Sin kO 
a9 af = a 





a = S§ 


where Q¢ , G@¢ are the deviation angles for the pointers 
on the null indicators for coarse and fine control, re- 
spectively; Up, is the maximum output voltage of the 
detector element; @ is the mismatch angle; R, and R, 
are the resistances in the signal- winding circuit of the 
null indicators; S is the slope of the characteristic for 
the null indicators; k is the transfer number between the 
detector elements for fine and coarse control. 

The semiautomatic control of a servodrive using 
two null indicators is based on using an appropriate de- 
vice (a semiautomatic control device) to achieve a con- 
stant alignment of the pointers of both null indicators 
with their zero positions, The basic difficulty in control 
of this type and in monitoring the tracking accuracy re- 
sides in the necessity for simultaneous observation of the 
pointers of two meters. In the analysis below we study a 
network for summing the fine and coarse signals; this 





network makes it possible to control the servodrive and 
to monitor the accuracy of the automatic tracking using 
a single null indicator. 

The operation of the network (Fig, 1) is based on a 
nonlinear volt-ampere resistance characteristic; the fine 
control signal and a portion of the coarse control signal 
are limited by the nonlinear resistor NR formed by two 
solid-state rectifiers at the input of the null indicator, 
For a small mismatch angle the pointer of the null indi- 
cator is deflected basically by the fine control signal; 
for a large mismatch angle it is reflected basically by 
the coarse control signal. In order to avoid reverse de- 
flection of the null indicator pointer when the phase of 
the fine control signal changes [i.e., when the value of 
the mismatch angle is in the range */k-20V/k (i= 1,2, 
3,...)], a branch with the resistor r, is introduced into 
the network. The resistor rs applies a voltage to the non- 
linear resistance for maintaining the phase of the mis- 
match, 


The Dependence of the Deflection Angle 
of the Null Indicator Pointer on the Mis- 
match Angle 











The dependence of the deflection angle a of the 
null indicator pointer on the mismatch angle @ can be 
derived by making use of the mesh current method, In 
order to facilitate the derivation of this relationship,we 
neglect the inductive reactance of the null indicator 
signal winding, since it is small compared to the resist- 
ances of the loops; we shall assume that the shape of the 
current sinusoid (at the carrier frequency w) is not dis- 
torted by the mismatch angle. Then 


Ug = 1, (ra tera + rat rs) + La (ria + 1) — La (rin + 71) + [7 
Us = 1, (ra +11) + Ie (ra +i +7) — Sara +7), 


U = Trg + 1s (rig + 74 + 7s) + Tria, 


U.—U; = — 1, (rig +171) — La (ria + 71) + Loria + La (rin + Pte +78 +71) 


where r is the static resistance of the NR and is deter- 
mined from the volt-ampere characteristic; tj, fj, are 








the resistances of the output windings in the detector 
elements, 
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Fig. 1. 
The determinant for the system of equations is 
A C rg —C 
C BOC 
A= ee D te = AB (DF — r3) — BFr: — 
— Cc —_ Cc Tig F 


—C*(AD + DF + BD — r2 — 72, — 2rsria) + 2DC* — 2Brzrinl, 


where 
A=ratn+tret+7s, B=ry¢+nit+r, 
C=ry+1, D = rig + re +a, F=ry t+ rig t+r5 +73. 


The current flowing through the signal winding of the null indicator is equal to 1=1,+Is, where 1,=4,/ A, I-A, /A, 
Us Cc rs — Cc | 
a Us BO —C|_ 
7% Uc 0 D Tig nt 
Uz — Us —_— Cc Tig F 
= U, (BDF — Br:, — CDF +- Crgrig + Cr?, — BDC + C*D — Brgrig) +- 
+ Uc (— Br3F — BCrig + C*rig + C*rg + BDE —C*D + Brzgris), 


no eee 
oy er ae er aa 
pti rs 0 U-. Tie a5 


—C—C Uce—Uf F 
= Us (ABriz — BCriz -+ C*ris oo BFr; — ACrig + CFr; — C*r, aed BCrs) oo 
+ U.(ABF + 208 — C*B — AC* — C?F — ABri, + r3C*? — BCrs + C*rig). 
Turning our attention to the deflection angle of the null indicator pointer instead of the current, we obtain 


a= SJ, or sn ig OO Os © 
Bm, + No ’ 





where 
m = DF — rz, — DC — rsrig + Arig — Crig — Frs + Crs, 


N =—_l oe CDF +- CrsFie 4. Cr, a C2D +- C*ri. — ACTiz +. CFr, —— C*rs, 
n= -— r;i* _ Cris 4- DC -|- Tsli2 -!- AF — Cc? — Arig — Crs, 

N, = C*rig + C*rs — C2D + 203 — AC?— OF +. 7,C? + Cre, 

mg = ADF — Ar? — Fr? —C?D — 2ryrinl, 

Nz = —C?AD — C2DF + C2r2 + Crd + 2DC? + 2rgrigC. 


Here mj, Nj are constant coefficients (i= 0, 1, 2). 
The expression for the coefficient B includes r (the static resistance of the nonlinear element). 
After replacing B with the expression Ry + r (where Ry = 14, + 1) and replacing U., and Uf with the expres- 


sions U,, sin 9 and U;, sink®, we obtain 
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This formula describes the dependence of the de- 
flection angle of the pointer on the mismatch angle be- 
tween the receiving and transmitting devices. The for- 
mula is inconvenient because it includes three variables: 
the mismatch angle @, the deflection angle a of the 
pointer, and the static resistance r of the nonlinear ele- 
ment, 

In formulating the functiona=f(@) we shall as- 
sume that the shape of the current sinusoid flowing 
through the movable winding of the null indicator is not 
distorted by the variation of the mismatch angle (or by 





sin k® [(2y + r) m + N) + sin ® [(2y +r) my + Ny) 
(i +r) me + Ne , 


(2) 


the variation in the position of the operating point on 
the volt-ampere characteristic of the nonlinear element, 

In order to eliminate one of the variables from the 
latter formula, we must determine the function r = f,(@) 
orr = f,(a), Obtaining either of these functions in ex- 
plicit form presents a complex mathematical problem. 
However, it is not mandatory to have these functions in 
explicit form to formulate the function a = f(@). In or- 
der to find r=f,(@),we shall make use of a function that 
expresses the current I, through the nonlinear resistance 
in terms of the mismatch angle 6: 


— 4s 
I, = nv 
where 
A Us rs - C 
pom Cc Us 0 —Ci|_ 
st rs Ue D Tia 
—C &—Uf Tig F 
= Us(— CDF + Cr. + ADF — Ar}, — Fr +C*D+ crs — ADC) + 
+ Ug (C*rig — C*rs — ACrig + CF rs — C*D — Cr? + ADC — Crorig) = 
= U;Q ao U,w. 
Then 


U;Q “f Uw 


a= (Ay +r) ms + N° 


(2) 





As we know, the relationship between the current 
and voltage for nonlinear solid-state rectifiers is given 
by the formula 

By 
) 


where U,, ig are the instantaneous values of the forward 
voltage and current for the rectifier; Ugy, Ig, are the 
forward voltage and current for the rectifier when its 
class is specified [1]. 

Turning from the voltage to the static resistance of 
the rectifier in the forward direction, we obtain 


ta 
1 


Ue= Ven ( 





U, U., i& aie 

ly Madkec ee a 
where es Van 
Cand” 





In the case under study, ne , where P=8,-1, 
Substituting into Formula (2), we obtain 


Vx - U,0+tw 
a 


(Ai +r) me + My 





or 
P P 
WV r+ lV r =Um(QsinkQ + wsin9), 
where 
= i, »t= 7 (Rum + Ns). 
VM VM 


Specifying the value of r, we find the value of 6; 
then these two values are used in conjunction with 
Formula (1) to determine the deflection angle a of the 


pointer, 
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In Fig. 2 we have shown the dependence of the de- 
flection angle a on the mismatch angle @ for two dif- 
ferent cases (for a closed circuit and an open circuit at 
the point a in Fig. 1). 


A Simplified Method for Computing the 
Network 





In performing the computation (selecting the param- 


eters),we assume the following parameters are known: 

1, The maximum voltages applied to the detector 
elements, and the transfer number between detector ele- 
ments, 

2. The null indicator sensitivity and the resistance 
ts of the signal winding. 

3, The volt-ampere characteristic of the nonlinear 
element, 

4, The resistances of the output windings in the de- 
tector elements, 

5. The mismatch angle @ within whose limits the 
null indicator pointer is deflected predominantly by the 
fine control signal. 

6, The mismatch angle 6, - 6,, within whose lim- 
its the pointer is deflected predominantly by the coarse 
control signal, 

1, The deflection angles a, and a, of the pointer 
which correspond to the angles 6, and 63. 

The angle 6, must be equal to the deflection angle 
of the pointer up to the pointer stop. 

The voltage drop AU across the nonlinear resist- 
ance is made equal to,or greater than,the null indicator 
voltage drop corresponding to a pointer deflection angle 
Os. 





where 


The resistance r, is computed from the formula 
tf, = 50 —ty The derivation of the formulas for de- 
1 
termining the resistors r; and r, is achieved on the basis 
of the following conditions: 


U , sin k8, — AU U,, sin 8, — AU AU 





n+l ra + MW ‘eq a 
U,,. sin U 
os oa 2 , (4) 
st "ts retry 


Formula (3) is the condition for deflection of the 
pointer through the angle a, when the mismatch angle 
is @, (excluding the effect of the current I,), Formula 
(4) is the condition for elimination of the possibility of 
reverse deflection of the pointer when the fine-control 
detector element is rotated through approximately 270° 
from the matched position. 

In Formulas (3) and (4), 


9 ae) . (rs +42) r’ 
—" em ° CQ UT Ptet+fst+r’ ’ 
where r’ is the resistance of the nonlinear element at 
the knee of the volt- ampere characteristic, 
After a series of simple transformations we obtain 


%) 2B mm — 8 ary, + Ae B, 
= ed += Ai — Fig, 








A, = Umsinké, — AU, B, = U,,sin6, — AU, C, = Py $+ Tia, 


AU 


e = U,sin6,, d = -——. 


The resistance r, is obtained from the formula 


U,,8in 0, — U, sin 01 
a“ Qs — a { 





rs 


The network described above was checked under 
laboratory conditions, Here the network parameters 

















a a 
a, — 
°, a, 
Fig. 2. 


eq 





were as follows: r,= 100 kohm, t,=1 kohm, r= 20 kohm, 
t5= 10 kohm. 


The null indicator was a ferrodynamic instrument of 
the type “IF-2"; the nonlinear resistance was a selenium 
bar of the type “AVS-18-18" with shorted outer leads 
(the voltage drop across it was ~1,1 v), The signals 
were generated by a receiving device with "SS-405" 
selsyns, The transfer number between the selsyns was 
k=20. The deflection angle of the null indicator pointer 
was &,=~25° and corresponded to a mismatch angle of 
~3°; the deflection angle a,=45° corresponded to the 
mismatch angle @,= 12°, 

Thus our network for summing the signals permits 
automatic assurance of a high sensitivity to small mis- 
match angles (this is the function of a fine-control null 
indicator) and a low sensitivity to large mismatch angles 
(this is the part played by a coarse-control null indicatod 








The fact that it is not necessary for the operator to 
shift his eyes from one null indicator to the other per- 
mits an improvement in the quality of semiautomatic 
control and makes it more convenient to monitor the 
accuracy of the automatic tracking. In addition, the 
construction of the mounting of the device is simplified, 


and more advantageous conditions are created for its 
protection from shock. 
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The paper performs a quantitative investigation of the process involved in switching the mag- 
netic state of sufficiently thin plane cores with rectangular hysteresis loops (for example, mag- 
netic ribbons with a thickness of up to 3-5), For thicker ribbons the qualitative picture is 
described, A criterion is established for the intensity of the magnetic skin effect; this criterion 
is a relative increase in the switching coefficient for the core, The results of the computations 





are compared with known experimental data on the magnetic switching of ribbon cores, 


In various fields of automation and computing engi- 
neering magnetic elements in the form of ribbon and 
sheet cores are widely used. Such cores are used, for 
example, in registers, in logic and memory cells, in de- 
coders, etc, It is of interest to analyze the physical proc- 
esses and the electrical computation method for the vari- 
ous types of pulsed magnetic circuits while taking into 
account the effect of certain basic factors (the magneto- 
viscous properties of the core material, and eddy cur- 
rents) which lead to closer coincidence between the com- 
puted results and the observed phenomena. 

The problem of studying the effect of eddy currents 
for switching the magnetic state of cores with a rectan- 
gular hysteresis loop requires the solution of nonlinear 
differential equations. 

For particular cases involving ribbon cores we use 
numerical integration methods [1]. The results are in 
good agreement with experimental data, In this paper 
we obtain an analytical solution of the problem. In gen- 
eral form such a solution may be obtained for the case 








jj 


Fig. 1. Diagram illus- 
trating the problem. 





of a weak surface effect which is of special practical in- 
terest, 

The presence of the surface effect leads to a pro- 
longation of the magnetic switching process and to a 
change in its form. We derive an expression which is 
convenient for determining the prolongation time. The 
paper shows the dependence of this time on thickness and 
properties of the ribbon core material. 

The problem is solved by the method of successive 
approximations used earlier for an analysis of the surface 
effect in ferrite ring cores. We shall cite the basic as- 
sumptions of this method and shall apply them to cores 
shaped in the form of a plate (Fig. 1). At the initial in- 
stant t= 0 the plate is magnetized in the direction of the: 
y axis to an amount equal to the remanent flux density 
B,. The magnetic switching is performed by an external 
magnetic field H= He(t) which is applied along the y 
axis, a 

The magnetic switching process in the surface layers 
of the plate is determined exclusively by the magneto- 
viscous properties of the ferromagnetic. The change in 
magnetic state can be expressed in terms of a pulse Q of 
the switching magnetic field [2): 


‘ 

Q= \(H—H,)at, (al) 
0 

B= B,th (ze —9). (3) 
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Here B, is the saturation flux density; r,,, and Hg are 
experimentally determined parameters that characterize 
the dynamic magnetic properties of the core material; 

B 
= Arth a * ; 
Y=Ar By 
For the surface layers (H=He, Q=Qe), 


Q, = \(H.— Ha)dt. 


0 


The graphs corresponding to Expressions (2) and (3) 
are plotted in Fig. 2,where the quantities B and Q, are 
given in relative units. They describe the magnetic 
switching process B=Be(t) in the surface layers of the 
plate exactly, Analogous assumptions are also the basis 
of the computation performed in the widely known paper 
[3]; however, the authors of this paper made use of a 
coarser approximation of expression (2) in analyzing the 
surface effect; 0B/8Q=const*. Such an assumption can 
be justified for a strong surface effect (in cores with a 
ribbon thickness of the order of 100 microns). 

The magnetic switching in the inner layers is relat- 
ed to the external magnetic field pulse in a more com- 
plex fashion, It is necessary to take into account the 
screening effect of the eddy currents, The eddy currents 
induced during the magnetic switching of the plate pro- 
duce an additional pulse Qjy for the inner layers, and 
thus ,Q=Q, + Qin. 

In order to take the eddy currents into account we 
make use of the integral form of the electrodynamics 
equations, Assuming that in a first approximation the 
value of the flux density B is constant over the cross sec- 
tion for the case of a weak surface effect, we proceed to 
solve the problem in the following approximation. From 
the first Maxwell equation we obtain the first approxima- 
tion for the intensity of the induced electric field: 


0B 
E=2£>5-. 





Fig. 2. Basic graphs; 
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The induced field is directed along the z axis. It 
causes eddy currents whose density J (taking the dis- 
placement currents into account) is determined by the 
electrical properties of the material; the conductivity 0 
and the permittivity €: 

0 \ 6B 
J =x(0 -f- 0 3) ar 

The magnetic field is computed from the total-cur- 
rent law, while taking the boundary conditions (H= He) 
at the surface of the plate into account: 


H=H, — + (bh —24)(c + e& a) ar 


where b is one-half the thickness of the plate whose mag- 
netic state is switched, 

Thus,in a first approximation the magnetic switch- 
ing pulse produced by the field depends on the depth of 
the layer in the plate, The magnitude of the pulse can 
be derived on the basis of (1) by integrating the expres- 
sion above for H with respect to t: 


The correction for the first approximation consists 
of a correction proportional to the conductivity of the 
ferromagnetic; 


= — + (P—2(0 + ee, 5)(B + B,). (4) 


n 


Qo =—+-(—2o(B+B), ©) 


and a correction which depends on the permittivity. 
The latter component can be written in the following 
form if we take (2) into account; 


Bt 


Q.-— ap eaneli— 2 


)(H — Ho) 





(6) 


* In terms of the notation used in [3], 8B/dQ = f(B). 
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Thus,both corrections to the pulse are expressed as 
a function of the flux density B. Expression (5) in Fig.2 
corresponds to the straight line 2 whose angle with the 
flux density axis is proportional to the conductivity of 
the ferromagnetic ,and increases as the layer under study 
is further removed from the surface of the plate. For 
the surface layer the angle which the straight line 2 
makes with the B axis is equal to zero (x?=b*). For the 
axial layer of the plate it is a maximum (x*= 0). 


The graph of the function (6) is a quadratic parabola 
(Curve 3 in Fig. 2) which intersects the axis of ordinates 
at the points +1, The maximum value of this correction 
increases with the intensity of the magnetic switching 
field. An interesting feature of this correction is its rap- 
id tendency toward zero at the beginning and end of the 
magnetic switching process, This explains the fact that 
the effect of the permittivity on the time required for 
completing the magnetic switching process is weakened. 


The simplicity of expressions (5) and (6) permits 
easy plotting of the corresponding graphs for the chosen 
value of x. By the same token it is easy to plot the 
graph of the total correction (4) for a definite value of 
x (Curve 4 in Fig. 2). 

An important problem is the appropriate choice of 
the values of flux density in Formulas (4), (5), and (6). 
The simplest approach is to choose this value equal to 
the magnitude of the surface flux density B=Be. 


Since the magnitude of the flux density actually de- 
creases with the depth of the layer, such a computation 
in a first approximation yields values of flux density that 
are too low. The computation will be somewhat more 
exact if we choose a value of flux density equal to the 
flux density in the axial layer for use in Formulas (4), 
(5), and (6). Here the result obtained in the first approx- 
imation will yield values of flux density that are too 
high for the individual layers. 

We chose a somewhat improved form of the approx- 
imation : 


B=B,th (le, + Q(B) — o| 


which also resulted in values that were too high, In ad- 
dition to achieving a certain refinement, the approxima- 
tion form above simplifies the subsequent analysis sub- 
stantially. The expression shows that in order to plot 

the graph for the variation of flux density in a layer with 
a specified value of x in the approximation under study, 
it is sufficient to shift the points on the original curve: 


B, = B,th (e.—¢) 


along the Q, axis by an amount equal to the correspond- 
ing value of Qin (Fig. 2). 





The error of the approximation can easily be esti- 









: BL—B\2 
mated. It does not exceed the value 50 ( : ) %, 
B, are B, 





where the values of the flux densities B and Be are taken 
fot one value of time (Fig. 2). 

Figure 3 shows the corresponding computation for 
the case of surface effect in a thin Permalloy sheet of 
the type "79-NM”" which was 5 p thick, [See Appendix 
for "79-NM” characteristics,] In that case the surface 
effect is sufficiently weak, and the approximation under 
study assures an accuracy of approximately 10%. 

The analysis of the surface effect is simplified for 
the case of metal ferromagnetics in which the displace- 
ment currents can be neglectedscompared to the conduc- 
tivity currents,right up to frequencies in the optical spec- 
trum [4]. In that case the correction of the pulse is de- 
termined basically by the linear component of (5). This 
case also includes the example of a Permalloy plate 
cited above. The magnetic switching curves for differ- 
ent layers in a metallic ferromagnetic (as obtained in 
the approximation under study) have an interesting prop- 
erty. The tangents to these curves always intersect at 
one point lying on the straight line B= ~—B, (Fig. 3) when 
they are drawn to points with an arbitrary fixed value of 
flux density. A definite point corresponds to each value 
of flux density B(-B,=B=B,). This property makes it 
possible to compute the magnetic switching rate 6 B,/ 
/0Q, for a layer with the coordinate x in terms of the 
known function of the magnetic switching rate in the 
surface layer 0B,/0Q¢. 

For a flux density By in the layer x, 








aB, = eB, »— 2s \7 
w, ~9e (tag FH2) 


where the quantity 0Be /OQ, is determined by substitut- 
ing Be =B, into expression (2), 

In order to compare the computations above with 
experimental data, it is necessary to make the transition 
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Fig. 3. Surface effect in a Permalloy ribbon. A) Mag- 


netic switching of the surface layer; B)magnetic 
switching of the central layer. 
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to the experimentally observed average values of flux 
density, since the emf in the metering winding is pro- 
portional to the rate of change of the average flux den- 
sity. In the approximation under study we have 


OB. r 
av_} — m 
max (3) = iene ’ (8) 
m°"3_ 


with an accuracy of no less than 12.5 (tmob*)"%;i,e., the 
maximum rate of change of the average flux density is 





6B r 
max ( art ) = -—— (He. — Hy). (8a) 
1+ mS 7 
The maximum value of the magnetic switching rate 
is achieved for an average core flux density 


r B { 
Bay.m= ae nf * =? (9) 
1 + 49 a 





The approximate expressions (8) and (9) can be used 
for the independent determination of the parameter rp 
from the experimental data for (@Bgy/0Q¢) or Bay, m 
that are obtained when a plate made of a known metal- 
lic magnetic material is subjected to magnetic switching, 

In order to check expressions (8) and (9) we make 
use of the results obtained from an investigation of rib- 
bon cores made of Permalloy "79-NM"* ribbon 10 and 
20 p thick [1], The surface effect in such cores is quite 
pronounced, Nevertheless, as we shall demonstrate be- 
low, Formulas (8) and (9) are in good agreement with 
experimental results, We shall write the known charac- 
teristics for "79- NM" Permalloy: 


B, = 7000 gauss; B= 4900 gauss; o= 1.8 10* 1/ohm -cm. 


In order to select the value of r,, we shall make use 
of the results obtained by V. L. Dyatlov, who performed 





an exact numerical integration for the problem. In [1) 
the value r,, = 950 ohm/cm was obtained, The results of 
the experiments and computations are cited in the table, 
In conclusion we shall estimate the most important 
characteristic of ribbon cores —the duration of the mag- 
netic switching~in the approximation used above, When 
there is no surface effect,this duration is determined 
completely by the magnetoviscous properties of the core. 
The theoretical duration of the magnetic switching given 
by expression (3) is infinite. Therefore in practice,we 
use the concept of the duration required to reach a 
steady state at some level, We often choose the average 
value of the flux density over the cross section Bay =B, 
as this level, The time t, required for the magnetic 
switching process to reach this level is assumed to be the 
duration of the process, Expression (3) permits easy de- 
termination of the field pulse Q, corresponding to this 
time (on the assumption that there is no surface effect); 
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When the surface effect is present.the magnetic 
switching process in the core terminates at the level 
Bay =B, for a field pulse which is increased by the 
amount —Qin(B,): 


Qn — —g P*B,. 


Thus, the relative increase in the pulse required for 
reaching completion of the magnetic switching at the 
specified level is 





IQ, pe 2B, "mn? 
Q 3 B @ ° “ 


Expression (10) determines the relative increase in 
the switching coefficient for the core. For magnetic 
switching of a core by means of rectangular current 
pulses (here it is possible to assume that He ~ Hp* const), 
relationship (10) also determines the relative prolonga- 
tion At, of the pulse, In that case, 
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Expression (11) shows that the relative prolongation 
of the switching time due to eddy currents is in a first 
approximation independent of the magnetic switching 
rate and directly proportional to the square of the plate 
thickness, 

For the case of ribbon cores with the parameters 
cited above, computation according to Formula (11) leads 
to the following values for the relative prolongation of 
the magnetic switching time: 28% for a thickness of 5 p; 
110% for a thickness of 10 p. 

These results can be used in designing magnetic ele- 
ments for computer- decision units, 

Expression (11) permits us to compute the allowable 
thickness of ribbon cores which assure the maximum 
rate of recording and readout of information. 

The recommended plotting of the graph for the de- 
pendence of the average magnetic flux density in the 
core on the magnetic switching pulse can be made the 
basis of the electrical design for pulse networks with rib- 
bon cores (for example, magnetic registers), 


APPENDIX 

We shall prove certain propositions cited in the 
paper. 

1, We shall derive Formula (7). Assume that the 
curve B. =Be(Qe) is given. We shall study a new curve 
B =B(Q,) which is obtained from the first curve by shift- 
ing its points along the axis O, by the amount f= {(B). 
Here f(B) is a specified function of B. The curve B(Qe) 
is obtained by replacing the argument Q, by the argu- 
ment Q=Qye - f(B) in the curve Be(Q,): 


B(Q,) = B, 12, —f (B)}. 


Differentiating this equation, we obtain 
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0Q, is 9(Q,—7()) 0Q, x 
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From the equation above we obtain 
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0Q, ee: 8Q, oB 9Q, 
where the derivatives are taken for the identical fixed 


b? - x? 





value B=B,. The substitution f(B)=-Qo= 
X 0 (B+B;) leads to (7). 
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2. We shall show that for a weak surface effect the 
the flux density , averaged over the cross section of 
the coresis approximately equal to the flux density in a 
certain layer—the layer located at a distance x=b/V3 
from the axis. In a linear approximation, 


dB 
B (Q2) = B(Q:) + (Q2— Q1) 90’ 


where Q, and Q, are the values of the external field 
pulse at any instants which are sufficiently close to- 
gether. Assuming Q1=Qe~ Qin and Q,=Q,, we arrive 
at a convenient expression which relates the values of 
the flux density in the surface layer and in the layer cor- 
responding to the specified value of the eddy current 
field pulse: 

aB, aB, 
B(Q,) = B(Q, — G+ in a0, ~ B,(Q.) + Gn 2, 

(12) 


Equation (12) can easily be averaged if we keep in 
mind the face that only the term Qj, , which is propor- 
tional to (b?~x*)/2 according to (4), is dependent on x: 


dz = 





B 
iT OB, 
Kay (Q,) = + B, (Q,) + QndQ, 
0 
= B,(Q,)— = (o + et, =) (B, + B,). 
(13) 
Expression (13) is equal to (12) for Qin = Qin (x) 
wheu x=b/V¥3; this is what we set out to prove. 
Assuming B.=0 and Qin =Qo in (13), we arrive at 
the approximate expression (9). 


TRANSLATOR'S APPENDIX 

The characteristics of the permalloy "79-NM”" are 
the following: thickness 0.01 mm; initial permeability 
H» (in gauss/oe) not less than 15,000; maximum perme- 
ability Upax (in gauss/oe) not less than 65,000; coer- 
cive force H, (in oe) 0,07; saturation flux density Bs 
(in gauss) not less than 7,500. (See also "Characteristics 
of Magneto Materials for Use in Magnetic Amplifiers 
Operating at Audio-frequencies," Automation and Re- 
mote Control, Vol, 18, No, 10, pp, 927-933.) 
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This article presents the theory of photoelectric fluxmeters (PEF), i.e., amplifiers with photo- 
electric optical amplification, Various types of PEF with and without feedback are compared, 
It is shown that the greatest accuracy of integration and the greatest speed can be achieved 
with a PEF circuit with negative feedback with respect to the derivative with the smallest 
time constants in the open circuit, The proposed theory makes it possible to determine the 
integration accuracy in dependence on the signal shape directly with respect to the curves 


recorded at the output, 


At the present time, the so-called photoelectric op- 
tical amplifiers (PEOA)* [1] are being used to an ever 
increasing extent as dc amplifiers in various measuring 
and automation circuits, 

The principle of photoelectrical optical amplifica- 
tion consists in using the galvanometer of the magneto- 
electric system at the input and in converting the dis- 
placement of the galvanometer light spot into voltage 
by means of an optical system and a photocell bridge [2]. 

The extremely low input noise level (~10°°-107**y) 
and the large amplification factor (up to 10° for voltage) 
make it possible to use PEOA in devices where ordinary 
electronic amplifiers cannot be utilized. 

A particular example of such a device is the photo- 
electric fluxmeter (PEF) [3]-an integrating amplifier 
based on PEOA-for measuring magnetic fluxes or elec- 
tric charges. The considerably greater sensitivity and 
accuracy, as well as a longer allowable measurement 
time in comparison with ordinary fluxmeters, render the 
PEF a valuable instrument in scientific-research work 
and automation. 

The face that different PEF circuits are recommend- 
ed in literature (see survey [3]) and a lack of compara- 
tive evaluations of their efficiency create a pressing prob- 
lem, which is to be solved by theoretical analysis of the 
existing variants for the purpose of selecting the variant 
most suitable for practical application. 

The PEF theory presented in this article is a develop- 
ment of the method described in [4], which is based on 
operational calculus [5}, 


PEF Circuits and Their Relation to Inte- 





grating Amplifier Circuits 
The block diagrams of PEF are identical with the 
basic integrating amplifier circuits [6], which have been 





developed independently and mainly before the PEF cir- 
cuitst . 

Figure 1 shows an integrator circuit with a negative 
series feedback with respect to the derivative. This inte- 
grator pertains to the overwhelming majority of the PEF 
circuits described in literature. The wiring diagram of 
such PEF is given in Fig. 2 and has been described in [3] 
and [7). The time constants of the amplification circuit 
(without taking into account the feedback) of this PEF 
variant are rather small [8], and, in the first approxima- 
tion, this circuit can be considered as an ideal amplifier 
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Fig. 1. Circuit of the RC-integrator 
with a negative series feedback. 


* There are now greater opportunities for the application 
of PEOA, since the industry has started the mass produc- 
tion of amplifiers of this kind (F17 and F117). 

T The actual differentiating circuit in the feedback loop 
is of no importance, For the sake of clarity, we shall 
use an RC circuit throughout, 

The RC, circuit mainly serves for ensuring sta- 
bility. 

The "+* and *~ * signs in Figs. 1-5 conventionally 
indicate the presence or absence of a 180 deg phase shift 
of the output voltage with respect to the input voltage in 
the amplification circuit of the open-loop system. 








with a constant amplification factor [9]. In [3] and (4), 
PEF of this type were called PEF of the first group. 

The prototype of a PEF without feedback, which has 
been proposed by S. P. Kapitsa [10], has an integrating 
amplifier (Fig. 3). This device uses the principle of inte- 
gration by means of a RC, circuit with a decoupling 
amplifier A, where the usual resistance is replaced by a 
higher-the dynamic resistance R of the amplifier A, out- 
put tube, The introduction of a dynamic resistance in- 
creases the time constant T,=RC;, of the integrating cir- 
cuit and brings the operating conditions of the circuit 











section, which is the input section with respect to the C; 
capacitor, closer to the operating conditions of the cur- 
rent generator, In Kapitsa's PEF, PEOA plays the role of 
the A, amplifier, and the dynamic resistance of vacuum 
photocells (p * 10° ohm, Tp* pC, * 1000 sec) plays the 
role of the tube dynamic resistance. 

The application of negative feedback with respect 
to the derivative in Kapitsa's model leads to a circuit 
which is identical to circuits of PEF of the first group 
(Fig. 2), but is characterized by a larger capacitance Cy 
(almost by three orders of magnitude), 
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Fig. 2. PEF wiring diagram which is similar to that in Fig. 1. G) Galvanometer; L) il- 
luminator which projects the light spot image on the photocells(PC) or the photoresis- 
tors (PR); A) electronic power amplifier; ©) magnetic flux to be measured. 


The integrating regime in the amplification circuit 
can also be secured by other means [11]. According to 
[3] and [4], we shall designate PEF with feedback and 
integration in the amplification circuit as fluxmeters of 
the second group. 

PEF circuits with combined feedback have been 
described in a number of papers [3]. The circuit in 
Fig. 2 (the block diagram is shown in Fig. 1) is supple- 
mented with rigid positive feedback with the coefficient 
Ag=t/ (ty + ts) (Fig. 4). 

The most widely used circuit for analog computer 
integrators is the circuit with parallel feedback (Fig. 5). 
By using Kirchhoff*s laws, it can be readily shown that, 
in contrast to the circuit in Fig. 1, the feedback of Fig.5 
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Fig. 3. Circuit of a RC-integrator using the dynamic 
resistance R of the amplifier A, output stage. 





does not reduce the input load current. Since this cur- 
rent lowers the magnetic measurement accuracy, the cir- 
cuit given in Fig. 5 is not used in PEF designs, 

The circuit shown in Fig. 2, supplemented with posi- 
tive feedback in the general case (Fig. 4), will be used 
as the basis hereafter. All particular circuit variants can 
be obtained by eliminating some of the circuit elementst 
or by changing the numerical value of individual param- 
eters, 


Basic Approximations of the PEF Theory 





In the proposed PEF theory, it is assumed that all 
directional elements,as well as amplitude distortions and 
loads,are absent. In particular, it is assumed that the 
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Fig.4. Wiring diagram of a RC-in- 
tegrator with combined feedback. 





3+ 





+e——_—_ 


























Yout 








+ Thus, for instance, Kapitsa's circuit without feedback 
is obtained by eliminating the feedback circuit. 
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connection of a measuring coil and an output instrument 
does not affect the PEF operation. 

The conditions which the parameters of the loads 
to be connected must satisfy if the presented theory is 
to remain valid,as well as the justification of the other 
assumptions, will have to be treated in a separate paper 


Notation and the Basic Equations 











The practical system of units is used throughout. 

The electromotive force Uj, induced in the meas- 
uring circuit or any voltage supplied to the instrument 
input is taken as the PEF input. The output voltage U 
of the amplifier A is at the same time the PEF output; 
f is the difference of light fluxes incident on the photo- 
cells of the photocell bridge; v is the bridge diagonal 
voltage; v, is the feedback voltage; vg is the same volt- 
age without feedback and without taking into account 
the capacitor C, leakage; t and @ are time; s is the 
differentiation operator. 

The operator conductances (transfer functions) of the 
elements are given below. 

1, Mirror galvanometer with an optical lever (see 
[11] and [12]). If we take into account the inductance 
of the measuring circuit, 


NS 


(tos + 1) (Hs* + Gs + 1) ' 





v,0= f=W,Yin (1) 


In practice, the T», H, and G parameters are close 
to the measuring circuit time constant (T,)~ L/R, where 
L is the inductance, and R is the resistance), to 1/ug, and 
to 2B/w) (wy is the galvanometer natural frequency, and 
8 is the damping coefficient), respectively. Moreover, 
in Eq. (1), N is the optical system conversion factor in 
lumeus per radian; Sy is the galvanometer sensitivity in 
radians per volt for the chosen operating conditions. 

As a rule, the galvanometer is overdamped* ® , and 
the measuring circuit inductance is neglected. Then 





Ww, (s) = is ’ 
(Tis + 1) (Tes + 1) (1a) 
Tit, = A, M+ =G, T1>tTs* 
R c 














Fig.5. Wiring diagram of a RC- 
integrator with parallel negative 
feedback. 


2. Photocell bridge (vacuum photocells FE or photo- 
resistors FS [7]). A stabilizing (correcting) Ry,C;, circuit 
is connected to the bridge circuit; 


Ts +1 


Wm (8) = OP Fo 4° v=W,,f. (2) 


Here, a is the FE sensitivity (amp/lu), p is their 
dynamic resistance, 


T= % + Ts, ™4= R,C,, ™3&= pc,. 


3. Direct current amplifier (the amplifier inertness 
is neglected); 


W, =k = const, U = W, ». (3) 


4. Differentiating element (without leakage and 
without taking into account the load [13)): 


Tas 
Wa(s)= tert 1,=W,U, (4) 
where Tq=RyCy, Ry is the circuit resistance, and C, is 
its capacitance, 
5. The feedback operator (Fig. 4): 


R 
W; (8) = Wa + G-— dom Wa—A. (5) 


where Rg, is the resistance of the C4 capacitor leakage, 
and Ag=r9/(t, +19) is the positive feedback coefficient, 
Equation (5) follows from Fig. 4. 

The open-circuit operator (without feedback) is ob- 
tained from (la), (2), and (3): 


Wols)=W WW =Kem(s),  U = Wellin (6) 
Here, 


Ky = S,Napk = const, 


e+ (7) 


Wo (8) = is +1) (ts $1) (tos +1) 





We shall introduce the PEF constants Cy, and Cy so 


that, with an accuracy to the basic errors , j Vi, 40 = C0 


for PEF without feedback (including PEF of the first group 
With an open cireuit, for which the integration error ex- 


ceeds Cy U)s also { U,, 40 = C,U when feedback is 
0 


** This is mainly connected with the requirement for 
greater vibration stability of the instrument during zero 
adjustment (with the Cq capacitor short-circuited), 
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present, 

We shall also introduce the creep time constants T 
T and T, for PEF with feedback and T, for PEF without 
feedback, It will become clear later on that 





T = K,C., Ci =%, + C.. (8) 
T 
T,=T—AK,’ (9) 
Ty = 13 + Ta (10) 
T 
Cy, = Ke: (11) 


By T,, we shall denote the time interval during 
which integration is performed. 


PEF Theory 
A. PEF without Feedback 
In this case, U=WgUj,. On the basis of (6) and (7), 








Ky Ts +1 1 
Tos + itis + 1T8+1° 





Wo (s) = (12) 
If the galvanometer is overdamped, it can be con- 

sidered that r,=0 in the first approximation. The cir- 

cuit parameters Ry and Cj, can be so chosen that [10] 


% = T}. (13) 


The first factor, to which the expression for Wo(s) is 
now reduced, coincides with the operator of the integrat- 
ing RC; circuit with an amplifier (cf. Fig. 3; in the given 
case, R= p + Ry): 


K 
Wols) = a547° (14) 
Hence [6], 
K t 
u =| Gnas (15) 


with the error in percentages [6] 


T, 
p= 50 . (16) 

Here, Ts is the creep time constant. 

For a more complete theory, we shall seek the rela- 
tion between the observed output voltage U and the volt- 
age Ujq of an ideal integrator of the same sensitivity as 
that of the actual one under consideration. For this pur- 
pose, we shall multiply every term of the identity 


ee oe egg | 
Te We Fat Faw! an 


by U after the preliminary transformations? t 


<= (¢-1)U = +0 + 





+{ B,U” + B,U’ + BU 


Tl m Bu’ |= Tw +P 


or 


i/i 
+ (G1) 0 = (es + 1) U + BU" + (Bs + BYU’ + 


1 
+ 10) =—TU + P. (18) 


Here, P is the portion of expression (18) in square 
brackets, 


By = Tats (Ti — Ta), By = (Ta + Ta) (T: — Ta), 


Bs = 1 — Ta, By = TaTs- 


As a result, we shall find that the ideal integrator 
output voltage is 


t 
Ke nae 
0 
Ug = Te\ Up O=U4-- +7. OM 
0 





The second term in (19) is the creep, and the third 
term represents the error which is due to the fact that 
the second and the third factor in (12) differ from unity 
("structural® error). 


B. PEF with Series Feedback 

A circuit with combined feedback represents the 
most general case. By substituting operators for figures 
in the well-known expression for the amplification fac- 
tor of an amplifier with feedback, we obtain 





Wo 


W (s) = i+ Ww, , 


U = Win. (20) 


TT Creep is the instrument reading deviation at any in- 
stant of time from the readings of an ideal integrator; it 
is determined by the finite magnitude of the largest time 
constant in the system, This definition is a natural ex~- 
tension to dynamic operating conditions of the concept 
connected with the inability of the needle of ordinary 
fluxmeters to remain in indifferent equilibrium if no in- 
put voltage is present, 

++ Transformation after the * sign is valid only if ™%<<1. 
For the evaluation of the error due to the ™%s+1~*1 substi- 
tution, see [4]. 
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By using the identity 


aa 
T4s Ww (21) 


{ ‘ 4 1 Ww 
mi+- - ~ —- 
+ s | Kota ' ( Waxy Kota ) + ( Ta ‘ )] : 








we obtain,in a manner similar to that used in deriving 
(19) after a number of transformations [including (18)], 
{ c 
Vig = Cc \ Uin dQ = U +f Er, oa Er, + Ers, (22) 


0 





where 
t 
Ud§ 
Pea ; (23) 
4 1 -—— Te 
is the creep, and 
P 
Er,  — T (24) 


is the error due to the difference between the amplifica- 
tion circuit and the ideal circuit (when W»= 1); this error 
will be called the “structural” error. The quantity 


t 
Ta 


Wa 
Ers = ate (UU) a0 — eu" (25) 
0 


is the error caused by the differentiation error in the 
feedback circuit; this will be called the “differential” 
errors 

C. Effect of Disturbances and Distortions 

Let a disturbance, equivalent to the Up voltage, 
arise at the output of the element (circuit with feedback) 
whose output signal U, is related to the input signal by 
the equation 





By finding the equation relating Up and U (for Ujn= 
=0), we obtain the instrument output error 
Ws 


Er, =-UJ0= T+WwW, Uy, W\W,= W,. (26) 


The transformation 


Pel. thee call ct be 
i+-WwW, =~T+Ww, Ww, 


reduces the theory of disturbances to the utilization of 
the results obtained in SubsectionB [if (1/W,)Up is sub- 
stituted for Uj,], i.e., to a variant of Eq. (22): 

i ( i 


where analogs (23)-(25) enter £ Er after substituting Er, 
for U. 


Due to the possibility of zero setting between indi- 
vidual measurements and the limitation of the output 
instrument transmission band, the frequency spectra for 
the disturbances to be registered and for the useful sig- 
nals are comparable***, Hence, LEr<<En, [since 
© Er<<U im(22)], and the term L Er in (27) can subse- 
quently be neglected, 

By applying transformations similar to (17) and (18) 
to the left-hand side of Eq. (27), we can readily obtain 
a sufficiently accurate expression for Er, in developed 
form. In particular, for W, = Wg, it follows from (19), 
after substitution of Uj, for U,that 


t 
Ud 
r. suo pw, 


Er, =p Upt+ p+ F at 





The last two terms differ from (23) and (24) only by 
the fact that U}, was substituted for U, and, therefore, 
they are small simultaneously with (23) and (24) if only 
Up =U. For Te>1 sec, the first term is smaller thanthe 
other terms for PEF of the first group and is considerably 
larger for PEF of the second group. Thus, approximately, 


T 
7p Uy < Bry <t Ubmax (29a) 
for PEF of the first group, and 
T 
Er =~ # U, (29b) 


for PEF of the second group, The right-hand side of @9a) 
was obtained from (28) by replacing the integral by its 
maximum value for Up = Up max = Const and t= T,, under 
the assumption that the basic errors are of the same or- 
der of magnitude. 

Equation (27), including its variants (28), (29a), and 
(29b), should be considered as the initial equation for de- 
termining the parameters characterizing the quality of 
individual elements and of the entire instrument; the op- 
timum sensitivity, the maximum allowable line current 
and the amplifier drift, the required feed stability, etc. 

The amplitude distortions of the open circuit are 
formally equivalent to disturbances, and their effect can 
be taken into account by using the derived equations, 

The above concerns only circuits with feedback. In 
circuits without feedback, the disturbance level is not 
lowered, which constitutes their main efficiency. . 


Operational Stability of PEF 





The Mikhailov stability criteria [13] were found for 
the most general case of a circuit with combined feed- 





*** The operation of zero setting unavoidably cuts the 
output signal frequency spectrum from below without 
discrimination between the useful signal and the disturb- 
ance; the same holds for the output instrument (for in- 
stance, a self-recorder), which can be considered as a 
barrier filter in the high-frequency region. 



































back, The galvanometer operator was taken in the form 
given by (1). A variant of the Mikhailov criterion, ac- 
cording to which the conditions for the sequence of roots 
of the real and imaginary parts are found from the equa- 
tion 1+ We(iw)W,(iw)=0 (w is the angular frequency of 
signal harmonics), was used, The problem is solved 
exactly, since it is reduced to the solution of two quad- 
ratic equations, It is sufficient to use the approximate 
solution (see table). 

The basic stability criteria follow from the require- 
ment that the roots be real; 


1 (1 / =? 


and from their sequence conditions: 
T T To. 
41> Ta+ To (31) 


The other conditions (not shown completely in the 
table) cannot be disturbed, and they do not have to be 
considered separately, 

It follows from (30) and (31) that the measuring cir- 
cuit inductance, the galvanometer natural frequency, the 
feedback circuit constant ty, and the value of T(Ko, ac- 
cordingly restrict each other. It is also obvious that the 
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Ts and T, constants cannot be reduced below the deter- 
mined level (30)-(31). 

It should be noted that stability can be secured also 
if the galvanometer operates under periodic conditions 
and in the presence of positive feedback, 

In contrast to the method used in [7}, the system be- 
havior can be determined with respect to the equations 
derived for the general case, and, thus, the stability mar- 
gin can be determined beforehand [13], 


Discussion of the Results 

A. Creep and the Role of Rigid Feedback 

The determination and nature of creep follow from 
expressions (19), (22), and (23), Assuming that the struc- 
tural and differential errors are equal to zero (this is 
equivalent to Wy=const and Wq=Tgs), we find that creep 
characterizes those properties of the instrument which 
permit the simulation of PEF by an integrating circuit of 
the RC type, which is described by an equation of the 
type (14). The creep is caused by incomplete compensa- 
tion of the input signal by the feedback signal with re- 











TTT It was indicated in [7] that the stability of PEF does 
not depend on the restoring moment W,, of the galvanom- 
eter tension wires, This is due to the fact that T, as well 


as 1/48 are proportional to 1/Wrp. 

















spect to the derivative or its formal analog (in PEF with- 
out feedback). If additional rigid feedback is present, 
the input quantity compensation can be increased by 
using the feedback signal portion which is proportional to 
dU/dt. Thus, it follows from Eq. (9) that this occurs for 
A>0O and A~1/Kp. In this case, creep will be at the 
minimum, This compensation can also be reduced, 
which takes place for A<0, for instance, due to Rg, leak- 
age. In this case, according to (5), creep becomes larger. 

Thus, the reduction of leakage through Ry, and the 
introduction of a rigid positive feedback are effective 
means for reducing creep. 

The necessary condition for applying the other meth- 
od is a sufficient smallness of amplitude distortions in the 
open-loop systemt++. In the opposite case, due to the 
instability of Ky in (9), the 0< 1-AK g=const << 1 rela- 
tion cannot be maintained throughout the entire measure- 
ment cycle, 

B. Justification of Securing Integrating Conditions in 
the Amplification Circuit (in Circuits with Feedback) 

The reduction of creep in the open circuit (in the 
amplification circuit), i.e., the operation of this circuit 
under integrating conditions (as recommended in [10] 
and [11}), is not justifiable. 

In the first place, the expedients whereby T is aug- 
mented due to an increase in Ty are of little effective- 
ness when they are not accompanied by an increase in 
Ka’, since they are reduced only to an increase in KgTq, 
which is considerably smaller than Ty. This is clear from 
an analysis of Eqs, (8) and (11) if we take into account 
that the use of feedback is justified only if Cy < Cy. 

Considering the effect of integration conditions in 
the amplification circuit on the structural error (24), we 
can readily see by taking into account (18) that such a 
regime tends to increase the error determined by the rate 
of change in the magnetic flux to be measured, i.e., the 
allowable speed of operation of PEF is reduced for a cer- 
tain given accuracy (due to the increase in T,), At the 
same time, this regime sharply reduces the efficiency of 
feedback with respect to noise suppression [see (29a) and 
(29b) }. 

The (24) and (29a) errors are the smallest for w9(s)=1. 
In practice, this is realized (approximately) in PEF of the 
first group if we provide elements with the smallest time 
constants in the amplification circuit,or if we correct the 
amplification circuit with respect to an ideal amplifier 
with Wy=const, In this, the restrictions imposed by sta- 
bility conditions must be taken into account. 


C. Role of the PEF Theory in Design Calculations 
and in Developing Methods for Estimating the Measure- 
ment Accuracy 

The given PEF theory relates in the most direct way 
the basic instrument errors, which are expressed in terms 
of the output quantity, to its parameters, eliminating the 
necessity of double spectral resolution, as is the case 
when Fourier transformation or a step function is used, 




















This makes it easy, for instance, to calculate beforehand 
the allowable shape and duration of pulses for a certain 
assigned accuracy and parameters,or to solve any reverse 
problem. 


SUMMARY 

1, The simplest,and at the same time sufficiently 
perfect, PEF variant (also with the greatest operating 
speed) is the instrument with a negative series feedback 
with respect to the output derivative, and with the maxi- 
mum approximation (as far as the stability conditions al- 
low) to the ideal character of amplification in the open 
circuit. 

2. The PEF variant with combined feedback, includ- 
ing rigid positive feedback, secures a considerable reduc- 
tion of creep; however, it can be effectively used only in 
the case of sufficiently small nonlinear distortions in the 
amplification circuit, 

3. In PEF theory, along with frequency methods, it 
is convenient to use the method where the basic errors 
are expressed in terms of integrals and derivatives of the 
output quantity. This method provides the simplest solu- 
tion in analyzing the instrument operation, in design, 
and in error determination. 

The method can also be useful in analyzing the op- 
eration of integrators of different types. 

The author extends his thanks to S, P, Kapitsa, 

R. R. Kharchenko, S. G. Rabinovich, and M, I. Levin for 
their help in clarifying a number of principles stated in 
the present article, 
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++ This requirement was not fulfilled in any of the pro- 
posed circuits with combined feedback [3]. It is under- 
standable that the expected lowering of creep did not 
take place, 

‘This situation occurs, for instance, when Ty increases 
due to an increase in C; in the operation model in [10). 
In [11], T increases proportionally with T, due to the 
simultaneous increase in the tube amplifier amplifica- 
tion factor, This can be brought about also without an 
increase in Ty. 
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This article describes methods for the synthesis of circuits whose operation can be described 
by means of a certain special class of time functions, The determination of this class of 
functions and methods for their investigation are given. 


Great interest was shown lately in problems connect- 
ed with the synthesis and analysis of circuits and devices 
whose operation can be completely described by means 
of a certain functional time relation between the input 
and output parameters, For the solution of such problem: 
in the case where the circuit under investigation can be 
reduced to its equivalent trigger circuit, various methods 
have been proposed (see, for instance, [1], [2], and [5}). 
Problems in the analysis and synthesis of a certain spe- 
cial class of circuits, not reducible to trigger circuits, 
will be considered in the present article. 

We shall consider the totality of sets 


(24, Za, « 


eee Ty t), (1) 


where x; (i=1, 2,... , ) can be equal to either 0 or 1, 
and t can assume the value of any integer from 0 to s~1, 
inclusively, It is obvious that the total number of differ- 
ent sets (1) for a certain fixed value of n is equal tos 2" 

Determination. A function determined by sets of 
the form (1), and assuming one of the two possible values 
0 or 1, is called a Boolean time function, 

The statement that there is a finite number of dif- 
ferent sets (1) is equivalent to the statement that any 
Boolean time function can be represented in the form of 
a finite table. 





Example 1. 
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Despite its “clarity of presentation,” such a table 
determination is unsuitable in the majority of practical 
cases because it is cumbersome. When written in ana- 
lytical form, the function under investigation is more 
suitable for practical work. 

Consider a certain Boolean time function 

¥ = P (21, Ta, -- +s Ty» ft) (Og tqQs— 1). (2) 

If t is now assigned a certain fixed value t=k where 

0<sk<s-1, function (2) assumes the form 
Vy = Py (Fa Te +--+ By). (3) 

Function (3) is now an ordinary logic algebra func- 
tion [2 and 6), If we lett assume all the allowable val- 
ues, we obtain the sequence of logic algebra functions 


Por Pr +++» Dey (4) 


Thus, for any Boolean time function, we can obtain 
a corresponding sequence of logic algebra functions. 

Example 2. For the Boolean time function from 
example 1, we have 


Po = 21%_\V 2%2 V 2472 V m= i, 
Pi = 717%2 V 2y%9= 7}. 


Let us now introduce the function T,,, which is de- 
termined by the relation 


0, ta 
w={ a (0<a<qs—‘). (5) 


For any s-1, the following equalities hold; 


s—1 


V t=—1, V t,t, =0, 
Am tay 1S (6) 


which signify that at the instant of time t=a, the T , func- 
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tion, and only Tq, is equal to 1, while all other r ;(i # «) 
functions are equal to 0 at this moment. 

By taking into account relation (5), any Boolean 
time function can be written in the following manner: 


@P (2, Tar +++s Bye t) = PoTo V PT VeeeV Pe—1To—1° 
(7) 
Formula (7) provides means for writing any Boolean 
time function in convenient analytical form. 
Example 3, By writing the function from example 1 
in analytical form, we obtain 


P (715 22, t)=T. V 24T1- 


Determination, If in relation (7) all @(i=1, 2, 
. . « , SL) functions are represented in disjunctive (con- 
junctive,accordingly) perfect normal form, we shall 
say that the given Boolean time function is represented 
in a disjunctive (conjunctive) perfect normal form. 

From the theory of logic algebra functions [6], it is 
well known that any logic algebra function can be writ- 
ten in one of the two perfect normal forms. This justi- 
fies us in stating that any Boolean time function can be 
written in one of the possible perfect normal forms. 





Relation (7) and the above-formulated determina- 
tion indicate that the nature of Boolean time functions 
is close to that of ordinary logic algebra functions. In 
solving problems in the analysis and synthesis of circuits 
whose operation is described by such functions, this 
makes it possible to use an apparatus similar to that al- 
ready available for logic algebra functions. 

One of the most important problems in the logic 
algebra method of circuit synthesis is the problem of 
minimizing the given function, In the case of Boolean 
time functions, this problem can be solved by using any 
of the methods described in [2], [4], and [6] after certain 
preliminary modifications, For an illustration of this 
statement, we shall consider the example of minimiza- 
tion of the Boolean time function from example 1 by 
means of the modified method of minimizing charts. 

Example 4, The disjunctive perfect normal form for 
our function is 


@ = (Zat2 \V Zire \V xy%2 VV 2472) To V (F172 VV 7172) Tr. 


Let us apply the method of minimizing charts to 
this function: 
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According to this minimizing chart, the minimum 
form for our function is given by @ = 2z:\/to. 

As can be seen from this example, the proposed 
minimization method makes it possible to eliminate in 
the given function the terms which de not depend on 
time, which is useful in the synthesis of time circuits, 
However, for large values of n or s , the corresponding 
minimizing charts become too cumbersome and imprac- 
cal, In these cases, the method of approximate minimi- 
zation can be recommended, According to this method, 





Pig = V%- (8) 


It can be readily seen that any Boolean time func- 
tion 


Y= @P(Z1, T+ ++ Ty t) (ON tQae—1) 


can be additionally determined until a periodic function 
with the period q>s is obtained, On the basis of this 
and relation (7), we can state that, in synthesizing any 
circuit whose operation is described by a certain Boolean 


minimization is performed separately for all 4% functions _ time function, we obtain a circuit of the type shown in nec 

in expression (7). This method {s especially efficient for Fig. 1, In this circuit, switch S successively connects cit- pre 

small n and large s values. cuits which process the %,... , %-; functions, Afters mi 

Example 5, For the function of example 1, the ap- units of time, the cycle will be repeated, For the time ste 

proximate minimization yields the following expression: unit in the operation of switch S, we can take any quan- pli 
tity (for instance, the physical time in which one cycle 

P = To V 247s. in the operation of a certain device is performed). “ 


Determination. A Boolean time function is said to 
be periodic with the period q if, for any t, the following 
condition holds; 





In such a diagramatic realization of Boolean time 
functions, the elimination of terms that do not contain 
Ty, in the function to be synthesized corresponds to the 
case where the circuits processing these terms are con- 
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nected to the circuit output,bypassing the switch S. The 
presence of terms consisting only of Tg in the function 
minimum form corresponds to the case where, at the in- 
stant of time t=a, a voltage comparable to unity is sup- 
plied to the circuit output. 

In conclusion, let us consider the synthesis of a cir- 
cuit that processes an assigned Boolean time function. 


Example 6. 
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The disjunctive perfect normal form for this func- 
tion is 


= (x 2028 V% oe V 142972) 
To \V(airars \V Titers \/-e 172% \/ Zyx9r9 \V Trae) Th. 


As a result of applying the method of minimizing 
charts, we obtain 


P = 222t3 \V 2yFaTo VV (F1 V Zor) Ti. 


The corresponding functional diagram is given in 
Fig. 2. Its realization by means of electronic tubes 
yields the functional diagram shown in Fig.3, The nota- 
tion in this diagram is the same as in [4]. 
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The paper deals with systematization of various principles of designing multiplying devices 
with regard to increasing their accuracy and speed of operation. 

It contains some results of development of new multiplying devices carried out in the 
Laboratory of Analog Computers in the Institute of Automation and Remote Control of tne 


Academy of Sciences of the USSR. 


Introduction 

In solving many scientific and technical problems, 
it is now more often necessary to use computers which 
perform nonlinear mathematical operations, in particu- 
lar, multiplication and division, 

The design of an analog multiplying (or dividing) 
device which would satisfy modern technological stand- 
ards is a very complicated problem. This is indicated 
by the great diversity of the proposed principles on which 
the multiplication and division operations are to be based. 
However, regardless of the great attention paid to this 
problem, even the most perfect designs of multiplying 
devices are inferior to linear computer components with 
respect to the basic technical specifications, accuracy 
and speed, 

The contradiction between the accuracy and the 
speed requirements, which is especially pronounced in 
nonlinear computer components, makes it possible to 
classify the multiplying devices according to the scheme 
given in Fig. 1. Although such a classification is some- 
what arbitrary, since it is impossible to make sharp dis- 
tinctions with respect to errors and the frequency pass- 
band for each group, it is very convenient in selecting 
the most suitable multiplying device design for a certain 
problem, in developing new designs, and, especially, in 
combining new and old designs in systems where high ac- 
curacy and good dynamic properties are required. 

The classification (Fig. 1) is based on principles 
which make it possible to determine the limits of techni- 
cal characteristics of the devices in each group. These 
limits are ultimately determined by the principles used 
in different designs, and, therefore, the diagram in Fig,1 
in many ways reflects the chart composed with respect to 
the characteristics determining the method of obtaining 
the product [1]. 

The aim of the present article is the systematization 
of the various principles used in the design of multiply- 
ing devices with respect to possibilities of their further 
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development in order to improve their accuracy and op- 
erating speed,as well as the consideration of methods for 
improving these characteristics, Along with this, this 
article deals with certain results obtained in designing 
new multiplying devices at the IAT AN SSSR Laboratory 
for Mathematical Analog Computers during the past few 
years, 


l, High-Accuracy Multipliers with a Nar- 





row Transmission Band with Respect to 





One Input and a Broad Transmission Band 





with Respect to the Other Input 





With respect to dynamic properties, circuits where 
the transfer constant changes stepwise are the most prom- 
ising circuits for multipliers of this group. A circuit for 
a multiplier of this type is shown in Fig. 2[2). The up- 
per part of this diagram consists of a closed-loop follow- 
up system where such admittance yj, is tracked that the 
voltage Us at the output of amplifier 1, equal to 





v,-| U.yu + fo — Uyyis |, (1) 


where yy, Yig» and yyg are the admittances of the ampli- 
fier input circuit, y, is the admittance of the circuit 
parallel to the amplifier, Ux is the input signal, and U, 
and U, are constant stabilized voltages, does not exceed 
the absolute value of the actuation voltage Us» of the 
command supply device 5. For Us gy, 0, under steady- 
state conditions, the admittance y,, is equal to 


U yy 
U 


Uys 
+ U, , (2) 





Yu. = 
k 


In order to obtain the product, an additional admit- 
tance yg3, equal to y4»,is formed in the circuit (Fig. 2), 
and this admittance is commutated by means of an ad- 
ditional system of relay contacts or electronic switches, 
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which are centrolled by means of counter 4, If the bi- 
nary counter 4 has a sufficiently large number of places, 
and if y;= yy3 = Y2, the voltage Uz at the output of am- 


plifier 2 in Fig. 2 will be approximately equal to ey, 
k 


The deviation of voltage Uz from the voltage equal 
to the ideal product value under steady-state conditions, 
if all other conditions are ideal, obviously depends on 
the number of places in counter 4; 


bys 20 y 
Vv ¥s 2° —1 


AU, =U 


where AUz Is the absolute value of the error, Ayj, is the 
admittance increment per single counter unit, and n is 
the number of binary places, 

If we take the ratio of the error absolute value to 
the maximum output voltage for Uy = Uy max for the ac- 
curacy criterion, we obtain 


2U, max . 


6U = : 
U, max (2”—1) 





or, for Ux max = Uy max = Uk= Uz max’ 


2 
OU, max I= 4" (3) 





Fig. 2. 
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For a qualitative estimate of the dynamic properties 
of the circuit under consideration, we shall calculate the 
error in multiplying the sinusoidal voltage Uy = Ux max * 
X sinwt by the constant voltage Uy =Uymax 48 4 func- 
tion of the frequency f=u/2m of signal Ux for n places of 
counter 4, If we consider that the command supply de- 
vice 5 is practically inertialess, we find the maximum 
rate of change in admittance yy = yqy: 


dys hat 
— = foAtes a —Ooe 
( a aan igAvs 2” --4' 


where f, is the pulse repetition frequency of generator 6, 
The rate of change in voltage Uz for Uy =const is 


dt adt| “ys v Ys at 
therefore, the maximum value is 
(20s) _2U=man's 
dt /max a*—4 (4) 
If the rate of change in Ux» 
aU. 
ar = Ux max ® So,” (5) 


dU, 


does not exceed the found value max, the multiply- 


ing device will track the Ux signal with an error not ex- 
ceeding the static error, which is determined according 
to Eq. (3). 

For signal U, frequencies higher than the limiting 
= a =) ' the error sharply increases, In 
this case, the output voltage of amplifier (2) (Fig. 2) will 
change according to a curve which coincides with a sinu- 
soid with an accuracy corresponding to the unit of the 
first place of the binary counter only over certain por- 
tions, while changing linearly over the remaining por- 
tions (Fig. 3). 


frequency f,= 
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The sinusoid and the straight line begin to diverge 








27k 2m 
at the instances of time t,+ 7 and ty+ = , when the 


sinusoid slope is equal to the maximum rate of change in 
the output voltage. 
From Eqs, (4) and (5), we find 





2f 
@t,; = arc cos [- — 6 | 


2"—1)ol° (6) 


If we know t,, the equation of the straight line issu- 
ing from the point Uz, can be written in the form 


(t¢ — ty) (7) 


aU, 
U,=U,,- (a). 


or 





U,=U Usmax| V tox ee ay — 580m). 


At the moment when the slopes of the straight line 
and the sinusoid are equal, for tr, =@1/w)-t;, the error 
is at the maximum and is equal to 


bu = eR are 008 = 


zmMax, (2" — 1) af (2" — 4) xf, 
(8) 
—2)V 1-4, ey Sot 


Figure 6 (curve a) represents the diagram of error 
variations 5 Uz (f) for an eleven-place binary counter, 
which is supplied with pulses at the repetition frequency 
fg=1 kc. Up to the frequency f=f,* 0.16 cps, the error 
is somewhat lower than 0,1%: however, for a frequency 
of 3 cps, it already attains 60%, 

Thus, the circuit under consideration, while secur- 
ing a rather high static accuracy, has very low dynamic 
properties (with respect to the Ux signal), and the possi- 
bilities of improving these properties are very limited, 
since the binary counter filling frequency and the upper 
operating frequency of the switching elements which 
control the yy, and y,, admittances (Fig. 2) must exceed 
the multiplier transmission band by more than two orders 
of magnitude, 

The dynamic properties of the circuit can be im- 
proved by using the method of place-by-place filling of 
the binary counter, which is sometimes used in convert- 
ing a continuous quantity into a discrete quantity, A 
variant of the structural diagram for such a device is 
shown in Fig. 4. In comparison with the preceding cir- 
cuit, this circuit contains an additional binary counter 
8, and a command switch 4, for which a binary counter 
with scaling to a number of units equal to the number of 
places of counter 7 can also be used, 





In the initial position, the command switch 4 does 
not transmit a pulse from generator 6 to the first place 
of counter7 as usual, but to its upper nth place, where- 
by it is transferred from position 0 to position 1, After 
this, the error voltage Us at the output of amplifier 1 is 
measured by the command supply device, and, in depend- 
ence on the magnitude and sign of the error, either the 
command for switching the nth place to position 0 or 
the command for its remaining in position 1 is given, At 
the same time, 1 is added in the counter of the com- 
mand switch 4, The next cyclic pulse from generator 6 
transfers the (n~1)th place of counter 7 into position 1; 
the Us voltage is then measured again, etc., until all of 
the n places are filled, 

Before the beginning of a new cycle, the contents of 
counter 7 are first transferred by a cyclic pulse to count- 
er 8 and are then tripped out to zero, Thus, the contents 
of counter 8, which controls admittance y,;, change only 
after the (n+ 1)th cyclic pulse from generator 6, but 
they change immediately by the number of units corre- 
sponding to the rate of change in signal Uy. In connec- 
tion with this,also,the character of the dynamic error 
changes, For the case where the Uy =U x max sin wt volt- 
age is multiplied by Uy =Uy max, the increase in error 
in comparison with the static error begins at a frequency 
which is n+1 times less than that used in the circuit in 
Fig. 2, i.e., at a frequency of only 0,15 cps; however,it 
rises much more slowly, 

If, in the first approximation, we assume that the 
nodal points at the end of each operating cycle are ac- 
curately tracked, the maximum error will change in di- 
rect proportion to the rate of change in Uy, i.e., 


_ anf (n + 1) 
OU, max * - reas os (9) 


Up 
*U, 


-U, 
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In this case, a 10% error will occur only at the 
1.35 cps frequency, while in the preceding case, such an 
error occurred at the 0,18 cps frequency (see curve b in 
Fig. 6). 

The deficiencies of the described method are: a con- 
siderable amount of additional equipment and more 
stringent requirements imposed on the transmission band 
of amplifier 1, which must transmit a rectangular signal 
with the frequency of the cyclic pulses from generator 6. 

A more perfect method for increasing the operating 
speed of circuits where the transfer constant changes step- 
wise is the method of automatically changing the num- 
ber of places in the counter, The essence of this method 
becomes clear if we consider the circuit in Fig. 2, Actu- 
ally, if, for n places in the counter in this circuit, a con- 
siderable increase in the dynamic error begins with the 
signal U, frequency equal to f,, then,for the (n-1)th 
place, this increase starts at the 2f, frequency, while the 
error before attaining the 2f, frequency also increases 
twofold. In the successive disconnection of lower count- 
er places in the circuit, for Ux change rates equal to the 
maximum rate of change in Uz (4), the maximum dy- 
namic error in following the sinusoidal voltage will vary 
in jumps at the limiting frequencies 





h= aot (end 


In the frequency range from fj to fj4;, it willremain 
constant and equal to 


(10) 


2 
8U ,max = grt _ 4° 


The rate of change in voltage U, in the circuit 
shown in Fig, 2 can be perhaps most suitably estimated 



































with respect to the magnitude of the error voltage Us, 
which is demonstrated by the structural diagram shown 

in Fig. 5. For the disconnection of lower places in 
counter 4, the circuit is provided with anadditional com- 
mand supply device. 7, and a reversive counter 8, with 
a number of places securing the storage of a number of 
units equal to the number of places to be disconnected, 

The command supply device 7 starts to fill counter 
8 when the error voltage Us is somewhat larger than 
|2U56| ,and to subtract units from it when |Us|< |Usl- 
For each addition of a unit to counter 8, the admittance 
y; increases twofold,and the supply of commands is 
switched to a higher place of counter 4, 

The automatic variation of the admittance y, se- 
cures the correspondence between the actuation voltages 
of the command supply devices 5 and 7 (which remain 
constant during operation) and the number of places in 
counter 4 (Fig. 5). 

The described method ensures the multiplication of 
the sinusoidal signal U, by a constant quantity at fg=1ke 
and n=11,with an error less than 6.4% at frequencies of 
up to 10 cps,and with an error of less than 12.8% at fre- 
quencies of up to 20 cps, etc, (Fig. 6c). Moreover, the 
application of this method does not involve more strin- 
gent requirements for the quality of amplifier 1, and the 
use of two multiplace counters is not necessary, 

Finally, another way of improving dynamic prop- 
erties of the circuit in Fig. 2 is the use of the time- pulse 
interpolation effect [2], 

When the circuit in Fig. 2 operates under self- oscil- 
lation conditions, and the number contained in the bi- 
nary counter alternately increases or decreases by unity, 
the yy2= yg; admittances change linearly ,and not in steps, 
as functions of Ux; therefore, the over-all number of 
counter places can be reduced without diminishing the 
accuracy, However, in this case, pulsations with the fre- 
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quency of self-oscillations whose amplitude sharply in- 
creases with a decrease in the number of places arise at 
the device output; therefore, the possibilities of broaden- 
ing the transmission band by this method are limited. 
The frequency-response curve of the circuit in Fig. 2 
under self-oscillation conditions for n=8 and f= 1 ke is 
shown in Fig. 6d. 


2. Accurate Narrow-Band Multipliers 











The most effective method for eliminating the in- 
fluence of intrinsic tube parameters in linear decision 
elements is the application of circuits with deep negative 
feedback. 

Multipliers based on the negative feedback princi- 
ple also secure a comparatively high accuracy; however, 
at the present time, the application of this principle is 
possible only in circuits where the output signals are first 
transformed into the amplitude or duty ratio of rectangu- 
lar oscillations [14 and 13], in circuits where different 
methods of signal modulation by sinusoidal voltage are 
used [14 and 22), and also in certain circuits with spe- 
cial electron-beam tubes [56 and 61}, 

As a rule, in replacing input voltages by signals of 
different shapes and higher frequencies, the production 
of final results requires also reverse transformation,which 
necessitates the use of detectors, frequency or phase dis- 
criminators, circuits for producing the average voltage, 
etc; therefore, the transmission band of multipliers with 
preliminary transformation of signals isessentially limited 
to the region of audio-frequencies and, in many cases, to 
regions of considerably lower frequencies, 

Up to now, the best technical characteristics of time- 
pulse devices have been secured by using circuits with 
stabilized electronic keys [7], for which, according to 
the data provided by the authors, the accuracy was 0, 1% 
for a transmission band of up to 200 cps. However, such 
circuits are extremely complicated, and they require a 








large number of decision amplifiers. In connection with 
this, an investigation with the aim of finding ways for 
improving the quality of time-pulse devices by using 
simpler means has been conducted at the Institute of 
Automation and Remote Control, Academy of Sciences, 
USSR. Figure 7a shows the block diagram of a multiplier 
with combined connection of time-pulse dividers, The 
structure of this circuit is similar to the circuit proposed 
by A. A. Fel'dbaum and A, L Manukhin [6]; however, 
the electronic keys are connected simultaneously to the 
input circuits and to the circuits parallel to amplifiers 1 
and 2, The upper part of the circuit in Fig. 7a consists 
of a closed-loop follow-up system where the factor Uy 

is transformed into relative operation times of keys Ky 
and Ky. Trigger Tr is switched from one steady state to 
another with a frequency equal to the frequency fsg of 
sawtooth oscillations. The relation between times t, and 
t, (Fig. 7b) for each of the trigger states depends linearly 
on voltage U,.due to changes in the duration of the total 
U1 + Usg voltage, which is smaller than the trigger switch- 
ing voltage Us,,and larger than Usp. 

The keys are controlled by the trigger in phase op- 
position in such a manner that, when key K, is open, Ky 
is closed, and conversely, The average value of the Uys 
voltage at the key Kg output, to which voltage U, is sup- 
plied, is proportional to its relative open-state time 
a=t,/T (Fig. 7b), f.e., Uke ay = ~Up%. Correspondingly, 
Uk av = Uy (1-a). 

For a sufficiently large amplification factor of am- 
plifier 1, such a voltage U, is obtained at its output that 


U, (i—a) U.o 
Ri+R,, ~~ Rit Rp 


a Ry + 
ia "sit K)* 


or 











































































































For the lower portion of the circuit, the following 
relation holds: 


a Roath 
Ys Toa ER 


By substituting here the +s value from (11),we 
have 


Instead of obtaining the duty ratio € of rectangular 
pulses, which is usually determined by the equation 


t _ 
Es as , the voltage U, is here transformed into the 


value of the fraction ——— 2 . If the range in which 


the duty ratio € changes is limited by definition by a 
number less than 1, the fraction — can change from 


zero to infinity. In practice, the duty ratio maximum 
value € may is limited by a number less than 1, and the 


value of the fraction ee is limited by a finite quan- 


tity, since the time (Fig. 7b) cannot change smoothly if 
it becomes smaller than the sum of the leading and 
trailing trigger front times ty+t,. If ty and t are given, 
€ max Will linearly decrease with an increase in the saw- 
tooth oscillation frequency fsg. since 


t her 
g = -tnex “amin 1 —2tmingg (12) 


and the upper limit of changes in 7 is related to 
the fg frequeticy by the equation 


(rw ® — Fe (13) 





By comparing (12) and (13), we can readily see that 
the value of (5) 
1-& 


max 

greater than € max for a large fsg frequency. For equal 
product scales, this makes it possible to reduce the trans- 
fer constant of the output decision amplifier, by means of 
which the error caused by the keys can be considerably 
reduced. If the fsg frequency is increased, filters for a 
broader transmission band can be designed, The simul- 
taneous connection of filters in input circuits and circuits 
parallel to amplifiers 1 and 2 (Fig. 7) also contributes to 
the improvement of the circuit dynamic properties due 
to the forcing action of filters Rf, Cfp, Re, and Rg, Cfs, 
Rg. 

An experimental check of the mock-up of a multi- 
plier circuit with combined connection of pulse dividers, 
which was based on EMU-8 decision amplifiers, indi- 
cated the possibility of securing the following character- 
istics after careful circuit adjustment; 

1) a maximum static error of 0.2 v if the factor 
voltages Ux and Uy and the output voltages Uz change 
in the interval from 0 to 100 v; 

2) drift of less than 40 mv over a period of 10 min 
when U, =100 v is multiplied by U, =0; 

3) a drift of not more than 70 mv over a period of 
10 min when U,=100 v is multiplied by Uy =0; 

4) a drop in the amplitude response in both chan- 
nels of not less than 5% for a frequency of up to 400 cps 
and for a commutation frequency of 10 kc, 

The wiring diagram of the mock-up is not given 
here, since it was constructed only for checking the pos- 
sibility of realizing the advantages of the considered 
structural scheme. 

As another example, we shall consider the logarith- 
mic time-pulse multiplier. The block diagram of the 
device for the case of single-quadrant multiplication is 
shown in Fig. 8a. Here, the transformation of the U,, 
Uy, and Uz voltages into relative durations of rectangu- 
lar pulses at the output of triggers Trl and Tr2 is secured 
by comparing them with exponentially changing voltages, 









































































































































canbe chosen to be considerably greater 
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The pulse generator PG switches triggers Trl and Tr2 and 
actuates the exponential voltage generator EG1 whose 
output voltage ue is compared with the U, voltage 
by the comparison unit CU1. At the moment when these 
two voltages are equal, which occurs after the time t,= 


=-rinZ& , & pulse which triggers EG2 appears at the 
0 


CUl1 output. A pulse at the CU2 output, which returns 
the Trl trigger into the initial position, appears after the 


time t,=-Tln be & . 
0 


The duration of the positive pulse at the Trl output 
U,U 
is equal to the sum of the times t,+t,=-TIn = . 
0 


A pulse of negative polarity with an amplitude 
equal to the amplitude of the pulse supplied by trigger 
Trl is tapped off the output of trigger Tr2; therefore, the 
Uz voltage at the amplifier output is such that its dura- 
tion is ty=t,+ t. 

Uz _ UU 
qe a. 

If multiplication in four quadrants has to be secured, 
all circuit units have to be doubled and adapted for work 
with negative voltages, 

The ease with which one can switch over to division 
operations and the considerable simplification in the case 
of the multiplication of several voltages constitute the 
advantages of this circuit. 

With regard to accuracy, this circuit occupies a mid- 
die position between closed-loop and open-loop devices. 
With respect to open-loop circuits, the error of which 
usually exceeds 1%, the gain in accuracy is obtained by 
replacing the sawtooth voltage by an exponential voltage, 
which can be secured with very high accuracy. 


Since tg=-TIn 































































































Uy my 
U, sinwt 2; - Uy 
U; mM, 
a — 
Uy cosest Up 
| 2; PhD y > 
m, |G “e 
Up : + 











As in other time-pulse devices, the transmission 
band remains lower than the frequency of PG cyclic 
pulses by more than one order of magnitude. 

An experimental check of the multiplier mock-up, 
which was constructed in accordance with the structural 
diagram and where standard elements used in computer 
units were utilized,yielded the following characteristics: 

1) The accuracy was 0.5% with respect to a maxi- 
mum output voltage of 100 v; 

2) the frequency response remained horizontal for 
from zero to 20 cps for a commutation frequency of 1 kc. 

As a rule, circuits based on different methods of sig- 
nal modulation [14-22] are considerably more compli- 
cated than time- pulse circuits; however, due to the possi- 
bility of selecting a frequency of sinusoidal carrier oscil- 
lations which is considerably higher than the frequency of 
rectangular oscillations (from 100 to 1000 kc), these cir- 
cuits have a broad transmission band (up to 1-10 kc), 
Attempts further to increase the carrier oscillation fre- 
quency lead to a sharp deterioration of accuracy, since 
the mutual influence of individual units becomes more 
pronounced, The design of, and comparison between, 
different variants of such circuits are given in [17] and 
[18]. 

Figure 9 shows a simplified variant of a multiplier 
circuit with amplitude modulation of signals; the circuit 
operating principle is based on phase relationships. Two 
high-frequency voltages with a phase shift of 90 deg are 
modulated by voltages Ux, Uy, and Uz by means of mod- 
ulators m,, m,, and mz. oe 

The phase of the vector sum of voltages Ug=U,sinwt 
and b, =m, U,coswt is compared by the phase detector 
PhD with the phase of the vector sum Dy+ Uz=myUy x 
X sinwt + mz,U,coswt (Fig. 9b). If the phases of these 
two sums do not coincide, an error voltage, whichcauses 
the output voltage to change until it assumes a magni- 
tude proportional to the product of Ux and Uy, arises at 
the amplifier input. 

With respect to accuracy and operating speed, such 
a circuit does not have any advantages over other devices 
with signal modulation, but it is considerably simpler 
and less expensive than the latter,and it secures a direct 
multiplication of factors with alternating signs, The fact 
that no stringent requirements are imposed on the ampli- 
fier transmission band, since the amplifier amplifies the 
already detected signal, also constitutes an advantage of 
this circuit. 

The deficiencies of this circuit are the necessity for 
an accurate adjustment of circuits in the m ,, my, and 
m,, modulators and the requirement fer a high carrier fre- 
quency stability. 

An experimental mock-up of the multiplier con- 
structed according to the diagram in Fig. 2 secured an 
accuracy of 0.5% for a transmission band of up to 5 kc. 
The carrier frequency was 465 kc. 

By using a quartz oscillator, the device accuracy 
can be additionally improved, 
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3. Medium-Accuracy Multipliers with a 














Medium Transmission Band 

This group comprises indirect- action multipliers, 
where the multiplication operation is replaced by sum- 
mation and functional conversion operations [23-50], 
and multipliers based on the Hall effect in semiconduc- 
tors (52 and 55}. 

The possibility of improving the accuracy of indi- 
rect- action circuits is limited by difficulties involved in 
constructing accurate nonlinear converters, Up to now, 
the best results (0,5- 1%) were obtained with circuits 
based on the principle of sectionally linear approxima- 
tion. The transmission band of indirect- action devices 
is limited by the parasitic parameters of functional con- 
verter elements. By using special corrective circuits, it 
can be brought up to a few tens and, in individual cases, 
even to a few. hundreds of kilocycles [27, 37, and 38}. 

In circuits where the Hall effect in semiconductors 
is used, an accuracy of 0.3-1% can be secured for a 
transmission band from 1 to 10 kc [52-55]. 

Multiplier circuits with quadratic functional con- 
verters are most widely used at the present time, which 
is explained by their relative simplicity and the possi- 
bility of application for solving a relatively large num- 
ber of problems. 

The following circuit types are used for squaring: 
diode circuits [25, 28, 30, 31, 35, 36, 37, and 39), non- 
linear semiconductor resistors of the thyrite or willemite 
type [32, 34, 46, 47, 48, 50, and 51], circuits with saw- 
tooth voltage detection [23, 24, 33, and 42], electron- 
beam functional converters, etc. The quadratic sections 
of the volt- ampere characteristics of semiconductor or 
vacuum diodes, the quadratic dependence of the anode 
current on the grid voltage in three-electrode tubes [25, 
29, 43, and 44}, etc, can also be used. 

Different variants of multiplier circuits based on 
quadrators as well as surveys of these circuits with indi- 
cations of past and possible future improvements, have 











been sufficiently treated in literature (see, for instance, 
{1]}), and, therefore, we shall not consider them here. 

In logarithmic multipliers (15, 63, and 65], the prod- 
uct is obtained as a result of finding the antilogarithm of 
the sum of logarithms of the input quantities, Such de- 
vices are usually constructed according to the diagram 
shown in Fig. 10. The determination of antilogarithms 
is here secured by connecting the logarithmic functional 
converter +In to the circuit parallel to the decision am- 
plifier. Voltages U, and Up are added to the U, and Uy 
factors in order to secure the multiplication of quantities 
with alternating signs, since the logarithms are deter- 
mined only for positive argument values, In this, the U; 
voltage at the output of amplifier 3 is set in such a man- 
ner that 


— WU, + U,)—W(U, + O,) = In Us, (14) 


i. 


Us =(U, + U,)(U,+U,). 


For Ua> |Uxmaxl and Up> |Uymax| . all loga- 
rithms will be determined; however, in order to obtain 
the final result, the UxUp, UyUa, and UgUp components 
must be subtracted at the circuit output. 

The character of error variations in the circuit is of 
theoretical interest. If we assume that all logarithmic 
converters have errors whose absolute values are Aj, (Uy), 
Ain(Uy), and A;,(Us), and the other operations are per- 
formed without errors, Eq. (14) should be rewritten in the 
following form: 


—In(U, + UU, + U,) = 


= In Us + Inexp[Aj, (UV) + Aig (Uy) + Ain (7 ))- 
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Fig. 10. 








Hence, 


Ug= 


=(7, + UU, + U,) exp (Ay, (U5) + Ayn (Oy) + Stn (Us)I- 


The Ayn(U,), Ain(Uy), and Ajn(Us) errors are small; 
therefore, 


Us (Ug + UU, + ,) 


(1 + Ain (Og) + Ain (Oy) + Aiy (Us))- es 


After subtracting the UgUp, UaUy, and UpU, com- 
ponents, we obtain 


8(U,) = 
Uy —Uzg «UU 
al =U; 





= (Ain (Ug) + Aig (Uy) + Sin (Ua) 
(16) 


Uz max 


Here, Uz; is the ideal product value, and Uza is the 
actual product value. 

If we assure that the errors of functional converters 
are constant and equal and that they are independent of 
the argument values, the relative product error will be 


(U,U,+U,U,+ U U0) 
UG. max 


U,U, 
t= Uz max 4m + 





3Ain- (17) 


Thus, the error depends linearly on the magnitude of 
the product itself and on each of the factors, and it con- 
tains a constant component, Physically, this is explained 
by the change in the transfer constant of amplifier 3, 
which is a consequence of connecting the amplifier to 
the functional converter feedback circuit. 

In the case of single-quadrant multiplication, when 
Ux and Uy have only one sign, the U, and Up voltages 
can be assumed to be equal to zero. Then, the errorwill 


decrease with a decrease in the product value, and its 
maximum value will not exceed 344). Actually, for 
small argument values, this error will be larger, since the 
slope of the characteristics of logarithmic converters isin 
this case larger, and their error sharply increases, 

The considered properties of logarithmic multipliers 
can be used when it is necessary to maintain a higher ac- 
curacy in a certain given range of argument changes, 

From among the useful properties of logarithmic de- 
vices, one should mention also the possibility of direct 
multiplication of several factors, division, squaring, evo- 
lution, etc. However, their widespread use is hampered 
by difficulties in producing simple and sufficiently accu- 
rate logarithmic functional converters, 

The following property of semiconductors is used in 
circuits based on the Hall effect; if a semiconductor 
plate is placed in a magnetic field and a current is trans- 
mitted through it, an emf will arise in the plate, the di- 
rection of which will be perpendicular to the current and 
the magnetic flux, i.e., E=k BL. 

A number of multiplier circuits based on this princi- 
ple are already in existence [52 and 53}; however, best re- 
sults are obtained by using closed-loop circuits, such as, 
for instance, the circuit shown in Fig. 11. Here, two 
equal elements, I and II, are placed in a magnetic flux. 
A steady current I, is transmitted through the first ele- 
ment, The Hall ernf in the first element is compared 
with the voltage Uy by means of a differential amplifier, 
whose output voltage causes the induction in the mag- 
netic circuit gap to change in such a manner that k,E,= 
=kek;BIg= U, , where k, is the Hall constant, and kg is the 


proportionality factor. Hence, B= 8 


Kykaly © 
BB woe UxU 
If the “shot the U, voltage is equal to —,1- 


The circuit structure is the same as for other circuits 
with a controllable transfer constants however, the con- 
trolled element is the variable emf. 
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The circuit itself is purely electronic, but without 
the preliminary conversion of input voltages, This 
makes it possible to improve its dynamic properties in 
comparison with circuits where input voltage modulation 
is used; however, up to now, it was possible to secure 
transmission bands of only up to 10 kc, The reason for 
the limited possibility of broadening the transmission 
band in the Ux channel is the presence of the time con- 


stant T = =( is the coil inductance, and r is the coil 


resistance plus the differential amplifier output resist- 
ance) in the differential amplifier feedback circuit. 

If the differential amplifier gain is sufficiently large, 
its output voltage will be related to the Uy voltage by 
the forcing element equation 


The proportionality factor kg depends on the Hall 
constant k;, the current I, the magnetic system charac- 
teristics, the differential amplifier gain with respect to 
the E, signal, etc, The enumerated factors make it pos- 
sible to vary ks in a sufficiently wide range; therefore, 
we shall take its value to be equal to 1 in order that 
Ug = U,, under steady-state conditions. Then, if U, 
changes, the voltage Up will additionally increase in pro- 
portion to the rate of change in U,. 

If the differential amplifier linearity is limited by 


the value of Up max, the maximum rate (“28 max can 
be determined by 





Up UO, (Pp) 
PU, (P) ax x a (19) 


If we assume that, in the transient process, Up max 
can exceed k times the maximum value of Up under 


steady- state conditions Up op = Ux max(Ug max kUp op). 
then, for determining the limiting value of the frequen- 


cy fy = sil of the sinusoidal signal Ux = Ux max sin wt, 
the discriminant of the equation 





w* (1 — sin* wt) — (+ - ey = 0, 


solved with respect to sinwt, should be put equal to zero. 
After this, we obtain 


i 
4, = tar VE—1. (20) 


The maximum value of the fj}; frequency depends 
on the possibility of increasing k as well as the possibil- 
ity of reducing T. 

With an increase in k, the operating range of 
changes in U, is narrowed, and, consequently, the accu- 
racy deteriorates, 


1010 


The possibility of reducing the time constant T ts 
limited for reasons of energy, since, if the rate of change 
= increases, also the power taken from the differential 
amplifier increases. 

Work proportional to the square of induction A=aB? 
is used up in creating the magnetic flux. Hence, the 
instantaneous power is 


aA dB 
ap = 2B a = Val. (21) 


For given values of induction Bop and Up max: 


(=) — Vemax!s max (22) 
dt) ax 2aB 


Consequently, the possibility of increasing the cir- 
cuit operating speed also by increasing ly max or, which 
is the same, by decreasing T, is limited by the maximum 
amplifier power. For these reasons, multipliers based on 
the Hall effect have been grouped together with devices 
with a medium transmission band. 

The accuracy of such devices is limited by instabil- 
ity and a considerable dependence of semiconductor prop- 
erties on temperature, 

At the present time Ge, Si, InSb, and InAs are used 
as Hall elements. Germanium and silicon elements have 
the largest Hall constant, and indium antimonide and 
indium arsenide make it possible to obtain a more power- 
ful output. In comparison with other elements, the prop- 
erties of InAs change slightly with temperature. 


4. Approximate Broad-Band Multipliers 

Such devices are provided with electron-beam tubes 
of a special design [56-61 and 66] or parametric linear 
circuits, where a current proportional to one factor is 
transmitted through a resistor which is controlled by the 
other factor [15]. 








U; Uy U, 


Fig. 12, I) Beam generating system; I) deflecting 
system 1; III) deflecting system 2; IV) deflecting 
system 3; V) anode system; DA) differential amplifier. 
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In circuits of the second type, usually the intrinsic 
characteristics of electron tubes are used, and, therefore, 
their accuracy is very low (5-10%). 

A high*accuracy (up to 0.5- 1%) can be secured in 
circuits with electron-beam tubes, Here, it is extremely 
difficult to reduce the error, since it depends on the pre- 
cision in manufacturing and assembling the tube elec- 
trodes, 

As an example, Fig. 12c shows the circuit of a mul- 
tiplier whose design is based on an electron-beam tube 
with a hyperbolic deflecting field [58]. In contrast to 
the uniform electrostatic field, the action of a hyperbolic 
field on the electronic beam is proportional, not only to 
the field strength, but also to the beam displacement 
with respect to the tube axis, 

In the circuit in question, the electron beam is first 
deflected in the horizontal direction proportionally to 
the Ux voltage by means of an ordinary deflecting sys- 
tem, after which it enters the hyperbolic field, which is 
created by the factor Uy. The resulting beam deflection 
in the horizontal direction is proportional to the UxUy 
product, while the direction of beam deflection corre- 
sponds to the sign rule in multiplication. The hyperbolic 
field is created by four electrodes, whose sections with a 
vertical plane yields equal hyperbolic branches. 

In order to obtain results in voltage form, a third de- 
flecting system and an anode system are provided in the 
tube. When the beam is deflected with respect to the 
cut in the middle of the anode system, the equality of 
currents tapped off each of its halves is disturbed,and a 
voltage Uz, which returns the beam into its initial posi- 
tion, arises at the differential amplifier output. As a re- 
sult, we obtain the following relation; Uz, =kUxUy. Ac- 
cording to the authors’ data, the accuracy of such a de- 
vice is 0.5%. The transmission band is 200 kc. 

In another circuit, which is shown in Fig. 13 [56], the 
following well-known principle is used; an electron mov- 
ing rectilinearily at the speed V in a magnetic field with 
the strength H is acted upon by a force equal to Hx ¥, 
Voltage Ux is applied to plates 1, The voltage Uy iscon- 
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Fig. 13, I) Beam generating system; II) deflecting system 1; 


Ill) deflecting system 2; Ph) photomultiplier. 





verted into current I,, which passes through the coil and 
creates the magnetic field, The greater the deflection of 
the beam by plates 1, the stronger the action of the mag- 
netic field. The resulting action on the beam is propor- 
tional to the product of U, and H. One-half of the screen 
in the tube is darkened, and, therefore, different signals 
are received by the photocell, in dependence on the beam 
position. As in the preceding case, there is a control cir- 
cuit, due to which the Lorentz field, which acts on the 
beam electrons, is compensated by the amplifier A out- 
put voltage, which is equal to Uz=kU,Uy. The error of 
such devices is equal to approximately 2%, 


In the multiplier proposed in [60], an electron-beam 
tube with an ordinary electrostatic deflecting system is 
used, A plate, cut into four, as shown in Fig. 14, is 
mounted in the forward part of the tube. A current pro- 
portional to the surface area of the corresponding beam 
portion is tapped off each of the four plates, The currents 
are then summed according to the signs indicated in the 
figure. If the current density of the entire beam is uni- 
form, and the factor voltages Uy and Uy are supplied to 
the deflecting plates, the total current will be propor- 
tional to their product, and the rule of signs will be ob- 
served, The fact that the beam cross section is not square 
has no importance, since the sum of segments comple- 
menting a square to form a circle is always equal tozera 
The nonuniformity of the beam current density affects 
the device accuracy, which is low; it is equal to 3%, 
According to the authors’ data, the transmission band is 
equal to 200 kc, 


There is still another electron-beam tube variant 
[59], whose operating principle is illustrated in Fig, 15, 
which is intended for use in circuits with quadrators, 
The tube has a pair of deflecting plates, which direct a 
rectangular beam to a profiled parabolic screen. The 
beam current tapped off such a screen is proportional to 
the square of the voltage supplied to the plates. The 
beam operates with voltages of either sign, and, there- 
fore, there is no necessity for separating the moduli in 
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Comparative Table of Multiplier Technical Characteristics 








S tatic Transmission 
Operating principle of the device accuracy, % | band, cps 
Very accurate, narrow-band with respect to one 
input,broad-band with respect tothe other input 
Circuits with potentiometers which are controlled by a follow- up system 0.1 0.5 
Circuits where the transfer constant changes stepwise 0.1 0.2 
The same circuit with the application of the time-pulse interpolation 
method 0.1 1,2 
The same circuit with place-by-place filling of the counter - 01 1,3 
The same circuit with automatic changing of the number of counter 
places 0.1 10 
Accurate narrow-band multipliers 
Time-pulse multipliers [6] 0.4 12 
Circuits with stabilized electron keys [7] 0.1 200 
Logarithmic time-pulse devices 0.5 50 
Tithe-pulse devices with combined connection of pulse dividers 0.2 400 
Multipliers with amplitude modulation of input voltages [16] 0.1 10000 
Multipliers based on phase relationships 1 5000 
Devices of medium accuracy with a medium 
transmission— band 
Multipliers based on diode quadrators 0.5 10000 
Multipliers based on thyrites 1 3000 
Multipliers based on the Hall effect ‘ 0.5 1000 
Approximate broad-band devices 
Multipliers based on Ohm's law 10 2-105 
Devices with electron-beam tubes 0.5 10° 
Combined devices with parallel channels 
Circuits with amplitude-separation of channels 0,02 100 
Circuits with frequency-separation of channels 0,02 10000 
and over 
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Fig. 15. I) Electron gun; II) deflecting plates; Ill) profiled screen; 
IV) collector. 
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using such a tube, According to the authors’ data, the 
accuracy in squaring is 1%, and the transmission band 
attains several megacycles per second, 


CONCLUSION 

1, The above classification and survey of multi- 
pliers indicates that the use of any single principle in the 
design limits the maximum capabilities of these devices, 
either with respect to accuracy or with respect to operat- 
ing speed (see Table of Technical Characteristics), 

2. The considered well-known methods of increas- 
ing the operating speed of circuits where the transfer con- 
stant changes stepwise make it possible to broaden their 
transmission band only to a slight degree (from 0,2 to 1.3 
cps), The method of automatically changing the number 
of binary counter places, which has been proposed here, 
makes it possible to improve the dynamic properties of 
a digital follow-up system up to 10 cps for a constant 
value of the product of accuracy and the transmission 
band, 

3. In analyzing the possibilities of improving the 
accuracy and the operating speed of narrow-band multi- 
pliers, several new circuits have been proposed; 

a) a circuit with combined connection of time-pulse 
dividers, which, for an accuracy of 0,2%, secures the 
broadening of the transmission band of time-pulse cir- 
cuits to 400 cps; 

b) a logarithmic time-pulse circuit, which secures 
an easy changeover to division operations and the possi- 
bility of direct multiplication of several factors; 

c) a multiplier with signal amplitude modulation, 
which is based on phase relationships and which is char- 
acterized by the simplicity of its circuit and the possibil- 
ity of direct four-quadrant multiplication, 

On the basis of the obtained results, we can reach 
the conclusion that the design of a broad-band and, at the 
same time, accurate multiplier should not be based on 
any single principle, but on a combination of different 
circuits into a system which would make it possible to 
use the good dynamic properties of one principle and the 
high accuracy of another. 

Combined multipliers will be considered in the next 
article [Avtomatika { Telemekhanika 21, 12 (1960), 
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Chasmar R, P, “An electrical multiplier utilizing 
the Hall Effect in indium arsenide,” Electronic 
Eng. 30, No, 369 (1958), 

Lofgren Lars, “Analog multiplier based on the 
Hall effect," Actes, Journees Internat, Calcul A- 
nalog, Brusels, 1955 (Brussels, 1956), [1I-151, Dis- 
cuss 115, 

Lofgren Lars, “Analog multiplier based on the 
Hall effect," J, of Appl. Phys, 29, No, 2 (1958). 
Hilsum C,, “Multiplication by semiconductors,” 
Electronic Eng. (November, 1958), 

Deeley E, M. and Mackey D. M., “Multiplication 
and division by electronic-analof methods," Na- 
ture, 163, No, 4147 (1949), 

Miller I, A., Solter A. S,, and Scott R, E,, "Wide- 
band analof function multiplier,” Electronics 28, 
2, 11 (1955). 

Schmidt W., “Die Hyperbeldrore, eine Electrones- 
trahlrohre zum multiplizeeren in Analogie-Rechen 
gerafen,” A. omgew. Phys., No, 2, (1956), 

Soltes A, S,, "Wide-band square-law circuit ele- 
ment," Trans, IRE ED-2, No, 2 (April, 1955), 
Angela E, I,, "An electron - beam tube for analog 
multiplier,” Rev. Sci, Instr, 25, No, 3 (March, 1954), 
Cundlach F, W., "A new electron-beam multiplier 
with an electrostatic hyperbolic field," Journees 
Internat, Calcul, Analog (Brussels, 1955), Actes 
1956, 
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LETTER TO THE EDITOR 
J. Maty45 


Translated from Avtomatika i Telemekhanika, Vol. 


pp. 1433-1434, October, 1960 


In its issue No, 1 of 1960, the journal Avtomatika 
i Telemekhanika published the article,"On methods for 
the simulation of rational- fraction transfer functions 
without differentiating elements? by B. Ya. Kogan [1]. 
Before the appearance of this article, no such surveys 
and comparisons of different methods for the solution of 
this problem were available in literature, Since cer- 
tain methods for the simulation of transfer functions 
have special advantages in various practical applica- 
tions, the article by B, Ya, Kogan must be recognized 
as being of interest to persons engaged in solving prac- 
tical problems by means of simulating devices, 

At the same time, the journal Slaboproud¢ Obzor 
No, 1, 1960, published my article," Methods for pro- 
gramming differential equations with constant coeffi- 
cients by means of simulating devices” [2], The above 
article treated also some other methods for the repro- 
duction of programming block diagrams by means of 
simulating devices, which were different from those 
proposed in[{1], I would like to give here a brief pre- 
sentation of these methods, thereby supplementing B. 
Ya, Kogan’s article, 

First, a few remarks concerning article [1]. 

1, The circuit in Fig. 1 [1] represents a circuit for 
dividing the input voltages by the coefficient a3, Con- 
sequently, this circuit is stable for all a,>0 (and not 
only for as>1, as was stated in[1]), For 0<a;<1, the 
resistor with the admittance a 3-1 is negative, but it 
can be connected to another feedback resistor, Thus, 

a resistor with admittance ay will be obtained in the 
feedback, For very small values of the a, coefficient, 
we shall obtain a circuit with a large amplification 
factor, which will be stable for all ag>0, (The above 
circuit is unstable only for ag=0,} For a3<0, the initial 
equation is multiplied by -1; for as=0, we shall obvi- 
ously obtain an equation of a lower order, The case 
where a = 1 is most often encountered in practice (gen- 
erally, a,,= 1). 

2, In the method of decomposing the initial dif- 
ferential equation into a system of first-order equations 
[1], the new coefficients a, (see Fig, 6a and the systems 
of equations (15), (16), and (16a) [1]) represent the 
initial conditions for the pulse transfer function 
w(t) = L'{ w(s)} of the system under consideration, 
which is given by 


Oty _p= wD (0), 
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21, No. 10, 


This statement has been proved in [2]. 


Composition of Block Diagrams for the 





Programming of Rational-Fraction Func- 





tions by Means of Simulating Devices 





without Auxiliary Computations 

Programming diagrams which contain all coe ffi- 
cients of the initial differential equation in nontrans- 
formed form have great advantages, In the above-men- 
tioned article by B. Ya, Kogan, only the direct integra- 
tion method and the method of combining the deriva- 
tives for b,= © have this property, The diagrams obtained 
as a result of using these methods are inverse (adjoint) 
with respect to each other, 

For m<n, other diagrams which have this property 
can be obtained by various methods, All such diagrams 
can be completely described by means of feedbacks and 
couplings “which supply a forward signal," Any such 
diagram can be composed in a purely formal manner 
without any auxiliary calculations by using the following 
rules and conditions: 

Rule 1, Any coefficient a, on the left-hand side of 
the initial differential equation corresponds in the dia- 
gram to a feedback loop through n-k integrators with 
the open-loop feedback gain -a,. 

Rule 2, Any coefficient bj, on the right-hand side 
of the initial differential equation corresponds in the 
programming diagram to a loop supplying a forward 
signal through k integrstors (passing from the system 
input to the system output through n-k integrators) with 
the gain b;. 

Condition 1, All loops supplying a forward signal 
must either go out from the system input (see Fig. 2[1)]) 
or converge to the system output (see Fig, 7[1)). 

Condition 2, Each integrator or summator in the 
diagram between the beginning of any feedback loop 
and the system output has only one input (no feedback 
loop or a loop supplying the forward signal must termi- 
nate in it), 

Every programming block diagram set up according 
to rules 1 and 2,and satisfying conditions 1 and 2, repre- 
sents a diagram which contains all coefficients of the 
initial differential equation in nontransformed form, 
and, conversely, any system which contains all the 
coefficients in nontransformed form can be composed in 
a purely formal manner according to the above rules 
and conditions, 





























As an example, the figure shows a programming 
block diagram for the case where n=5 and m=1, (The 
gains of all integrators and summators are assumed to be 
equal to -1.) 

Programming diagrams which contain all coetfi- 
cients of the initial differential equation in nontrans- 
formed form have the following important advantages; 


a) The programming block diagram can be obtained 
without any auxiliary calculations in a purely formal 
manner according to the above- formulated rules 1 and 2 
and conditions 1 and 2, 


b) The change in any coefficient of the initial equa- 
tion can be effected independently of other coefficients 
by means of a suitable potentiometer. 


c) Among these programming block diagrams, there 
are such diagrams which provide derivatives of the out- 


put quantity at the outputs of some of its integrators up 
to the n-m-1 order (see figure). 

It follows from the above that programming block 
diagrams which contain all coefficients of the initial 
differential equation in nontransformed form often have 
important advantages over other diagrams in solving 
practical problems. They are described in detail! in [2]. 


LITERATURE CITED 
1, 8B. Ya, Kogan, “On methods for the simulation of 
rational- fraction transfer functions without differen 
tiating elements," Avtomat, { Telernekh. 21, No. 1 
(1960).° 
’ a4 a4 4 
2. J. MatyaS, *Programovani linearnich diferencialnich 
rovnic s konstantnimi koeficienty na dif, analys4- 
torech,” Slaboproudy obzor, No. 1 (1960), 


* See English translation. 
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WORK OF THE ALL-MOSCOW SEMINAR ON MAGNETIC 


AUTOMATION COMPONENTS 
G. V. Subbotina 


Translated from Avtomatika i Telemekhanika, Vol, 21, No. 10, 


p. 1435, October, 1960 


During the first half of 1960, the All-Moscow 
Seminar on Magnetic Automation Components 
continued its work. The Seminar is being held at the 
Institute of Automation and Remote Control, Academy 
of Sciences, USSR, under the guidance of Dr. Tech. 
Sci. Prof. M. A, Rozenblat. 


At the Seminar session which was held on 
February 24th, M. A. Rozenblat and G, V. Subbotina 
submitted a report entitled ” Stability of multistage 
magnetic amplifiers,” The report dealt with a method 
for determining the transfer function of a multistage 
summing magnetic amplifier and an investigation of 
stability criteria of magnetic amplifiers with internal 
feedback (self-magnetization) by using frequency 
criteria and considering the lag in load circuits. The 
authors showed that the obtained stability criterion 
makes it possible to synthesize a stable two-stage 
magnetic amplifier with an assigned error for a 
certain given gain. 


In their report "Reversive dc magnetic amplifiers 
of greater efficiency? which was submitted on March 
23rd, M. A, Boyarchenkov and M, A. Rozenblat 
presented the general principles for the design of 
highly efficient reversive amplifiers for work with dc 
loads, A number of new transformerless reversive 
amplifier circuits have been proposed. 


On April 13th, N. P. Vasil’eva and I. S. 
Gashkovets submitted a report entitled “Stability of 
long circuits with logic elements,'* which dealt with 
the requirements for the “input — output” character- 
istics of logic elements for securing stable operation 
of annular or long circuits composed of identical logic- 
al repeaters or identical inverters,as well as circuits 
satisfying these requirements. 


On April 27th, R. A. Lipman and A. I. Moskalev 
reported on the results of investigations of magnetic 
amplifiers with self-magnetization and a dc output 
with a circuit for voltage doubling by means of cap- 
acitors. This circuit secures a twice greater amplifica- 
tion with respect to voltage in comparison with ‘the 
ordinarily used two-core circuits of this type. The 
circuit sensitivity to the feed voltage shape, which 
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must be sinusoidal, constitutes a deficiency of this 
circuit. 

In his report, "Magnetic transistor pulse amplifiers 
with de feed, “which was submitted on May 11th, L. N. 
Kiselev analyzed a circuit with both magnetic amplif- 
iers and transistors, operating under switching condit- 
ions, where the magnetic amplifier plays the role of 
the duty regulator. Since the magnetic amplifier rep- 
resents the input stage, the transistor invertors for feed- 
ing and controlling the magnetic amplifier must secure 
only a power sufficient for feeding the low- power mag- 
netic amplit..r and for controlling the output stage 
triodes. The characteristics of a 2-kw amplifier with 
an output stage based on four P- 207 transistors con- 
nected in parallel were given. 

On May 25th, M. A. Ol'shvang reported on the 
development of magnetic semiconductor amplifiers 
with ac outputs. A magnetic amplifier is used for the 
first stage, and a push-pull amplifier with transistors 
which operate under switching conditions is used as the 
second stage. The semiconductor stage gain is consid- 
erably increased by providing positive feedback with 
respect to the load current by means of a feedback 
transformer. Such amplifier circuits for the control of 
two- and three-phase asynchronous motors ,as well as a 
reversive circuit with a half-wave dc output,were 
described. 

The last session of the Seminar in the first half of 
this year was held on June 8th. At this session, Yu. A. 
Avakhom submitted the report, "Contro] of cores with 
a rectangular hysteresis loop.” The author proposed a 
new method of controlling cores for magnetic amplifi- 
ers with internal feedback with respect to the quasi- 
load characteristic which is close to the actual charac- 
teristic of magnetic amplifier control. The simulation 
of operating conditions for small and large control 
circuit resistances was secured by means of two differ- 
ent measuring circuits. 

The participants at the Seminar noted that the sub- 
jects of the reports are important and of interest from 
the theoretical point of view, as well as for practical 


application. The reports were recommended for publica- 
tion. 
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LIST OF LITERATURE ON MAGNETIC COMPONENTS 


FOR AUTOMATION, REMOTE CONTROL, 
AND COMPUTER TECHNOLOGY FOR 1959 


Compiled by G. V. Subbotina and I. 8. Trefilova 
Translated from Avtomatika i Telemekhanika, Vol, 21, No, 10, pp. 1436-1450, 


October, 1960 


GENERAL PROBLEMS 
(Terminology, bibliography, standardization 
problems, etc.) 

Measurement of magnetic parameters, Domestic 
and foreign book and periodical literature for 1955-1958; 
13 titles (No, 1767), Composite Bibliographic Index of 
Technical Papers, 1958 (Gos, Publich. Nauchno- Tekh, 
SO AN SSR, Moscow, 1959) p. 171, 

A. D. Nesterenko, On the Terminology in the Field 
of Magnetic Measurements, Problems in Magnetic 
Measurements ( Izd, AN UkrSSR, Kiev, 1959) p, 5. 

I, 1. Pekker, “All-Union seminar on contactless 
magnetic components for automation, remote control, 
and computer technology, Moscow, October, 1959", 
Izvest, Vyssh Ucheb, Zaved,, Elektromekhanika, 11, 
137 (1959), 

G. V. Subbotina, “Work of the All- Moscow seminar 
on magnetic automation components in", Avtomat.i 
Telemekh, 20, 12, 1698 (1959)? 

G. V. Subbotina, " List of literature on magnetic 
amplifiers and contactless magnetic components for 
19% (supplement)”, Avtomat. i Telemekh. 20, 3, 376 
(1959), x 

G. V. Subbotina and I, S, Trefilova, “ Magnetic 
components for automation, remote control, and com- 
puter technology". Annotated literature index for 1957 
(Sovetskoe Radio Izd., 1959) p. 71. 

G. V. Subbotina and I, S, Trefilova, ™ List of litera- 
ture on magnetic components for automation, remote 
control, and computer technology for 1958 (list of do- 
mestic papers)" Avtomat, i Telemekh, 20, 9, 1902 
(1959),° bad 

Yu, P, Filimonov and M, I, Kaplan, “ Literature on 
contactless magnetic components with a rectangular 
hysteresis loop", Radioelektr, Promy., 15, 55 (1959). 


FERROMAGNETIC MATERIALS AND CORES 
a) Magnetic Materials 
V. G. Borisenko, "On suitable methods of deter- 


mining the magnetic properties of cold-rolled trans formi- 


er steel”, Zavodskaya Lab, 25, 12, 1422 (1959)® Fig. 2, 
Table 3; Bibl, 4, 


V. G, Vasil'ev and V. A, Zverev, “Electronic 
simulation of hysteresis characteristics of ferromagnetic 
materials", Izvest, Vyssh, Ucheb, Zaved., fiektrome- 
khanika, 9, 3 (1959)* Bibl, 3. 

V. I, Prozhzhina and N, V, Erofeeva, “On the struc- 
ture of the family of symmetric hysteresis loops in ferro- 
magnetic materials", Izvest, Akad, Nauk SSSR, Sev. 
Fiz, 23, 3, 304 (1959)® Fig. 2, Bibl, 5, 

R. K, Kalinin’, “Ratio of magnetic losses in ferrite 
cores with an open-loop magnetic circuit," Tr, Inst, 
Fiz, AN LatvSSR, 11, 73 (1959) Bibl, 8, 

L. V. Kirenskii, et al,, "On technical magnetiza- 
tion of ferromagnetic materials,” Doklady Akad, Nauk 
SSSR 128, 2, 288 (1959) Bibl, 10, 

I, Kraus, “On the problem of securing a rectangular 
hysteresis loop in manganese ferrites,” Czechoslovak 
Physics Journal 9, 3, 408 (1959), 

N. I, Lapkin and V, S, Mes’kin, “On the effect of 


“high-temperature annealing on the magnetic properties 


of magnetically soft alloys," Izvest, Akad, Nauk SSR, 
Ser, Fiz, 23,.3, 280 (1959) Fig. 6, Bibl, 5. 

G. Matveev and V, Khomich, “Ferrites: new mag- 
netic materials," Radio, 8, 52 (1959), 

"On the determination of the properties of electrical 
steel” [Résuméd of the discussion in the journal Zavod- 
skaya Laboratoriya 25, 4 (1959)] Zavodskaya Lab, 25, 
12, 1415 (1959)* Bibl. 9. 

A, I, Pirogov, “Effect of temperature on the mag- 
netic reversal process in ferrite cores," Inzh, Fiz, Zhur, 
(Minsk, 1959), 

K, A, Piskarev, “An oxide magnetic material," 
Author's Certificate, Class 21g, 303, No, 123259, Byull, 
Izobre., 20 (1959), 

I. M, Puzei, “Temperature-stabilization of the 
magnetic properties of Permalloy," Doklady Akad, Nauk 
SSSR 120, 4, 768 (1958), 

L, I, Rabkin, "A magnetic core", Author's Certifi- 
cate, Class 21d", 48, 117963, Byull, Izobret, , No. 3 
(1959). 

L, I, Rabkin, “Certain problems in magnetic di- 
electrics physics,” Izvest, Akad, Nauk SSSR, Ser, Fiz., 
22, 10, 1276 (1958) Bibl, 6, 
































L. 1, Rabkin and Z, I, Novikova, “Electrical prop- 


erties of magnetic dielectrics and ferrites,” Izvest, Akad, 


Nauk SSR, Ser, Fiz, 23, 3, 388 (1959) Fig, 15, Bibl, 2, 
N, S. Siunov and V, M, Pavlinin, "Simultaneous 
magnetization of steel by currents of different frequency 
(50 and 200 cps), “Izvest, Vyssh. Ucheb, Zaved., Elec- 

tromekhanika, 3, 7 (1959) Fig, 12, Bibl, 4. 

L, V. Sobolev, "Voltampere characteristics and the 
temperature dependence of the conductivity of nickel— 
zinc ferrites, " Vest, Mosk, Ser, Unt,, Matem,, Mekhan,, 
Astron,, Fiz., Khim 1, 117 (1959) Bibl, 4, 

R. V. Telesnin, T. V. Dmitrieva, and S, A, Mate- 
vosyan, “On the temperature dependence of the mag- 
netic ductility of iron and some iron—nickel alloys," 
Nauchnye Doklady Vysshei Shkoy Fiz, i Matem, 5, 205 
(1958) Fig. 6, Bibl, 7. 

"Standard ferrite cores, " Radio, 8, 55 (1959), 

Ya, S, Shvartsbart, “Magnetically soft materials 
produced by the Elektrostal’ Plant," Collection: Mag- 
netic Amplifiers (Mat, Nauch,— Tekh, Goveshchariya, 
Khar’kov, 1958 Moscow, 1959) pp, 72-77, Table 4, 


b) Magnetization Processes and the Dy- 
namic Characteristics of Ferromagnetic 
Materials (Losses, Eddy Currents, Viscos- 
ity, etc.) 

V. V. Bardizh and V, V. Kobelev, “Calculation of 
magnetic reversal curves for ferrite cores," Izvest, ITM 
and VT AN SSSR (1959) p, 26, Fig. 15; Bibl. 3. 

V. F. Belyavskii and K, M, Polivanov, “On the dy- 
namic characteristics of ferromagnetic cores," Nauch, 
Doklady Vyssh, Shkol., Elektromekhanika i Avtomatika 
2, 14 (1959), Fig. 12, Bibl, 12, 

V. F. Belyavskii and K, M, Polivanov, “Cores made 
of anisotropic ferromagnetic strips,” Izvest, Vyssh. Ucheb. 
Zaved., Elektrotekhnika, 10, 3 (1959) Bibl, 1. 

V. Gaiko and I, Daniel—Sabo, “Effect of variable 
magnetic fields on the residual state of magnetically 
soft material,” Czechoslovak.Physics Journal 9, 1, 37 
(1959) Bibl, 5, 

E. I, Gurvich and E, I, Kondorskii, “Effect of micro- 
scopic inhomogeneities on the dynamic characteristics 
of magnetically soft alloys in weak fields," Izvest, Akad, 
Nauk SSSR, Ser, Fiz, 23, 3, 324 (1959) Fig. 3, Bibl, 12, 

V. L. Dyatlov, “Consideration of eddy currents and 
viscosity in the magnetic reversal of strip cores in large 
fields," Nauch, Doklady Vyssh, Shkol, Elektromekhanika 
i Avtomatika, 2, 3 (1959) Fig. 7, Bibl, 10, 


V. G, Elkin and A, I, Pirogov, “Dynamic character- 
istics of ferrites and their application in designing cir- 
cuits with ferrite cores, © Tr. VIRTU, Minsk, No, 15 
(1959). 


E, I. Klepfer and G, M. Tikhomirov, "Problems in 
the theory of ferrite cells," Izvest, Vyssh. Ucheb.Zaved., 
Elektromekhanika, No. 12 (1959). 
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A, 1, Pirogov,” Dynamics of pulse magnetic reversal 
in ferrites with a rectangular hysteresis loop,” Collec- 
tion; Tr. VTs AN SSSR (1959), 

A. I. Pirogov, Pulse Magnetic Reversal in Ferrites 
with a Rectangular Hysteresis Loop, Candidate's Thesis 
(Mosk, Energ. Inst., Moscow, 1959). 

M, A, Rozenblat, “Effect of eddy currents on the 
characteristics of magnetic amplifiers with feedback in 
the case of low-resistance control circuits,” Avtomat, 

i Telemekh, 20, 12, 1659 (1959)* Fig. 8, Bibl. 4. 

N, A. Smol’kov and V, F, Belov, “Certain properties 
of ferrites under pulse conditions," Izvest, Akad. Nauk 
SSSR, Ser, Fiz, 23, 3, 357 (1959) Fig. 5, Bibl, 15. 

R. V. Telesnin and E, F, Kuritsyn, "On the mag- 
netic reversal rate in ferrites," Izvest, Akad, Nauk SSSR, 
Ser. Fiz, 23, 3, 352 (1959) Fig. 6, Bibl. 9. 

R. V. Telesnin and A, G, Shishkov, “Effect of mag- 
netic viscosity on the frequency characteristics of fer- 
rites,” Izvest, Akad, Nauk SSSR Ser. Fiz, 23, 3, 343 
(1959) Fig. 8, Bibl, 8, 

Yu, M, Shamaev, V. L. Dyatlov, and A, I, Pirogov, 
“Dynamic characteristics of ferrites,” Nauch, Doklady 
Vyssh, Shkol, , Elektromekhanika No, 1 (1959), 


c) Core Design and Production Technol- 
ogy 

I. S, Babayants and A, I, Ivanov, “ A device for 
automatically maintaining a constant rotation speed of 
the bobbin in toroidal winding of cores,“ Author's Cer- 
tificate, Class 21g, 1m. No, 123627. Byull, Izobret., 
21 (1959). 

A. M, Bamdas, V, A, Somov, and A, O, Shmidt, 
"Certain design variants of single-phase and three-phase 
transformers which are controlled by shunt magnetiza- 
tion,” Tr. Gor'kovskogo Politekhnicheskogo 15, 4, 27 
(1959) Fig, 15, Bibl, 7, oe 

N. I. Vorob’ev, "The use of enamel as insulation 
material between sheets in magnetic cores made of 
electrical strip steel,” Izvest, Vysshikh Ucheb. Zaved., 
Elektromekhanika, 1, 125 (1959) Fig. 6, Bibl, 2, 

I, K, Levochkin, "The present state of development 
of a new line of magnetic amplifiers with magnetic 
circuits made of Il-shaped plates," Collection: Mag- 
netic Amplifiers (Scientificotechnical Conference Trans- 
actions, Khar’kov, 1958; Moscow, 1959) p, 40, Fig. 5. 

A. V. Netushil, “Effect of insulation thickness on 
the magnetic properties of foliated cores," Tr. MéI, 30, 
142 (1958), 

L, I. Rabkin, “A magnetic core," Author's Certif- 
icate, Class 21g, 3193, No, 120621, Byull, Izobre, 12 
(1959). 

B.G. Sokolov, “A line of 30- to 3000-w magnetic 
amplifiers with magnetic circuits consisting of three- 
leg plates," Collection: Magnetic Amplifiers (Scientific- 
otechnical Conference Transactions, Kharkov, 1958; 
Moscow, 1959) pp. 78-85, Fig. 4. 
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V. M. Kharlamov, “A method for the rejection of 
unsuitable magnetic materials,” Author's Certificate, 
Class 21e, 3749, No. 118902 (594056 of March 7, 1958). 
Byull. Izobret., 7 (1959). 

T. K. Tsogoev, “An adapter to the winding device 
for fastening the wire end at the lead-out of components 
(in particular, chokes)," Author's Certificate, Class 21g, 
1g, No. 123626. Byull. Izobret., 21 (1959). 

3. General Problems in the Theory of 
Nonlinear Magnetic Circuits (Calculation 
of Circuits with Steel Components, Ferro- 











resonance, etc.) 











L. A. Bessonov, “On the theory of oscillations in 
one of the ferroresonance circuits,” Izvest. Vysshikh 
Ucheb. Zaved., Elektromekhanika, 2, 43 (1959) Fig. 12, 
Bibl. 2. 

B. K. Bul’, “Investigation of the field in the vicinity 
of the air gap and the calculation of permeance,” Vest. 
Elektropromy. 9, 66 (1959). 

E F. Galteev and V. G. Morozov, “Compound 
ferroresonance voltage stabilizer for ac generators with 
permanent magnets,” Nauchnye Doklady Vysshei Shkoly, 
Elektromekhanika i Avtomatika, 2, 151 (1959) Fig. 

8, Bibl. 3. 

V. P. Glukhov, “Design of a circuit consisting of a 
resistor and a choke with steel,” Vestn. Elektropromy. , 
10, 42 (1959). 

B. Z. Zil*berman, “Condition for the maximum 
accuracy of ferroresonance voltage stabilizers," Radio- 
tekhnika 14, 4, 75 (1959). 

V. D. Krochakevich, Transient Processes in Three- 
Phase Circuits with Rectifier Bridges. Candidate's 
Thesis (Lenin Komsomol Riga Aircraft- Engineering 
Military Academy, 1959). 

M. S. Libkind and G. I. Kugushev, “Graphs for deter- 
mining the harmonic composition of magnetizing cur- 
rents in superimposing a steady field on a variable field,” 
Elektrichestvo, 7, 55 (1959) Fig. 5, Bibl. 1. 

E. M. Pevzner, “On the investigation of transient 
processes by means of dc currents in electric circuits 
with steel components,” Izvest. Vysshikh Ucheb. Zaved., 
Priborostroenie, 1, 63 (1959) Fig. 1, Bibl. 1. 

S. P. Pivovarov, “Design of ferroresonance voltage 
stabilizers," Vest. Elektroprom.,6, 60 (1959) Fig. 3, 
Bibl. 5. 

B. A. Piontkovskii, “A ferroresonance voltage stabi- 
lizer," Author's Certificate, Class 214", 35,,, No. 
121485. Byull. Izobret., 15 (1959). 

G. M. Tikhomirov, "Representation of the hysteresis 
loop equation in the form of an integral with a variable 
upper limit,” Izvest..Vysshikh Ucheb. Zaved., flektro- 
mekhanika, 6, 3 (1959). 


4. Magnetic Amplifiers, Theory, Cir- 

cuits, and Design 

a) Books, Monographs, and Dissertations 
M. A. Boyarchenkov and M. A. Rozenblat, Fast 


Reversive Electric Drives with Magnetic Amplifiers 
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